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ABSTRACT 

Aims. We present photometric properties and distance measurements of 252high redshift Type la supernovae (0.15 < z < 1.1) 
discovered during the first three years of the Supernova Legacy Survey (SNLS). These events were detected and their multi-colour 
light curves measured using the MegaPrime/MegaCam instrument at the Canada-France-Hawaii Telescope (CFHT), by repeatedly 
imaging four one-square degree fields in four bands. Follow-up spectroscopy was performed at the VLT, Gemini and Keck telescopes 
to confirm the nature of the supernovae and to measure their redshifts. 

Methods. Systematic uncertainties arising from light curve modeling are studied, making use of two techniques to derive the peak 
magnitude, shape and colour of the supernovae, and taking advantage of a precise calibration of the SNLS fields. 
Results. A flat ACDM cosmological fit to 231 SNLS high redshift Type la supernovae alone gives £l M = 0.211 ± 0.034(stat) ± 
0.069(sys). The dominant systematic uncertainty comes from uncertainties in the photometric calibration. Systematic uncertainties 
from light curve fitters come next with a total contribution of +0.026 on Qm- No clear evidence is found for a possible evolution of 
the slope {f>) of the colour-luminosity relation with redshift. 
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* Based on observations obtained with MegaPrime/MegaCam, a 
joint project of CFHT and CEA/DAPNIA, at the Canada-France- 
Hawaii Telescope (CFHT) which is operated by the National Research 
Council (NRC) of Canada, the Institut National des Sciences de 
l'Univers of the Centre National de la Recherche Scientifique (CNRS) 
of France, and the University of Hawaii. This work is based in part 
on data products produced at the Canadian Astronomy Data Centre 
as part of the Canada-France-Hawaii Telescope Legacy Survey, a col- 
laborative project of NRC and CNRS. Based on observations ob- 
tained at the European Southern Observatory using the Very Large 
Telescope on the Cerro Paranal (ESO Large Programme 171.A- 



0486 & 176.A-0589). Based on observations (programs GS-2003B-Q- 
8, GN-2003B-Q-9, GS-2004A-Q-11, GN-2004A-Q-19, GS-2004B-Q- 
31, GN-2004B-Q-16, GS-2005A-Q-1 1, GN-2005A-Q-1 1, GS-2005B- 
Q-6, GN-2005B-Q-7, GN-2006A-Q-7, GN-2006B-Q-10) obtained at 
the Gemini Observatory, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under a cooperative 
agreement with the NSF on behalf of the Gemini partnership: the 
National Science Foundation (United States), the Particle Physics and 
Astronomy Research Council (United Kingdom), the National Research 
Council (Canada), CONICYT (Chile), the Australian Research Council 
(Australia), CNPq (Brazil) and CONICET (Argentina). Based on ob- 
servations obtained at the W.M. Keck Observatory, which is operated 
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1. Introduction 



Since 1998 dRiess et al Jl998tlPerlmutter et al.ll999h . surveys of 
cosmologically distant Type la supernovae (SNe la) have shown 
that the expansion of the Universe is accelerating, distant SNe la 
being fainter than expected in a decelerating Universe. With the 
assumption that the Universe can be described on average as 
isotropic and homogeneous, this acceleration implies either the 
existence of a fluid with negative pressure, usually called dark 
energy, a cosmological constant, or modifications of gravity on 
cosmological scales. Several other cosmological probes such as 
cosmic microwave background (CMB) anisotropies and baryon 
acoustic oscillations (BAO) confirm this result, but SNe la obser- 
vations are currently the most sensitive technique to study dark 
energy or its alternatives, since they can be used to directly mea- 
sure the history of the expansion of the Universe. 

Recent results from high re dshift SNe la s urveys 
dAstier et all l2006t iRiess et al.l 12007b IWood-Vasev etail 120071: 
iFreedman et al.l l2009t) and combinat i ons of SNe la samples 
dKowalski et al.1 120081: iHicken et all l2009bl) give consistent 
measurements of the effective equation of state parameter 
of dark energy (w, the ratio of pressure over density) in the 
redshift range accessible with SNe data. Values consistent with 
w = — 1 are found, as expected for a cosmological constant, 
with uncertainties of order of 10%, including s ystematics. In 
contrast with those analyses, Kes sler et alj {2.009) report a much 
larger systematic uncertainty on w. They obtain a discrepancy 
Sw = 0.2 when using two different techniques to estimate SNe 
distances. We discuss this issue in 94.21 

This paper is the first of a set of three papers which present 
the cosmological analysis of the Supernova Legacy Survey 
(SNLS) after three years of operation. The SN la sample is 
about four t i mes a s large as the first year sample presented in 
lAstier etaTl d2006l hereafter A06). Since systematic uncertain- 
ties are close to the statistical limits, large efforts on all aspects 
of the analysis have been made to identify, reduce and propa- 
gate them to the final cosmological result. A key feature of the 
work presented in this set of papers is that the full analysis has 
been performed twice with significantly different techniques, in- 
cluding for photometry, detection of candidates, spectroscopic 
identification, calibration, and distance estimate. The step by 
step comparison of the results obtained with these different tech- 
niques made it possible to pin down systematics often ignored 
in previous works. This first paper describes the photometric 
data reduction and the estimation of light curve fit parameters 
for this sample. It also reports on cosmological constraints that 
can be obtained from SN LS supernovae alone. The second pa- 
per (IConlev et al.llin prepl hereafter C10) presents cosmological 
constraints obtained by combining the SNLS 3-yr sample with 
lower and higher redshift SN la data and the 3rd paper presents 
cosmological constra ints when combining S N data with other 
cosmological probes dSullivan et alJlin prepj) . 

After an overview of the survey in ^2) the photometric data 
reduction is described in ^together with the photometric tech- 
niques used. The t wo light curve fitti ng techniqu es used in 
this a nalysis, SiFTO dConlev et al.ll2008l) and SALT2 dGuv et alJ 
|2007|) . are presented in §|4] Their results are compared and com- 
bined in $5] to constrain Qm for a flat ACDM cosmological 
model. We conclude in $6] 
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the University of California and the National Aeronautics and Space 
Administration. The Observatory was made possible by the generous fi- 
nancial support of the W.M. Keck Foundation. Mark Sullivan acknowl- 
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2. Overview of the Supernova Legacy Survey 

The Supernova Legacy Survey (SNLS) uses data taken as part 
of the deep component of the five-year Canada-France-Hawaii 
Telescope Legacy Survey (CFHT-LS). CFHT-LS is an optical 
ima ging survey using th e one square degree MegaCam cam- 
era dBoulade et al.1 120031) on the CFHT The deep component 
conducted repeat imaging of 4 low Galactic-extinction fields 
(see A06 for the field coordinates). The survey ended in August 
20 08. A detailed d e script ion of the telescope and camera is given 
in iRegnault etai] d2009l) (hereafter R09). The data are time- 
sequenced with observations conducted every 3-4 nights in dark 
time. Four filters were used (denoted gMfMiMZ/ull similar to 
those of the Sloan Digital Sky Survey (SDSS). This data allows 
us to obtain high-quality multi-colour SN light curves. Since the 
first year analysis in A06 the cadence and exposure time of the 
Zm observations were increased to improve the signal-to-noise in 

the light curves of the most distant SNe. 

For the "real-time" SN searches (see IPerrett et alj lin prepj 
for details) the SNLS data wa s reduced by the CFHT-devel oped 
Elixir data reduction system dMagnier & Cuillandrell2004l) and 
processed through two independent searchpipelines, from which 
combined candidate lists were generatecfl Both pipelines gen- 
erated difference images, subtracting from the search images 
deep references constructed from previous observations, leaving 
only sources which had varied since the reference epoch. The 
subtracted images were searched using various automated tech- 
niques, with all likely candidates visually inspected by human 
eye. SN candidates were prioritised for spectr oscopic follow 
up us ing a photometric identification technique dSullivan et al.l 
120061) which identifies likely SNe la from 2-3 epochs of real- 
time photometry, predicting redshifts, phases, and temporal evo- 
lution in magnitude. New candidates were checked against 
a database of existing variable sources to exclude previously 
known AGN or variable stars. 

Spectroscopic follow-up was used to confirm SN types 
and measure redshifts. The survey was allocated 60 hours 
per semester on the Gemini North and South telescopes 
using the GMOS instruments dHook et all |2004|) and 120 
hours per year at the European Southern Observatory Very 
Large Telescopes (VLT ) using the FORS-1/2 instruments 
dAppenzeller et al.| [l998). Spectroscopic ti me was also obt ained 
at the Ke ck telescopes usin g the LRIS dOke et al.lll995h and 
DEIMOS dFaber et al.ll2003l) instruments (PI: Perlmutter). Some 
additional types and redshifts were obtained as part of a sep- 
arate detailed stu dies p rogram also undertaken using Keck- 
I/LRIS dEllis et al]l2008l) . 

Data from the first year o f the Gemini program were pre- 



sented in Howe 



inlBronder et al 



1 et al. | d2005l) . from the second and third years 



(2008) and data from May 2006 to July 2006 in 



IWalker et alJ ( in prep J) . Spe ctra obtained at VLT are published 
in lBalland et al.l 00091) (see lBaumont et alJl2008l for the reduc- 



tion technique used), and Keck s pectra from the detailed studies 
program are in Ell is et al ] d2008l) . All spectra were analysed and 
uniformly typed. Two sub-classes are considered in this work 
according to the confidence level of the spectroscopic identi- 
fication: certain SN la are denoted "SN la" (corresponding to 
th e confidence indice s CI 5 and 4 in the classification scheme 
of lHowell etai]|2005l) . and probable SN la (CI 3) are labeled 



Supplementary « M -band data was also acquired as part of CFHT- 
LS; these data were not time-sequenced and not used in the SN light 
curves. 

2 Candidate lists are available from 
http : / /legacy . astro . utoronto . ca/ 
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"SN la*". For most SNe, two identification techniques were 
used and their results cross- che cked (based on the a nalysis pre- 
sented in lHowell et al.ll2005l and lBalland et al.ll2009l) . 

3. Photometry measurement 

Elixir pre-processed images were retrieved from the Canadian 
Astronomy Data Centre (CADCfl This process performs ba- 
sic image "de-trending" (bias subtraction, flat-field correction 
and fringe removal in z'm and zm bands). The flat-fields and 
fringe maps are more precise than the ones used for the real- 
time SN search as they are constructed from an entire queue run 
of data (including non CFHT-LS data), from median stacks of 
twilight and science exposures for the flat-field and fringe maps 
respectively. The Elixir pipeline also attempts to produce images 
with photometric uniformity (i.e. a constant zero point across 
the mosaic) by constructing photometric correction frames from 
dithered observations of dense stellar fields. However, this pro- 
cess is not perfect and some radial trends in the photometric zero 
point remain. This is shown in detail in R09, who found a cen- 
tre to edge variation of 0.02 magnitude, and derived refined 
photometric correction frames. 

The subsequent treatment of images including sky back- 
ground subtraction, astrometry, and photometric correction has 
been performed in two independent pipelines. T he first one fol- 
lows A06. SExtractor dBertin & Arnoutsl[l996l) is used to pro- 
duce an image catalog and a sky background map that is sub- 
tracted from the image. Second order moments of the objects 
are derived (using an iterative Gaussian-weighted fit), and point- 
like sources are identified to estimate the image quality (here- 
after IQ, given by the FWHM of the point spread function). A 
weight map is derived using the spatial sky variance, the bad 
pixel map provided by Elixir, and a map of cosmic rays hits and 
satellite trails identified with dedicated algorithms. The astrom- 
etry of each CCD image is obtained using a match of the image 
catalog to an astro metric catalog derived from observations of 
IStone et ail d 1999 ) astrometric calibration regions and USNO-B 
dMonet et all 120031) or SDSS (lAdelman-McCarthv et ai]|2007t) 
star catalogs. The photometric corrections (more precisely the 
ratio of the ones derived in R09 to those applied to Elixir images) 
are not applied to the images but to the fluxes of objects in the 
catalogs. The other pipeline performs the same reduction steps 
but uses different software (for instance the Image Reduction and 
Analysis Facility packageQ). 



3.1. Measurement of the supernova fluxes 

Two photometry techniques are considered. We present both and 
compare their outcome. The goal of this study is to cross-check 
the methods, and evaluate photometric systematic uncertainties 
from the level of agreement reached in this comparison. This 
also allows us to select the best technique for deriving light curve 
parameters of the SNe. 



3.1 .1 . Simultaneous fit of galaxy and supernova fluxes 

This first method (hereafter method A) was described in A06. 
It consists first in resampling all the images of a field in each 



pass-band to the pixel grid of the best IQ image (hereafter called 
reference image) for which a PSF (point spread function) model 
is derived. Convolution kernels of the reference image to each 
aligned image are then derived for subsequent use. The photo- 
metric ratio of the two images is simply given by the integral of 
the kernel. The photometric fit then consists in fitting simultane- 
ously an image of the host galaxy (at the sampling and the im- 
age quality of the reference image) together with the position of 
the supernova and its fluxes in all images, convolving both the 
galaxy and PSF models using the kernels obtained previously. 
The supernova flux is forced to zero in images where its flux is 
negligible. When we calculate the pixel uncertainties, we delib- 
erately ignore the contribution of the supernova and host galaxy 
fluxes, so that PSF inaccuracies affect the flux of bright and faint 
sources in the same way (this is further discussed in 33.2l i. Note 
that the supernovae considered here are faint enough for the de- 
parture from statistical opt imality to be n egligible. 

Whereas DAOPHOT dStetsonl [l987h was used in A06 to 
model the PSF, we have since then developed an independent 
PSF modeling code for this analysis. Also, some improvements 
have been made in the selection of images for alignment. The 
reference frame is larger than the actual reference CCD image 
to cover gaps between CCDs in the MegaCam focal plane. For 
large dithers, images from other CCDs overlapping the refer- 
ence image are included so that different fractions of the same 
image can be aligned on different reference frames. This allows 
us to measure fluxes for all SNe in the field, including those that 
sometimes fall in the gaps between CCDs (depending on the ac- 
tual pointing). 

Some cuts are applied to the images: an IQ better than 1 .5" is 
required, the x 1 per degree of freedom of the convolution kernel 
fit has to be better than 1.5, and the image correlation coefficient 
at a 50 pixel lag has to be lower than 0.08. This latter cut allows 
us to get rid of images with large fringe residuals (in z'm and zm 
bands) or background subtraction problems. 

The ratio of images that pass the cuts, the number of expo- 
sures and total exposure time for each field and band for this 
dataset is shown in Table [TJ along with the average IQ of the 
selected images. The relatively low efficiencies of about 80% in 
z'm band for the fields D2 and D3 are due to issues with fringe 
subtractions for some MegaCam runs. 

For a known fixed position of the SN in the images, the pho- 
tometry fit is linear and we do not expect any bias for a perfect 
PSF model. However, when the position of the SN is fit simul- 
taneously with the fluxes, we have (for a Gaussian PSF, see § [B] 
for a proof) the following fractional bias in flux: 



Sf E[f]-f 
f f 



Varif) 



(1) 



3 http://www.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/cadc/ 

4 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 



where Var(f) is the variance of the light curve amplitude /, ob- 
tained from a fit combining observations at different epochs. 

Figure [TJ shows the expected bias of the observed SNe as 
a function of redshift in ^m^m^mZm bands. For the ym and z'm 
bands, the biases are respectively of 1 mmag and 0.6 mmag at 
z — 0.8, and 5 mmag and 2 mmag at z — 1. The gM and zm light 
curves have a low signal to noise ratio (S IN) at high redshift 
so we do not fit the position for those bands and use instead the 
weighted average of the ones obtained in the tm and z'm band fits. 
With an astrometric precision better than 10 milliarcsec for the 



5 The dimensionless image correlation coefficient at a lag A is given 
by 2 (ft ~ fi+&f I 2 jf, where /; is the flux in the pixel i (the average 
value of fj is close to zero as the images were previously background 
subtracted). 
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Table 1. Total exposure time and average IQ of the SNLS 3 year 
data set. The selected images are those that pass the quality cuts 
(see text for details). 
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transformation of coordinates from one image to another (esti- 
mated from the RMS of the residual match of the catalogs), and 
an average IQ of 0.9", we do not expect a bias due to the trans- 
formation of coordinates larger than 0.3 mmag (see Eq. IB. 21 ). 
As a consequence, this effect can be ignored, and the total bias is 
given by that obtained in iu which has the highest S /N (reaching 
a maximum value of 2 mmag at z = 1). Those biases are an order 
of magnitude smaller than calibration uncertainties for redshifts 
z < 0.8, and other sources of uncertainties at higher z (see 35.41 ). 
so we do not correct for them and ignore their negligible con- 
tribution to uncertainties in the following. The correlation intro- 
duced in this process between the magnitudes in different bands 
will also be neglected. 

As in A06, the photometry fit provides us with flux estimates 
for all exposures in a given night. Since we do not expect signif- 
icant variation of the SN luminosity on this time scale, we aver- 
age these measurements and use the additional scatter between 
measurements to account for the imperfect kernel evaluation and 
under-estimated statistical uncertainties due to pixel correlations 
(introduced by the resampling of images). A few outliers (about 
1% of the data) are rejected at this level, some of which may be 
due to unidentified cosmic ray hits. Since the measurements on 
the light curve are correlated since they share the uncertainty on 
the host galaxy flux and SN position (in rM and im bands), their 
covariance matrix is recorded for a subsequent usage in the light 
curve fits. The quality of the photometry fit and the fraction of 
outliers rejected at this level are the same as in A06. 

3.1.2. PSF photometry on image subtractions 

We now describe the second photometric method (method B). 
The PSF and IQ are measured for each image from a list of 
isolated stellar objects with 35-50 objects used for each CCD. 
Photometric alignment is performed using a multiplicative seal- 




Fig. 1. Statistical uncertainty on magnitudes (left panel) and ex- 
pected bias on PSF photometry with a simultaneous fit of the po- 
sition (right panel) as a function of supernova redshift. The dif- 
ferent gMfMiMZM MegaCam bands are shown respectively with 
green solid, red dashed, black dotted, and blue dashed-dotted 
curves. 



ing factor from a comparison to the set of tertiary calibrating 
stars, discussed in 33.21 Aperture photometry with a radius of 
4" is used for this alignment, and the average offset to PSF pho- 
tometry is recorded at the same time. 

SN flux measurement are performed on a CCD-by-CCD ba- 
sis; even though the telescope is dithered during observations 
on a given night, we measure the flux on each observation in- 
dependently. A series of deep reference images are constructed 
for each season by combining the data from every other season 
in bins of IQ. The calibrated images (and their weight-maps) 
entering a given reference frame are geometrically re-sampled 
to a common pixel coordinate system using a kernel conserv- 
ing flux. They are then combined using a weighted mean with 
5-cr outlier rejection, generating a deep reference and associated 
weight-map. Each field/filter/season combination typically has 
five or more statistically independent references each with a dif- 
ferent mean IQ and each at least 6 times deeper (or 36 times the 
integration) than the individual exposures from which the SN 
measurements are made. 

For each calibrated image containing SN light, the reference 
with the closest (but superior) IQ is geometrically re-sampled 
to the same pixel coordinate system, and the PSFs of the two 
images matched by degrading the reference image to match the 
image containing the SN light. This adjusted reference image 
is then subtracted from the science image. This process avoids 
any geometric re-sampling of SN pixels and introduces only a 
minimum correlated noise in the subtraction of the deep ref- 
erence. Flux measurement is then performed on the difference 
images using a custom-written PSF-fitting program using PSFs 
measured on the unsubtracted images, weighting the fit using the 
appropriate weight map. 

3.1.3. Flux measurement uncertainties 

The SNe data set considered consists of a sub-sample of all the 
SNe la detected and spectroscopically confirmed by the SNLS 
up to July 2006. They are listed Table[9] For those labeled "la*", 
one can not completely exclude Ib/c core collapse SNe. Light 
curves from 4 SNe out of 285 could not be obtained for the fol- 
lowing reasons: lack of observations (beginning or end of an 
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band Am e RMS for f ga , < 0.3 

g M -0.003 + 0.005 -0.046 + 0.007 0.028 
r M -0.006 ±0.004 -0.025 +0.003 0.023 
iu -0.009 ±0.004 -0.010 ±0.002 0.025 
z M -0.002 ±0.010 -0.028 ±0.005 0.054 

Table 2. Average difference of the fitted peak magnitudes ob- 
tained with methods A and B after the correction for the host 
galaxy over-subtraction (A), slope of the correlation with f ga i 
(e), and RMS for high S/N SNe in faint galaxies (see text for 
details). 



observing run, bad weather conditions, instrument failures), an- 
other variable object at the position of the SN (active host galaxy, 
or another SN in the same host), vicinity of a very bright star, or 
a specific location on the edge of the focal plane. 

The calibration of light curves on tertiary standard stars ob- 
tained with method A is detailed in 33.21 Light curves obtained 
with method B were calibrated in two steps, first with a 4"radius 
aperture photometry, then with an aperture to PSF photometry 
correction. We checked that both methods were equivalent, so 
that we can safely assume for this comparison that both light 
curve sets have identical calibration. 

In order to compare the flux scale of each individual light 
curve obtained with the two methods, both light curve sets were 
fitted with a fiducial light curve model and the resulting light 
curve amplitudes were compared. In this process, only the pho- 
tometric points present in both reductions were considered to 
minimise the statistical noise of the comparison (a residual sta- 
tistical noise is expected since different cuts for the selection of 
images within a night were applied and references for image sub- 
tractions vary from night to night for method B). 

From this comparison, a systematic over-subtraction of the 
galaxy was identified in the data set obtained with PSF photom- 
etry on image subtractions (method B; about 2% of the galaxy 
flux at the SN location, estimated with the PSF of the image 
with best IQ). This bias was traced to a systematic normalisa- 
tion offset of the convolution kernel: its integral was forced to 
unity after the current and reference images were photometri- 
cally aligned using 4" aperture photometry of field stars. This 
approach turned out to be incorrect. As the impact on SN fluxes 
is a linear function of the host galaxy flux at the location of the 
SN, we resorted to pursue the comparison, focusing on SNe in 
faint hosts, and correcting approximately for the identified bias. 
We account for this effect (for comparison purposes only) by 
fitting a linear relation between the difference of peak magni- 
tudes obtained with both pipelines and the galaxy to SN flux 
ratio {f g a i): m A - m B = Am + e x f gal . 

We exclude from this fit noisy light curves (S/N < 50 in 
gM^MiM and 20 in zm , where the S /N is given by the uncertainty 
on the fitted light curve amplitude) and assign a 20% uncertainty 
on fgai which is not well defined (in addition to the statistical 
uncertainty on the SN flux that enters in the ratio). The results 
are listed in Table [2] For §m^m^m bands, the residual dispersion 
of order of 0.025 mag can be attributed to imperfect PSF/kernel 
modeling in the pipelines. The larger dispersion in z^-band of 
0.05 mag is due to a much lower signal to noise ratio and fringe 
residuals in the images (note that we do not use exactly the same 
images in both pipelines so that the Poisson photon noise im- 
pacts the comparison). The average offsets from 0.002 or 0.009 
mag depending on the pass-band are marginally significant, the 
correction for the host galaxy over-subtraction of the order of 
0.01 magnitude on average for all bands being only approximate. 



band < f mm /f lma > 

g M -0.0007 ± 0.0008 

r M -0.0003 ± 0.0008 

i M -0.0014 ± 0.0012 
z M 0.0066 ± 0.0037 

Table 3. Ratio of measured flux before 22 restframe days (be- 
fore maximum) to the fitted lightcurve amplitude, averaged over 
events, for method A. 



Light curves obtained with the simultaneous fit (method 
A) are respectively 13%, 13%, 3% and 16% less scattered in 
gM^MiM and zm bands than those obtained with photometry on 
subtractions (method B) with an uncertainty of 2% for all bands. 
Those numbers are based on the standard deviation of residuals 
to a light curve fit per pass-band with a free amplitude, width and 
date of maximum light. About half of this difference can be at- 
tributed to the relative shallowness of the reference images used 
for the subtractions. 

In order to assess the accuracy of the host subtraction tech- 
nique of method A, we compute the average measured flux be- 
fore 22 rest frame days before maximum light (/ zero ) for each 
supernova and band, when available. Table |3]reports the average 
ratio of f zelo to the fitted light curve amplitude for method A, and 
does not display any significant bias. 

To summarise the comparison of photometric methods: 
Method B over-subtracts host galaxies by ~2%, which we later 
traced to the photometric alignment method prior to subtraction. 
Both methods agree on the amplitude of light curves (includ- 
ing the calibration transfer from field stars to supernovae) if one 
concentrates on supernovae hosted in faint galaxies. Method A 
is free from host subtraction biases at 1 to 3 mmag accuracy, 
depending on band. 

We therefore choose to use the light curve photometry from 
the simultaneous fit method (method A) in the subsequent anal- 
ysis. 

3.2. Photometric calibration uncertainties 

A catalog of calibrated tertiary stars has been obtained for each 
of the four CFHT- LS deep fields by R09 (those stars were cal- 
ibrated against the lLandoltl[l992l catalog of secondary UBVRI 
standard stars). 

It consists of a set of gMfM^MZM magnitudes for each star in a 
local natural system (average airmass of the observations, trans- 
mission function at the average position of the star in the focal 
plane). The photometric calibration consists in transfering the 
raw SN light curves to this system. For this purpose, a zero point 
{ZP = mag+2.5 log 10 (flux)) has to be associated with each light 
curve. It is simply obtained by the flux measurements of tertiary 
stars using the same photometry technique on the same set of 
images as the ones used to derive the SN light curve. 

The light curves obtained with the simultaneous photometric 
fit ( 93.1. U are calibrated as follows. The fluxes of tertiary stars 
on the same CCD as the SN are measured with the same PSF 
photometry technique as the one used for SNe, the only differ- 
ence being that a galaxy model is not considered. Least squares 
optimality dictates that the pixel uncertainties used for the fit in- 
clude the contribution of the star itself, but we do not do so: we 
do not want the relative weights of pixels to change between su- 
pernovae (faint objects) and tertiary stars (mostly bright objects) 
such that PSF innaccuracies affect both in the same way. So, the 
pixel weighting we adopt for bright stars is the one that is op- 
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Fig. 2. Residuals from the calibration fit ( PSF magnitude mi- 
nus aperture magnitude of the catalog) for CCD #13, band in 
field D 1 . Each point is a tertiary star. Those marked with open 
symbols were excluded from the fit based on cuts on RMS, mag- 
nitude, and number of observations. The fit is performed itera- 
tively with a 2.5 <x outlier rejection. 



timal for faint objects, and this choice does not adversly affect 
the statistical uncertainty of the flux ratio of SN to tertiary stars, 
since the later are brighter and more numerous. 

Since we obtain a flux measurement for each observation, 
this allows us to check for the photometric alignment, and dis- 
card variable stars. Before assigning a zero point to the SN us- 
ing the average fluxes of stars and their associated magnitudes, 
the photometric corrections for the reference image, provided by 
R09, are applied to the stars and SN measurements, depending 
on their location on the focal plane. At this level, the zero point 
uniformity within each reference CCD frame is checked. 

As an example, residuals from the calibration of a central 
CCD in field D 1 for r M band are shown in Fig.|2] The RMS of the 
differences between the PSF magnitudes and the aperture mag- 
nitudes (used in R09) for the tertiary stars are shown in Fig. [3] 
they are typically of 0.008 mag. With about 30 stars per CCD on 
average, the typical statistical uncertainty on the zero point de- 
termination of a light curve is of order of 0.002 mag. Since this 
number is averaged out when combining SNe observed on dif- 
ferent fields on different CCDs, it adds a negligible uncertainty 
to the cosmological analysis. 

Chromatic systematic effects were identified in the residual 
differences between PSF and aperture magnitude as shown in 
Fig. 3] These are expected because of the PSF variation with 
wavelength which induces a colour term in the comparison of 
PSF and aperture magnitudes (for an achromatic PSF model per 
pass-band). At level of a few millimagnitudes, one cannot ex- 
clude some chromatic effects in the aperture photometry, which 
is contaminated by light reflected in the MegaPrime optical sys- 
tem, forming haloes around stars in the images. 

In the rMiMZM bands, the colour dependent offset is of or- 
der of 0.002 mag. For such a small effect, we cannot separate 
it from the residual offset as a function of magnitude (see be- 
low) since there is a correlation between the magnitudes and 
colours of the tertiary stars. We will hence treat this offset as 
an additional source of systematic uncertainty, rather than try 
and correct for it. In gM band, the effect is however much larger. 
It is clearly related to the colours of stars. It induces a relative 
change of calibration of 0.015 magnitude for SNe at redshifts 
0.2 and 0.6 (at higher redshifts, g^-band light curves are not 
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Fig. 3. RMS of the differences between PSF and aperture mag- 
nitude of the tertiary stars for the ^m^m^mZm bands. In each 
histogram, there are 36 x 4 entries, each corresponding to a 
CCD/field combination, for which a zero point is determined. 



used to derive distances, since they correspond to observations 
at rest-frame wavelengths shorter than 3000 A, see ® . 

In order to account for this PSF variation with wavelength, 
we have modified the effective transmission of the instrument 
in gM-band with a multiplicative correction depending linearly 
with wavelength: corr(A) = 1 + 0.048(/* - 4979 A)/(1000 A). 
This correction was adjusted on the residuals of Fig. [4] us- 
ing synthetic ma gnitudes of stellar spec tr a obtained with the 
Phoenix model faauschi dt et ail Il997t iBaron & Hauschildd 
119981: lHauschildt & Baronlll999l and references thereinfl The 
resulting colour dependent effect is shown as a blue curve in 
Fig-EKforgAf-band only). 

As a last check of the calibration to tertiary stars, Figure [5] 
presents the average differences between PSF and aperture mag- 
nitudes as a function of the magnitudes of the tertiary stars, 
where PSF magnitudes were corrected for the colour dependent 
terms. The brightest stars which are saturated on a fraction of 
the images (depending on image IQ, exposure time, sky trans- 
parency and Poisson fluctuation) were excluded from this anal- 
ysis because of the potential biases introduced by the selection 
of unsaturated observations (excluding positive statistical fluctu- 
ations reaching the saturation level results in a bias on the aver- 
age). The dimmest stars were also discarded to avoid biases due 
to the preferential detection of positive Poisson fluctuations of 
the signal. The range of magnitude considered was obtained by 
selecting stars on the plateau of Fig. 16 in R09 which presents 
the number of observations of the tertiary stars as a function of 
magnitude. Discrepancies for the dimmest stars are visible: those 
are consistent with a residual background in the aperture of the 
tertiary stars. In R09 §4.2, a residual background of order of 
+0.06, -0.03, -0.23, and -0.04 Analog Digital Units (hereafter 
ADU) per pixel is found for the §m^m^mZm bands. Assigning a 
systematic uncertainty of 0.1 ADU per pixel on this background 
for gMfMZM, and 0.2 ADU for the /^-band, leads to the shaded 
error bars shown in Fig. [5] Accounting for this, we assign a 0.002 
mag systematic uncertainty on the photometric calibration trans- 
fer. 



6 The study presented here relies on version 2.6.1 of the Phoenix / 
GAIA spectral library, that can be retrieved from the Phoenix ftp server: 
ftp : // ftp . hs . uni-hamburg . de/pub/outgoing/phoenix/GAIA 
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Fig. 4. Differences between PSF and aperture magnitudes of the 
tertiary stars in the gMfMiMZM bands as a function of the colour 
of stars (gray dots). The black points with error bars represent the 
average deviation and its uncertainty in bins of colour. Typical 
SNe colours at maximum light are marked with red dotted ver- 
tical lines. The blue curve on the top panel shows the effect of 
the PSF wavelength dependent correction on synthetic magni- 
tudes obtained with PHOENIX stellar models. The colours of 
BD+17 4708 are marked by the black vertical lines. 



In this section, we have presented the photometry and cali- 
bration of the SNLS third year SN sample. One of the two tech- 



niques developed in the collaboration has been selected based 
on a comparison of their performances. We use this light curve 
data set in 35. 1 I to determine for each supernova the parameters 
needed to estimate its distance. 



3.3. Interpretation of the photometric calibration 

The calibrated data are formally the ratio of the target flux F$n 
as measured by the detector to that of a reference star F re f that 
would have been observed in the exact same conditions. In prac- 
tice, we have determined a flux F meas in ADU and a zero point 
ZPmeas which provides a normalisation coefficient. 

As we will see in the next section, the supernova light curves 
are fit with a spectral sequence model Msn{^, t, ...) in order to 
extract from the light curves a flux intensity and additional pa- 
rameters that characterise the diversity of SNe la. For this pur- 
pose, one has to be equipped with a model of the instrument 
response as a function of wavelength T(A) for each pass-band in 
order to compare the model to the observations. 

The light curve fit consists in comparing the measured quan- 
tity 



x ^Q-0MZP mm -m ref ) 



to the model 



D 



model 



FSN/Fref 



J T(A)M SN (A,t,...)dA 
jT(A)M ref (A)dA 



(2) 



(3) 



where M re f{A) is the spectrum of the reference star and m re f its 
magnitude. The integrals in Eqj3]have the dimension of counts 
in the detector per unit time (Le. in units of ADUs -1 ), as they 
model the actual observations^. However, Eq. [3] tells us we do 
not need to know the normalisation of the instrumental response 
function, even from pass-band to pass-band, as it cancels in the 
ratio. 

For SNLS, the magnitude system defined in R09 has been 
anchored to the star BD +17 4708 with a prescription for its mag- 
nitudes in this system. As we have photometrically aligned the 
SN light curves to the catalog of stars defining this system, we 
use BD +17 4708 as our reference star and consider for its mag- 
nitudes those provided by R09: g M = 9.6906, r M = 9.2183, 
im = 8.9142, and zm = 8.7736. Note that those magnitudes are 
conventional, one could have changed them in the definition of 
the magnitude system as long as the magnitudes of the tertiary 
stars and hence the values of ZP meas in Eq. [2] are modified ac- 
cordingly. We use for M re f(A) the spectrum of B D +17 4708 as 
measured on HST by iBohlin & Gillilandl d2004 : we retrieved 
the latest version of this spectrum on the CALSPEC database^. 

Concerning the instrumental model (T(A)), a complication 
arises from the fact that the transmission of MegaCam filters 
varies across the focal plane, following a radial pattern, with 
a typical variation of the central wavelength up to 5 nm from 
the centre to the edge. The magnitude system has been designed 
so that no colour transformations are required to translate our 
SN measurements to this system; the small residual airmass- 
dependent colour corrections cancel out on average and do not 
contribute significantly to the uncertainty of the measurements, 



7 As a consequence, if and M re f are spectral energy densities 
(in ergcirT 2 s' 1 A -1 for instance), then T(X) is the product of an effec- 
tive collection area (accounting for all the transmission and reflection 
efficiencies in the optical path), a dimensionless quantum efficiency, the 
CCD gain, and the inverse of the photon energy (in Aerg~'). 

8 ftp : //ftp . stsci . edu/cdbs/current_calspec/ 
bd_17d4708_stisnic_082 . ascii 
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Fig. 5. Average differences between PSF and aperture magnitude 
of the tertiary stars in §m^m^mZm bands as a function of the star 
magnitude (once corrected for the colour dependent terms). The 
back error bars represent the statistical uncertainty on the av- 
erage. The shaded areas represent the uncertainty on aperture 
magnitudes due to a systematic uncertainty on the residual back- 
ground in the images used for calibration (0. 1 ADU per pixel for 
gMfMZM, 0.2 ADU for the / M -band) 



and the variation of the effective pass-bands within the dithering 
pattern of the DEEP field observations are negligible. But the 
magnitude system has to be interpreted using a model of the re- 



sponse function of the instrument that follows this radial pattern. 
We use the transmissions provided by R09 to integrate the spec- 
tral model of supernovae, using the average position of the SNe 
observations on the focal plane. Note also that the magnitudes 
of BD + 17 4708 have been chosen to be the same everywhere in 
the MegaPrime focal plane despite the fact that the filter trans- 
missions vary. 

The SNe la light curves of the SNLS 3-year sample are 
available on-line at the Centre de Donnees astronomiques de 
Strasbourg (CDS). An example of the published data is given 
in Table [TO] 



4. Modeling the supernova light curves 

4.1. General considerations 

The only two light curve parameters that have been found so far 
to correlate with the luminosity of SNe la are the width of their 
light curves (first measured in Johnson B-band, see e.g. lPhillipi 
1 19931 and then extended to other bands) and their colour as mea- 
sured, e.g., by the difference of magnitudes (or ratio of fluxes) in 
rest-frame B and V bands. Once corrected for those, the absolute 
maximum luminosity has a dispersion of the order of 15%. All 
recent high-redshift cosmological analyses use those parameters 
(luminosity, light curve sha pe and colour) in a more or less obvi- 
ous way to derive distances dRiess et al. 2004, 20 071 lAstier et al.l 
120061: IWood-Vasev et al.l 120071) with the notable exception of 
the CMAGIC tech nique which relies on a colour-magnitu de di- 
agram of SNe la (I Wang et alJl2003bt IConlev et al.ll2006h . The 
challenge of a cosmological application is to derive those pa- 
rameters with a minimal redshift dependent bias. 



4.1.1. An empirical modeling 

The goal of the light curve fitting is then to evaluate for each 
SN those parameters from observations performed with a lim- 
ited set of observer-frame filters and a limited cadence of ob- 
servations. This requires a model of the spectral sequence of 
the SN in order to interpolate among observations. Despite the 
fact that there is a broad consensus as to the basic physical pic- 
ture of the explosion mechanism, it is extremely difficult cur- 
rently to make quantitative predictions for the observed signal 
based on a physical model. Indeed extremely precise 3D mod- 
eling is required in order to simulate the flame propagation in 
the SN progenitor. As a consequence, an empirical modeling 
of the observables is needed. Historically, light curve templates 
were built in a limite d set of filters from a sample of nearby 
SNe (see for instance iGoldhaber et al.ll200ll) . This required a 
correction of the obs ervations for redshi fted supernovae, usually 
called ^-correction dNugent et al.ll2002h . Those were performed 
using an average spectral sequence based on a set of spectra 
obtained at different phases (days after maximum light) of the 
SN. This method i s applied for the MLCS2k2 light curve fit- 
ter dJha et al.ll2007l) . with tabulated /^-corrections as a function 
of phase, redshift and colour. More recently, techniques based 
on an explicit modeling of the spectral sequence have been de- 
veloped. The data are not corrected to rest-frame filters but di- 
rectly compared to the integral of the spectra in a mode l of the in- 
strumental re sponse (SALT(2):lGuv et al.ll2005l I20071 hereafter 
G07, SiFTO: IConlev et al.ll2008l hereafter C08). The advantage 
of this approach is to keep track of the correlations between the 
light curve shape, colours and the spectral properties in the fit- 
ting process. 
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4.1 .2. Impact of a limited training sample: using high-z SNe 

This light curve fitting technique is a fundamental ingredient of 
the cosmological analysis. Especially, the assumed broad-band 
colour relations (i.e. the relative amplitude of the SN spectral 
model at a wavelength scale of order of 1000A, beyond a sim- 
ple colour tilt, in other words the curvature of the spectra) in the 
wavelength range of validity of the model have a direct impact 
on the derived distances. In order to illustrate this, let us con- 
sider two SNe observed in and im band at redshifts of 0.5 
and 0.8. Since those tm and Im observations correspond to the 
rest-frame B and V, and U and B bands respectively, the ratio 
of distances derived for those two SNe is directly a function of 
the (U - B) - (B - V) colouiQ difference of the model. Since all 
light curve models are empirically derived from a limited train- 
ing set, this latter colour has an uncertainty which introduces a 
redshift-dependent cor relation among the derived supernova dis- 
tances (see for instance lKnop et alj2003l) . Since many more SNe 
have been observed at high-redshift than at low redshift, high-z 
SNe must be considered in the training of the light curve models 
in order to overcome the statistical limitations of the nearby sam- 
ple. This has been done with the SALT2 and SiFTO models; it 
was possible since both techniques do not make use of distances 
in their training process. 

4.1 .3. Modeling of the near UV emission 

High-z SNe allow the observation of the rest-frame near UV 
emission from the ground without the need of space telescopes. 
The near UV is modelled in SALT2 and SiFTO using SNLS pho- 
tometric (in gM and bands up to a redshift of 1) and spectro- 
scopic observations (see references in Table [5). Using near UV 
data allows for a drastic improvement of the colour and hence 
distance estimate for SNe at redshifts of order of unity, where the 
sensitivity of the rest-frame B and V is limited by the quantum 
efficiency drop of MegaCam CCDs in the zm band. However, we 
still lack spectroscopic observations at early and late phases (the 
primary goal of the SNLS spectroscopic program was to provide 
an identification of the SNe which is easier at maximum bright- 
ness). 

4.1 .4. Diversity of SNe la colours: intrinsic variation or 
absorption by dust 

There is still much debate about the treatment of the SN colour 
parameter (generally anchored to B — V at maximum light). 
Whereas all cosmological analyses based on SNe perform a lin- 
ear correction of distance moduli (i.e. logarithm of distances) 
with the measured colour, the value of the coefficient used and 
its interpretation differ significantly from one analysis to another. 
In A06, this coefficient B is marginalised over in the cosmolog- 
ical fit, without any attempt to separate the reddening effect of 
dust absorption or a potential intrinsic variati on. On the contrary, 
the M LCS2k2 tec hnique used in ESSENC E (IWood-Vasev et al] 
l2007h. GOODS dRiess et alJ [20041 [2007h and SDSS surveys 
(IKessler et al 1 120091) . assumes that the derived (B — V) colour 
offset primarily comes from extinction by dust, and therefore 
that the B parameter should be identified with the Rb value of the 



ICardelli et al.l (Il989l) extinction law. When B is fit at the same 
time as cosmology, values ranging from about 2 to 3 are found 
depending on the technique used to derive it. Th ose values are 
system atically smaller than the value of 4.1 in the ICardelli et alJ 
(1 1989b extinction law. 

The large range of values obtained for B is likely to be a 
consequence of different assumptions on the uncertainties of the 
(B - V) colour estimates (and to a lesser extent the colour range 
of the SNe sample considered). For a given data set, the larger 
the assumed uncertain ties on (B—V), the large r the fitted B value. 
This issue is raised bv lFreedman e t al. (2003); we come back to 
it in 35.51 Whereas in previous papers (including papers from 
the SNLS collaboration: A06, G07, C08), low values of B = 2 
were found, we get larger values when accounting for an intrin- 
sic scatter in SNe colour relations in this paper (see 34.41 and 



9 The SNe colours mentioned in this section are always considered at 
maximum light. Note however that the estimates of the colours resulting 
from a light curve fit are actually an average of the difference between 
the data and the model with a weight that varies with the phase of the 
observations. 



35.5l i. Fitting for B o r not has some consequences. For instance, 
IConlev et al.l (120071) have shown that either we live at the cen- 
tre of an under -dense region of the Universe as proposed by 
iJha et al 1 (120071) . or the relation between SN colours and lumi- 
nosity does not follow the one expected for the Galactic extinc- 
tion and B < Rb- 

This low value of B points to either an unusual extinction law 
in host galaxies of SNe la or an intrinsic colour variation that 
dominates the effect of extinction. One hint is that the colour 
variation law (which describes how the SN flux varies with 
colour as a function of wavelength) can be derive d from the SN 
data th emselves, and differs significantly from the ICardelli et al.l 
(1 19891) extinction law in the near UV and U-ba nd, even for ex- 
treme values of R B (see iGuv et al.ll2005l 120071 and Fig. [6). In 
SiFTO, the derived relation between the (U - B) and (B - V) 
colours of SNe can not be explained with an extinction law ei- 
ther. While there is not yet a definitive proof that the colour vari- 
ation we observe is intrinsic to the SN, we still have to relax the 
ass umption that i t is pur ely due to dust extinction as modelled 
by ICardelli et al.l (Il989l) . This has some consequences for the 
cosmological analysis. Indeed, applying an incorrect correction 
to luminosity introduces a redshift dependent bias since the av- 
erage colour of SNe varies with redshift because of Malmquist 
bias (bluer SNe are brighter and hence dominant near the de- 
tection limit of a survey). This occurs at the highest redshifts of 
all surveys but also for nearby SNe that were observed by other 
means (see e.g. lConlev et a l. 2007). 

In the MLCS2k2 approach, a colour excess E{B - V) is mea- 
sured as the difference between the observed colour and that of 
the model. In this model, the intrinsic variability of SNe is ad- 
dressed with a single parameter (A), and any possible additional 
intrinsic variation is unaccounted for. As a consequence, one ex- 
pects that, at some level, the E(B - V) estimate resulting from 
the light-curve fit combines both dust reddening and a possible 
residual intrinsic colour variation. There are several examples of 
SNe la being clearly extinguished by dust. In contrast, there is 
no proof that part of the (B - V) colour variation is driven by 
intrinsic SNe properties. The only hints come from observations 
which point to a value of B < Rb and a colour variation law 
incompatible with a standard dust extinction law. Nevertheless, 
as the physical mechanism responsible for the SN la explosions 
obviously involves more than a single parameter (composition 
of the progenitor, ignition conditions, ...), it is very possible that 
some properties of the SN impact on the observed (B - V) colour 
beyond that addres sed by A. 

As an example. iKasen et"aT1 (120091) present a class of super- 
novae simulations with various ignition conditions, deflagration 
to detonation transition (DDT), and viewing angles, which ex- 
hibit intrinsic variations of colour. The models with the same 
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DDT criterion present a range of (B - V) colours (at maximum 
light) that are not fully correlated with the decline ratJ^l Such an 
effect, which is present in those simulations, has not been ruled 
out on real data. 

It is further possible that the slope of the colour-luminosity 
relation evolves as a function of redshift. A larger range of dust 
opacities is expected at higher redshifts where galaxies host 
younger stellar populations, so that the relative weight of dust 
extinction and intrinsic variation in the colour - lumin osity rela- 
tion should evolve with redshift. Kessleret al. (2009) found an 
evolution of the B parameter with redshift using S ALT2, but with 
an opposite trend with respect to the expectations (they found a 
lower B at higher z). We study this issue in 95.51 

4.2. Comparison of light curve fitters and distance estimate 

Using different light curve fitters to estimate distances allows us 
to quantify the systematic uncertainties associated with this step 
of the analysis. As shown in CIO and to a lesser extent in jj5.6l of 
this paper, the combined uncertainties on cosmological param- 
eters are dominated by systematics, so the statistical precision 
of light curve fitters do not play a major role in their selection. 
For a cosmology application, the most crucial aspects are those 
that may lead to a redshift-dependent bias: i) /^-corrections un- 
certainties, ii) any bias that could arise when fitting low signal to 
noise light curves, iii) biases associated with selection effects. 

Several authors have already performed this comparison of 
light curve fit t ers. T he largest systematic effects were found by 
iKessler et al.1 d2009t) (hereafter K09), who applied MLCS2k2 
and SALT2 to a sample combining low redshift SNe, HST data 
and the first release of SDSS-II, ESSENCE and SNLS data. They 
found a difference on w of 0.2 when using MLCS2k2 or SALT2. 
This difference is larger that any other source of uncertainty. The 
authors have looked for the sources of the observed discrepancy 
(see § 1 1 in K09 for a detailed discussion). For this purpose, light 
curves were fit with modified versions of MLCS2k2 designed to 
replicate some elements of the SALT2 model. They explain the 
discrepancies by two main differences between the light curve 
fitters: differences in the light curve templates, and the use of 
colour priors. We now review these two effects, and discuss af- 
terwards the differences in the treatment of colour. We eventually 
summarise this comparison section. 

4.2.1 . Systematic differences in the light curve templates 

IKessler et al.l (I2009I) provide evidence that SALT2 and 
MLCS2k2 predict differently the rest-frame [/-band. Excluding 
rest-frame [/-band light curves from the fit of SDSS SNe at 
z > 0.21 changes the MLCS2k2 and SALT2 distance moduli by 
^ 0.13 and 0.07 magnitude respectivel\FT. Note that the offset 
found for SALT2 is marginally consistent (at the 2<x level) with 
the calibration systematic uncertainties of order of 0.01 magni- 
tude for all bands reported both for the SDSS-II and SNLS first 
year releases (the SALT2 version used in K09 was trained on 
SNLS first year data). 

The SALT2 and SiFTO light curve fitters are more reliable 
than MLCS2k2 in rest-frame [/-band because they were cali- 



10 We selected models with the DDT criterion 3 in Table 2 of 
iKasen et"aT1 d2009h (supplementary material) and found a residual (B - 
V) dispersion of 0.03 after correcting for a marginal correlation with de- 
cline rate. The model spectra were kindly provided to us by the author. 

11 We consider here the relative change of distance moduli with re- 
spect to SNe at z < 0.21 



brated using high z supernovae. They benefit from the SNLS 
precise calibration and are not sensitive to the systematic un- 
certainties that are known to affect the observer-frame t/-band 
calibration of low redshift SNe. In particular, the effective U- 
band response function is poorly known for many SNe observa- 
tions because of important variations of the atmospheric trans- 
mission at wavelengths shorter than 350 nm. Also, when fitting 
with SALT2 or SiFTO (see S3 and COl l. we find a larger dis- 
persion of residuals in the U -band for low redshift SNe (0. 1 mag, 
RMS) than in the gM band at z — 0.4 where gM roughly corre- 
sponds to rest-frame U (0.05 mag). The systematic offset found 
by K09 could possibly be affected by some evolution of SNe 
properties in the UV. There is no evidence for such an effect, but 
even in this case, using SNe at z — 0.4 to calibrate the model rest- 
frame U -band is more reliable than using SNe at lower z. Indeed, 
the model rest-frame U -band is primarily used to estimate dis- 
tances of SNe at high redshift, and since a monotonous change 
of the average SNe UV emission with redshift is expected due 
to an evolution of metallicity or age of the progenitor popula- 
tion, U -band measurements done at the highest possible redshift 
should therefore be used. 

4.2.2. Using priors on SNe colours. 

In the MLCS2k2 approach, a prior is used that forces the extinc- 
tion Ay to be positive. This is equivalent to applying a prior on 
colour excess. However, since there is no evidence that the ob- 
served variation of SNe la colours is entirely and solely due to 
extinction by dust, applying such a prior is not justified. 

Furthermore, even if the model was qualitatively correct, ap- 
plying a prior is non optimal since any bias on this prior will pro- 
duce a bias on cosmology and hence artificially large systematic 
uncertainties. For insta nce, more than half of the difference on w 
between the analysis of IWood-Vasev et al]|2007l and that of K09 
is due to the choice of prior (see §10.1.4 in K09). Alternatively, 
not using a prior on colours does not bias the estimate; one only 
obtains larger and reliable statistical uncertainties. 

4.2.3. Estimating distances 

SALT2 and SiFTO use the colour in their distance modulus 
whereas MLCS2k2 uses a colour excess. This might be regarded 
as conceptually different, but, as shown below, the distance mod- 
uli used in the two approaches are mathematically equivalent, as 
long as the slopes of the brighter-slower and colour-luminosity 
relations are fit to the Hubble diagram. 

MLCS2k2 defines the colour excess as E(B - V) = (B - 
V) bs - (B - V) mo dei- At first order in A, the shape parameter, 
(B - V) mo dei — a + b x A and the distance modulus can then be 
written (at first order in shape and colour parameters): 

Umlcs = m B -M B + aA-R B E(B-V) 

= m B - M B + (a + R B b)A-R B (B - V) obs +R B a 
— m B — M' B + a' shape — BC 

where shape is any shape parameter and C — (B - V) /„ or any 
affine function of it. This last expression is exactly the one used 
in A06 and in this paper (Eq. [5J. As a consequence, choosing to 
use a colour excess or directly the observed colour to parametrise 
the distance modulus only modifies the meaning and values of 
M B and a coefficients, and does not change distances. 

So at first order, when no prior is used on colour and R B is 
treated as a free parameter (B), the MLCS2k2 approach is con- 
tained in the SALT2/SiFTO as both the colour variation law 
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(or colour relations, see 94.3.2b and the slope of the colour- 
luminosity relation of SALT2/SiFTO are parametrised in a way 
that allows the model to match a standard dust extinction law. 
Disentangling the contribution of intrinsic colour variation and 
dust extinction requires adding parameters beyond the simple 
one-parameter description of the SN intrinsic variation imple- 
mented in MLCS2k2, SiFTO and SALT2 (A, stretch, Xi). As a 
consequence, when using these fitters, the observed colour lu- 
minosity relation could evolve with redshift as it combines two 
effects whose relative weight is expected to change with redshift. 
We investigate this in 95.51 

4.2.4. Summary of the light curve fitter comparison 

We have shown above that using a colour or colour excess in the 
distance estimator is not the source of the differences between 
SALT2/SiFTO and the MLCS2k2-like approaches. MLCS2k2- 
like parametrisations assume a Cardelli-like colour variation law 
and in some cases a colour-luminosity relation, while the S ALT2 
or SiFTO approach allows one to derive these quantities from 
the data itself. The values obtained for these quantities differ 
significantly from the assumptions made in MLCS2k2. In ad- 
dition, colour priors are used in the MLCS2k2 approach, which 
we think are unnecessary and subject to additional systematics. 

To estimate systematic uncertainties arising for the light 
curve fitting process, we therefore chose to fit all SN light curves 
using both SALT2 and SiFTO models as described in the follow- 
ing section. 

4.3. Fitting with SALT and SIFTO models 

We use SALT2 and SiFTO to fit light curves and derive three pa- 
rameters for each SN: an amplitude conventionally described by 
the peak rest-frame B-band magnitude m* , a shape parameter, 
and a colour labeled "C" that roughly corresponds to the rest- 
frame B — V colour at maximum light. The two models that were 
fit on low-z and SNLS data minimise the modeling statistical un- 
certainty and cover a large wavelength range (300-700 nm) well 
suited for the SNLS data set. For both fitters, the light curves are 
obtained by integrating a spectral energy density (SED) varying 
with phase in a model of the instrumental response function. The 
effective colour parameter that is derived is due to a combination 
of extinction and intrinsic variation, and this has to be kept in 
mind in the subsequent usage of those parameters to derive dis- 
tances. Those two techniques differ substantially in their detailed 
parametrisation of observables and in the procedures considered 
for training and light curve fitting. As a consequence, in the fol- 
lowing, we use both in order to estimate the impact on cosmol- 
ogy of the modeling choices that were made. Before comparing 
the result of the fits obtained on the SNLS sample, we present 
both approaches, along with the minor modifications that were 
applied since the publication of the methods in G07 and C08. 
The two approaches are compared and combined in 95.11 

4.3.1. SALT2 

The SALT2 method consists in modeling the SNe la spectral en- 
ergy density (SED) variation with time and its diversity using a 
linear combination of several principal components multiplied 
by the exponential of a colour dependent function of wavelength 
(which we will call "colour variation law", although it can model 
a pure extinction law). This model is trained on a large sample 
of nearby SNe listed in Table [T2l and a sub-samble of the SNLS 
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Fig. 6. SALT2 colour variation law for C — 0.1 (thick solid 
curve), along with th e S ALT2 colour variation law of G07 and 
ICardelli etaO d!989t) extinction law for E(B - V) = 0.1 and two 
different values of Ry (1.0 and 3.1 for the dotted and dashed blue 
curves). 



SNe la listed in Table |9] We use only SNLS SNe with redshifts 
z < 0.7, with an unambiguous spectroscopic identification (ex- 
cluding the Ia-k) and a good light curve sampling (see 94.5b . 
Both light curves and spectra are used in the fit. Since signifi- 
cant calibration uncertainties are expected for most spectra (due 
for instance to slit losses), wavelength dependent corrections for 
each spectrum are included in the model. The model parameters, 
the calibration correction coefficients of the spectra, and the SNe 
parameters are fitted simultaneously. Only two principal com- 
ponents have been considered; one that represents the average 
SN la, and one which can be identified as the shape variation of 
light curves of SiFTO. This training sample is larger than the one 
considered in G07. A few technical modifications, detailed in 
Appendix|Aj have been applied to the training procedure: higher 
resolution for the components and the colour variation law, a new 
regularisation scheme, improved handling of the residual scatter 
about the model, and propagation of the model statistical uncer- 
tainties. The colour variation law is shown in Fig.[6]for C — 0.1 
(i.e. for a B - V colou r exces s of 0.1). It differs significantly 
from the ICardelli et al] ( fl989b extinction law for wavelengths 
A < 370 nm even when extreme values of Ry are considered. 
Adding more parameters to the colour law with respect to G07 
has resulted in a steeper variation at short wavelength. 

For the fit of SNe entering the Hubble diagram, all light 
curves for which the effective wavelength of the associated re- 
sponse function lies in the model validity range (300-700 nm 
rest-frame) are fit simultaneously. The actual fit is a simple least 
square minimisation taking into account both the covariance ma- 
trix of the flux measurements and a model uncertainty (described 
in 9A.2t . Since the uncertainty on the model depends on the pa- 
rameters that are fitted, the fit is performed iteratively with up- 
dates of the model uncertainties at each step until convergence 
is reached. The model parameters and the code to use them are 
public; see 9A.4l for details. 

4.3.2. SiFTO 

The SiFTO model consists of a SED sequence whose time evo- 
lution has been calibrated on a large SNe la sample combining 
low-z and SNLS data, starting from the spectral sequence de- 
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Fig. 7. (B - V) corrected for stretch as a function of(U-B) for a 
selected sample of SNLS SNe as measured by SiFTO (open cir- 
cles) along with the linear relations of Eq.|4]and Table [5] (a and 
c terms, solid lines) and the synthetised ones of SALT2 (dashed 
line) for the two colour scatter hypothesis. SALT2 does not con- 
tain explicitly colour relations as SiFTO, but the a, b and c coef- 
ficients can still be derived from the model. 



rived bv lHsiaoet ail (120071) . The light curve shape variability is 
modelled with a time-stretching of the SED sequence about the 
date of maximum light in rest-frame B-band, with a wavelength- 
dependent stretch factor (s) indexed by its value in B-band. The 
broad-band wavelength dependent calibration of the SED se- 
quence is not performed at this stage of the model building. 
Indeed, the derivation of the amplitude and colour parameters 
is obtained in two steps: i) Contrary to SALT2, the amplitude of 
each light curve is a fit parameter, leading to an observer-frame 
peak magnitude for each pass band, ii) These peak magnitudes 
are then used to adjust the SED at maximum-light with a smooth 
multiplicative function of wavelength so that the flux integrated 
in any rest-frame filter can be evaluated. Up to 5 rest-frame fil- 
ters are used - UmUBVR^l - resulting in 4 colour combinations 
of (U 2 -B), (U-B), (B-V) and (V-R), with the limitation that 
only the rest-frame filters within 450A of an observer-frame fil- 
ter are considered. The second step consists in converting these 
colours into an estimate of (B - V), denoted (B - V) pre d, using 
equations of the form 

(B - VOuBpred =a[(U-B) + 0.5] + b (s - 1) + c (4) 

with similar equations relating (B - V) with (E/02 - B) and (V - 
R). The data used to derive the linear relation of the (B - V) vs 
(U-B) relation is shown in Fig. [7] The estimation of the a, b, c 
parameters entering in these relations depends on the assumed 
intrinsic scatter of the SNe colours, this is discussed in 34.41 

The final SiFTO colour parameter C is determined by a 
weighted combination of each (B - V) pre d together with any ac- 
tual measurement of (B — V), degrading each (B - V) pre d by the 
dispersion measured for each relation. The uncertainties in the 
derived colour relations introduce correlations between SNe that 
are propagated in the analysis. 

Since C08, SiFTO has been retrained with a much larger 
sample of SNe including a sub-sample of the SNLS light curves 
presented in this paper (with stringent requirements on the time 
sampling) and many more at low reds hift. It also uses an up- 
dated version of the SED sequence from lHsiao et alJ (120071) that 
has been extended in the near infra-red. 



12 Ufa is an artificial filter defined as U blueshifted by z = -0.2 



As a conclusion to the description of the techniques, SiFTO 
does not contain an explicit colour variation law as a function of 
wavelength, nor a broad-band colour calibration of the SED se- 
quence as in SALT2. These two pieces of information, which are 
essential for distance estimate, are coded into the linear colour 
relations of Eq. [4] Also, the light curve stretching behaviour is a 
prior of the SiFTO parametrisation, whereas it is only approxi- 
mately realised as an outcome of the SALT2 training. Indeed the 
second principal component of SALT2 turns out to be close to a 
derivative of the first one with respect to a stretch factor, hence 
mimicking a time stretching at the first order. As a consequence, 
based on these differences, we do not expect to get exactly the 
same light curve shapes so that the peak magnitudes and colours 
obtained with both fitters may have small systematic offsets and 
a stretch dependence. Some additional dispersion due to the dif- 
ferences in the data weighting in the least square minimisation 
is also anticipated. 

4.4. Residual scatter 

The SALT2 and SiFTO empirical models do not fully account 
for the diversity of SNe la. As a consequence, some intrinsic 
scatter about the best fit model is expected on top of the measure- 
ment uncertainties. The wavelength dependence of this residual 
scatter has to be evaluated as it affects the determination of some 
of the model parameters. As an example, for SiFTO, the fitted 
slope of the relation between (U-B) and (B - V) depends on 
whether the scatter about the linear relation is attributed to the 
estimate of (U - E) or (B - V). The same effect applies to the 
determination of the S ALT2 colour variation law. Since a bias on 
this slope (for S ALT2, this slope is encoded in the second deriva- 
tive of the colour variation law in the UBV wavelength range) 
leads to bias on the average colour of SNe at high redshift where 
blue SNe are preferentially selected, we have to determine the 
amplitude of this scatter as a function of wavelength. 

This cannot be obtained in the SiFTO framework where only 
two colours are compared at a time to derive the model parame- 
ters. On the contrary, it is possible with SALT2 since it assumes 
a single relation between all the peak magnitudes of a given su- 
pernova (up to five for low redshift SNe: UBVRI). Note how- 
ever that with the SiFTO method, this unique relation between 
all observations is restored in the final colour estimate when the 
various (B - V) pie d are combined. The technical details of the 
determination of the residual scatter are given in 3A.2I Using 
two parametrisations of the scatter as a function of wavelength 
(the exponential of a polynomial, or a combination of sigmoids, 
labeled respectively hereafter "EXPPOL" and "SIGMOID", see 
Eq IA.2l i, one obtains the two estimates of the residual scatter 
shown in Fig. [8] This scatter model is translated into scatter in 
the reference bands used for SiFTO colour relations in Table [4] 
The values of the SALT2 residual scatter at the central wave- 
length of the U02BVR filters is used. The assumption that this 
scatter on magnitudes is uncorrelated between bands translates 
into correlated uncertainties in the colours used for the SiFTO 
colour relations. When accounted for in the determination of the 
colour relations, this leads to two sets of a, b, c parameters for 
each colour relation that are listed in Table[5] 

In this determination of the residual scatter, we identified an 
additional scatter of 0.1 mag in observer [/-band for low red- 
shift SNe (see also C10 for a discussion about this issue). This 
is significantly larger than the dispersion obtained in the same 
rest-frame wavelength range with SNLS data (given by gM band 
observations at z — 0.35, see TableHJi. We hence attribute this ad- 
ditional dispersion to calibration uncertainties or varying effec- 
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t/02 
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B 


V 


R 


EXPPOL 
SIGMOID 


0.179 
0.215 


0.063 
0.048 


0.029 
0.020 


0.019 
0.020 


0.025 
0.034 



Table 4. Residual scatter in SiFTO reference pass-bands for the 
two parametrisations of the scatter of Eq lA.2l 



EXPPOL 



SIGMOID 



(B- 


- V) vs (U 02 


-B) 


(B- 


- V) vs (U Q2 - B) 


a 


= +0.186 ± 


0.032 


a 


= +0.225 ± 0.032 


b 


= +0.371 ± 


0.134 


b 


= +0.443 +0.134 


c 


= -0.112 + 


0.015 


c 


= -0.130 + 0.015 


(B 


-V)vs(V 


-B) 


(B 


- V) vs (U - B) 


a 


= +0.734 + 


0.040 


a 


= +0.716 + 0.040 


b 


= +0.323 ± 


0.059 


b 


= +0.333 +0.059 


c 


= -0.109 + 


0.009 


c 


= -0.104 + 0.009 


(B 


- V) vs (V 


-R) 


(B 


- V) vs(V-R) 


a 


= +1.755 + 


0.114 


a 


= +2.011+0.114 


b 


= +0.411 ± 


0.061 


b 


= +0.473 ± 0.061 


c 


= -0.885 ± 


0.060 


c 


= -1.018 + 0.060 



Table 5. Parameters of the SiFTO colour relations on (B - V) 
obtained with the two residual scatter estimates of Table [4] 



tive filter transmissions from one U -band light curve to another. 
We de-weighted those light curves accordingly in the training of 
SALT2, and subsequently in the light curve fits. As for SiFTO, 
they were discarded for the fit of the (U - B) (B - V) colour 
relation and were not used for the determination of the colour 
parameters. 
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Fig. 8. S ALT2 training sample broad-band colour dispersion as a 
function of wavelength for two different functional forms (solid 
curve: EXPPOL, dashed curve: SIGMOID). 



Figure Qpresents the difference between SiFTO and SALT2 
colour relation between (U - B) and (B - V). Part of the differ- 
ences of slopes can be explained by the use of different samples 
to derive those relations, and part can be attributed to the intrin- 
sic differences in the methods used to derive those relations. Still, 
one can see that they have the same average estimate of (B - V) 
given (U - B), so we do not expect this to have a significant 
impact on cosmology, as we will see in the next section. 



4.5. Light curve sampling 

We investigate in this section the reliability of the light curve pa- 
rameters determined on poorly sampled light curves, in particu- 
lar those which lack photometric measurements before the date 
of maximum light (hereafter T max ). Any bias on T max induces a 
bias on the parameters needed to estimate a distance (magnitude, 
shape and colour), and such a bias may arise if the modeling of 
the decline of the light curves is not sufficiently accurate. This 
is not really an issue for the SNLS sample presented in this pa- 
per as only 2% of the SNe lack photometry in the rising part of 
the light curves (between -10 and -1 days with respect to T max , 
rest-frame), but it has to be studied for the selection of external 
samples (see C10) for which this fraction is much larger (es- 
pecially at low redshift because of a follow-up triggered by the 
spectroscopic identification). 

For the purpose of this study we select well sampled SNe 
data with high S /N from SNLS (using only those at z < 0.4) 
and from external samples: SNe at z < 0.1 from various surveys 
(see re ferences in C10) and SNe at 0.06 < z < 0.3 from the 
SDSS dHoltzman et alj|2008l) . We compare the light curve pa- 
rameters 9 E {m* B ,Xi,C,fi} derived from the full light curve fits 
- hereafter 9f - to those resulting from fits ignoring a fraction 
of the data points. For each SN, we derive several estimates of 9 
removing sequencially an increasing number of early points, and 
index each of those estimates by ry = (Tfi rst - T mnx )/(l + z) the 
phase of the first remaining data point. Of course for real data ry 
is not directly observable and we have to rely on ry, an estimate 
of Tf which results from the fitted value of T max . We study the 
average offsets A9 = 9{tj) - 9f as a function of ry. One first has 
to correct for a trivial bias (showing up even for an unbiased es- 
timator) that is due to the finite width of the tj distribution. For 
instance, if we consider an input distribution with tj > 0, we 
will find biased A9 for negative values of 77 simply because we 
are looking at events for which fy is systematically underesti- 
mated (and hence T max over-estimated). Within each bin around 
a given value of ¥}, we have an average bias on a parameter 6 
given by 



A,6(T f ) = 



jg* g £Z? = l ["-4] drOpiu^du 



where the sum is on the number of tests performed (indexed by 
Ty), and p(T/\Tf) the likelihood of the estimator ry knowing the 
true value Tf (e is half the bin size). d T = -(1 + z) d9/dT max 
is the derivative of 9 with respect to jy. Such a bias can be 
estimated for each value of ry and corrected for assuming a 
Gaussian likelihood (with a cr given by the redshift corrected 
measurement uncertainty on T max ). We will consider in the fol- 
lowing A' 9 - A9 - A,9 to study potential biases. 

The results for the rest-frame peak magnitude in B-band, the 
peak (B - V) colour and shape parameter are shown in Fig. [9] 
The potential bias on the distance modulus fi is also shown (us- 
ing best fit values for a and B, see Eq. [5] and 35.2b . For light 
curves without premax data but with ry < 5 days, we do not 
find significant biases on m* B and (B - V) but a small bias on X\ . 
The uncertainties on those biases are difficult to estimate, as the 
de-biasing is only approximate. 
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Fig. 9. From top left to bottom right, A'm* B , A'C, A'Xi and A'fi 
for SALT2 as a function of f} the estimated phase of the first 
measurement. 

For subsequent analyses, we will select SNe according to the 
following minimum requirements (based on the available phases 
t = {T bs - T ma x)/(1 + z) of photometric observations): 

(i) Measurements at four different epochs or more in the range 
-10 < t < +35 days. Three are mandatory to estimate both 
the date of maximum light and the shape parameter, one more is 
needed to get at least one degree of freedom. 

(ii) At least one measurement in the range -10 < r < +5 days 
(equivalent to tj < 5 days). 

(ii) At least one measurement in the range +5 < r < +20 days 
for a reasonable evaluation of the shape. 

(iii) At least two bands with one measurement or more in the 
range -8 < r < +10 days in order to evaluate the SN peak 
luminosity and colour with confidence. 

With this selection applied to the test data set described 
above, the estimated bias on distance moduli for the selected 
SNe without premax data is compatible with zero (A'p = 
-0.004 ± 0.004). However an additional uncertainty due to the 
de-biasing of order of 0.01 cannot be excluded. We use this lat- 
ter number as a systematic uncertainty on the estimated distance 
moduli of SNe without premax data. 

5. Light curve parameters of SNLS SNe la and 
measurement of n M 

We present in this section a fit of the SNLS Hubble diagram 
based on the results obtained with the two light curve fitters de- 
scribed above. The goals of this study are: i) to compare the 
SALT2 and SiFTO light curve parameters, ii) see how their dif- 
ferences impact distance estimates and cosmological fits, iii) 
propagate the uncertainties of the whole analysis chain to the 
cosmological results, iv) present the constraints on cosmology 
provided by a single survey of high redshift supernovae. 

5.1. Light curve parameters of SNLS SNe la 

We apply in this section the two techniques described above to 
determine the peak magnitudes, colours and light curve shapes 
of the SNLS SNe la. Those are the basic ingredients to deter- 
mine distances, they are used in C10, along with external SNe 
data samples, to build a Hubble diagram and constrain cosmo- 
logical models. Here we focus on consistency tests and compar- 
isons of the outcome of the two light curve fitters, any difference 



being considered as a systematic uncertainty that will add to the 
photometry and calibration uncertainties discussed in $3] 

In order to get reliable estimates of magnitude, colour and 
light curve shape for each SN, we apply the selection described 
in {34. 51 This set of sampling cuts discards 25 SNe so that we are 
left with 252 SNe la/la*, excluding the peculiar SNe 03Dlcm, 
03D3bb, 05Dlby and 05D3gy. The discarded SNe are identified 
by the label "(s)" in Tabled 

5.1 .1 . Internal consistency check of the light curve fitters 

In this section we try to quantify one of the primary virtues 
of the light curve fitters, which consists of deriving redshift- 
independent parameters (magnitude, shape and colour), from ob- 
servations in a limited set of filters. This is possible in SNLS 
thanks to the high quality light curves obtained in gM^M and 
im bands. For supernovae in the redshift range [0.2,0.7], this 
corresponds to observations spanning the rest-frame wavelength 
range 2900-6400 A. 

The fitters should on average give the same estimates of the 
parameters for fits based on either gM+fM or r^+iM light curves. 
The results of such a test are shown in Fig.[10]for the difference 
A(B - V) of the estimates of rest-frame (B - V) colour, dispersed 
as a function of redshift. The average values of A(B - V) in red- 
shift bins are compatible with zero (the maximum deviations are 
respectively of 0.024 ± 0.02 and 0.020 ± 0.012 for SALT2 and 
SiFTO). This tells us that we can reliably use rest-frame U and B 
observations to estimate the (B-V) colour at high-redshift where 
we do not have rest-frame V-band observations. Note however 
that the statistical uncertainties on the colour relations are ac- 
counted for: uncertainty on the c parameter of the colour relation 
for SiFTO (of order of 0.01 for (U — B)), and co variance matrix 
of the training for SALT2 propagated to the distance moduli as 
detailed in COl 



5.1 .2. Comparison of the SALT2 and SiFTO light curve 
parameters 

We compare in this section the output of the two light curve 
fitters. The results obtained with SALT2 and SiFTO are listed 
in Table QT| Since those empirical models were primarily de- 
signed for distance estimates, it is important to briefly review 
how the three resulting parameters (rest-frame magnitude, shape 
and colour) are combined for this purpose. As in A06, a distance 
modulus is defined by the following linear combination of the 
parameters: 



M + ffX shape -BxC 



(5) 



where m* B is a rest-frame B magnitude, the shape parameter is 
(s - 1) for SiFTO (where s is the stretch factor), Xi for SALT2, 
and C is the colour parameter. The absolute magnitude M, and 
the linear coefficients a and B are fitted simultaneously with 
the cosmological parameters (and also marginalised over for the 
cosmological constraints). This approach will be applied and 
discussed in detail in C10. We simply note here that any con- 
stant systematic offset between SiFTO and SALT2 parameters 
will have no effect on cosmology as they will be absorbed in the 
absolute magnitude, that any multiplicative factor between the 
SiFTO and SALT2 estimates of the shape and C parameters will 
be compensated for by the a and B coefficients, and finally that a 
shape dependence of the differences of m* B and C estimates with 
SiFTO and SALT2 are also accounted for by the a term. 
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0.2 




a* = -0.008 ± 0.005 

a° B = 0.013 ± 0.004 

a c = 0.997 ± 0.097 

a* = 0.002 ± 0.009 

a° = 0.035 ± 0.008 

af =0.107 ±0.006 

a° = 0.991 ± 0.006 

Table 6. Coefficients of the linear transformation of SALT2 pa- 
rameters to match those of SiFTO as defined by Eq.|6] 



Whereas the shape and colour parameters do not show any 
significant deviation from zero when averaged in redshift bins, 
there is a systematic difference on mg for some redshifts (-0.02+ 
0.005 at z 0.7 and +0.02 + 0.005 at z 0.9). This points 
to different effective ^-corrections between SiFTO and SALT2. 
Those differences can be attributed to some discrepancies in the 
sp ectral sequences u sed (that of SALT2 and the one derived 
by iHsiao et ail |2007|) . as the spectral distortions performed in 
SiFTO do not compensate for differences of fluxes on a wave- 
length range much shorter than the gaps between pass-bands (see 
the red curve on Fig. QT|top panel). However the choice of the 
empirical modeling of the SN diversity (global stretching or lin- 
ear additive components) can also play a role, as the differences 
in the training sample and technical details of the training pro- 
cedures. We will hence treat those differences as an additional 
source of systematic uncertainty. 



-0.2 



0.3 0.4 0.5 

SN Redshift 
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Fig. 10. Differences of estimated (B - V) colour at maximum 
from the fit of ymIm or gMfM light curves for both fitters 
(top: SALT2, bottom: SiFTO). 



As a consequence, in order to compare of the outcome of 
the two fitters, we will first transform the SALT2 parameters to 
match those of SiFTO with the following equations: 



i* B + «g Xi + a 



C = a c c C + a x c X x +a° c 
s' = a^Xi+a ^ 



(6) 



a x s and a° s are inevitable because of the fundamentally different 
modeling of the variability of light curve shapes. a° c and a c c are 
due to the fact that for SALT2, C is an internal parameter of the 
model that only approximately correspond to a (B - V) colour at 
maximum, with a constant term so that the average value of C 
on the training sample is zero, whereas for SiFTO C is exactly 
the (B - V) colour at maximum. a° B , and a* can be different 
from zero since the shape of the B and V bands light curves are 
not exactly the same for SALT2 and SiFTO. The values of the 
coefficients obtained by the comparison of the data of Table QT| 
are listed in Table [6] We apply this transformation of SALT2 
parameters to focus on differences between the two light curve 
fitters that lead to biases on distance moduli (Eq. [5jl. The values 
of the transformation coefficients themselves are simply related 
to differences in the definition of the parameters in both models 
and do not contain any meaningful information. 

The differences of those transformed magnitude, shape and 
colour parameters of SALT2 with those obtained with SiFTO are 
dispersed as a function of redshift in Fig.fTTI 



5.2. SALT2 and SIFTO based Hubble diagrams 

We present in this section fits of the Hubble diagram based on 
SALT2 and SiFTO parameters. SALT2 parameters are trans- 
formed to match those of SiFTO (see Eq.|6]and Tabled- We use 
the distance modulus fi given in Eq.[5]with the shape and C pa- 
rameters matching the SiFTO stretch parameter and the (B - V) 
colour at maximum. We consider here a ACDM model, com- 
posed of unrelativistic matter and a cosmological constant, and 
further assume a spatial flatness, so that we are left with a sin- 
gle cosmological parameter to fit: Q.m (or alternatively Qa = 
1 - Qm)- We apply the same procedure as the one described in 
A06. The estimation consists of minimising 



X 



S 



(M /^cosmo) 



(7) 



where // cosmo = 5 log 10 [d L (Q M ,z)/10pc], d L being the luminos- 
ity distance, and cr,„, is an additional dispersion which value 
is chosen in order to obtain a reduced y 2 — 1 ■ This additional 
dispersion accounts for all the sources of variability of SNe lu- 
minosities beyond their first-order correlation to the shape of 
light curves and the colours of SNe. The measurement uncer- 
tainties cr{p) are a function of the parameters a and B. We use 
Ho = 70 kms -1 Mpc~'. This is purely conventional as it only 
alters the value of the absolute magnitude M in Eq. [5] through 
the normalisation of di, and M is fitted simultaneously with the 
cosmological parameters. 

In addition to the sampling cuts listed in 94.51 we discard 
SNe with a peak rest-frame (B — V) > 0.2. Those red SNe are 
found only at z < 0.6 in our sample because they are fainter 
than the average, and hence undetected (or unidentified spectro- 
scopically) at higher redshifts. Discarding them minimises the 
potential bias on distance moduli due to an inadequate colour 
correction. Indeed the average colour correction coefficient (J3) 
we derive from the bulk of the SNe may not apply to those red 
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Fig. 11. Comparison of SALT2 and SiFTO parameters for a se- 
lected sub-sample of SNe of Table fTTI Linear transformations of 
Eq.|6]and Table[6]were applied to SALT2 parameters. Each dot 
is a SN, the light grey areas represent the RMS of the distribu- 
tion in redshift bins, and the dark areas the uncertainty on the 
average values. The red curve on the top panel for the rest-frame 
B band peak magnitude is the expected discrepancy due to the 
differences of the SALT2 and SiFTO spectral sequences. 



SNe that are more likely to be extinguished by dust in their host 
galaxy than the bluer ones. This cut, applied to both SALT2 and 
SiFTO samples, discards 14 SNe0 We mention in ^5A3l the 



13 The 14 SNe discarded because of their colour are at z < 0.71. For 
those redshifts, the average uncertainty on the colour is 0.03. 





CAT T1 




( \ 


n i -i- n cyxk 

U. LyO ± U.UjO 


u.zij ± v.vjj 


M 


-19.218 ±0.032 


-19.210 ±0.030 


a 


1.295 ±0.112 


1.487 ± 0.097 


P 


3.181 ±0.131 


3.212 ±0.128 


RMS 


0.173 ±0.008 


0.150 ±0.007 


0~int 


0.087 


0.087 



Table 7. Qm and the nuisance parameters with their statistical 
uncertainties derived from the Hubble diagram fits using either 
SALT2 or SiFTO. SALT2 parameters were linearly transformed 
(see Eq.|6]and Tabled; in particular, the a coefficient applies to 
stretch and not the original SALT2 X\ parameter. 



impact of a change of this colour cut on the cosmological fit. 
In addition to this, the peak magnitudes of 3 SNe could not be 
obtained with SiFTO because of a lack of observations in gM 
and zm bands. Finally, a 3<x clipping to both SALT2 and SiFTO 
Hubble diagram residuals was applied, which removes 7 more 
SNe (4 for SiFTO and 4 for SALT2 with only one in common, 
the others being consistent with a statistical fluctuation). 

A Malmquist bias correction is applied to the distance mod- 
uli. This correction is based on a detailed study of the survey 
detection and spectroscopic identification efficiency using fake 
SNe included in the im ages entering the detec tion pipeline. This 
analysis is described in lPerrett et al.l din prep.h . The correction is 
equal to zero for z < 0.4 and rises at higher z reaching a value of 
^ 0.025 magatz= 1. 

The results obtained with both light curve fitters are listed 
in Table [7] The differences in the light curve fitters lead to a 
systematic difference of 0.020 on that corresponds to about 
two thirds of the statistical uncertainty. The RMS of the Hubble 
residuals are similar, SiFTO residuals being slightly less scat- 
tered than the SALT2 ones. The Hubble diagram residuals from 
the corresponding best fit models are shown in Fig.[T2] The cor- 
relation of the residuals from the Hubble diagram obtained with 
both fitters is shown in Fig. [13] The value of the correlation co- 
efficient of 0.82 shows that most of the residual scatter is due to 
variability beyond that addressed by a single shape coefficient 
and a colour, irrespective of the way those parameters were de- 
rived. 

5.3. Combination of SALT2 and SiFTO results 

The aim of using two light curve fitters is to evaluate the uncer- 
tainties on distances associated with the choices made in the de- 
velopment of those empirical models. In the following, as we do 
not have any physical justification for choosing one or the other, 
we use as a central value the average of the two, and propagate 
the systematic differences of the two fitters to the uncertainties 
on cosmological parameters. 

We use for this purpose the method described in C10. In 
a few words, all systematic uncertainties are accounted for in 
the covariance of the distance moduli. The distance modulus of 
Equation[5]can be rewritten as yU = 6 T T]- M where 6 = {1, a, -B] 
is the vector of the nuisance parameters, and rj = {m* B , s - 1, C) 
is the vector of the SNe parameters. With H, the Jacobian matrix 
of the derivative of 77, with respect to the systematic uncertain- 
ties, and C sys the covariance matrix of those uncertainties, the 
covariance of the distance moduli of two supernovae is given by 

CovUa, fij) = 6 T (SijC v , + ///C,.,,//;) 6 (8) 
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Fig. 12. Residuals from the Hubble diagram of the best fit 
ACDM cosmological model obtained with the two light curve 
fitters SALT2 (top) and SiFTO (bottom). 



where C, h is the covariance matrix of the measurements r\-, (av- 
erage of SALT2 and SiFTO results, including the intrinsic cr„„ 
obtained in 95.21 with each fitter independently) assuming that 
the SALT2 and SiFTO results are fully correlated. 

For the systematic uncertainty addressed here, we assign to 
H the difference of S ALT2 and SiFTO parameters, averaged over 
a redshift range about each SIS0, and consider a variance of one 
for this uncertainty in the covariance matrix C svs . 

When fitting for a flat ACDM cosmological model, we ob- 
tain with this procedure Qm = 0.201 + 0.043, where the uncer- 
tainty includes a statistical uncertainty of 0.034 and an uncer- 
tainty due to the difference of the fitters of 0.026. If we now 
release the flatness prior, we obtain the confidence contours in 
the plane shown in Fig. [14] The combined confidence 

contour contains the ones obtained with SiFTO or SALT2 only. 
Those results do not include the full sources of systematic un- 
certainties, which will be detailed in Sects. [5~4l and l5.5l and final 
results will be presented in Sect. 15.61 

5.4. Systematic uncertainties from light curve fits 

We review here the list of systematic uncertainties that affect the 
determination of the light curve parameters and enter the covari- 
ance matrix C sys . 



14 We use for the averaging about each SN redshift a Gaussian filter 
with cr : = 0.05. 
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Fig. 13. Correlations of residuals from the Hubble diagram ob- 
tained with SiFTO and SALT2. The outlier is 04D2cc for which 
data in the rising part of light curves is missing; the statistical 
uncertainty on its distance modulus // is large with both light 
curve fitters (0.47 with SiFTO and 0.28 with SALT2, the dif- 
ference being due to differences in the model uncertainties, or 
"error snake"). 



5.4.1. Calibration systematic uncertainties 

The calibration systematics have been identified and evaluated 
in R09, they are listed in Table [8] All of them affect the light 
curve parameters. They can be divided in two sets: covariant 
magnitude uncertainties for each of the gMfM^MZM bands and 
uncertainties on the modeling of filters. The former magnitude 
uncertainties comprise uncertainties on the calibration transfer 
from Landolt stars to tertiary stars, from tertiary stars to super- 
nova light curves, and uncertainties on the magnitude to flux 
transformation given by the uncertainties of our flux standard 
BD+17 4708 magnitudes and its spectrum. Additional uncer- 
tainties come from the calibration of the PSF photometry as dis- 
cussed in 93.21 

The uncertainties on the transmission function of the instru- 
ment are approximated by uncertainties on the central wave- 
length of each effective filter (an effective filter combines quan- 
tum efficiency of the CCDs, transmission of the MegaPrime op- 
tics and filters, reflectivity of the primary mirror, and a model of 
the atmospheric transmission at the average airmass of the ob- 
servations, see R09); any higher order uncertainty on the effec- 
tive transmission has a negligible impact on the total uncertainty 
budget. 

In order to propagate these uncertainties, the derivative of 
the peak magnitude, colour and stretch (or X\ for SALT2) with 
respect to each of these systematics (H matrices in Eq. [8} have 
to be evaluated for each SN. For SALT2, this requires a full re- 
training of the model and new light curve fits for each kind of 
uncertainty. For the SiFTO model, calibration uncertainties only 
impact the linear colour relations. 

As an example, the impact of the calibration uncertainties on 
the values of m* B are shown in Fig.[l5]for SALT2; they are similar 
for SiFTO. The sensitivity of m* B to calibration offsets depend on 
the weight of each band in the fit which in turn depends on the 
assumed model uncertainties. 
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Table 8. Calibration systematic uncertainties 





gM 


I'M 


lM 


Zm 


MegaCam magnitude system 1 " 


±0.005 


±0.005 


±0.008 


±0.019 


Low S/N PSF photometry bias 


< 0.001* 


< 0.002 


< 0.001 


< 0.001* 


Calibration of the PSF photometry 


±0.002 


±0.002 


±0.002 


±0.002 


Total 


±0.006 


±0.006 


±0.008 


±0.019 


Central wavelength of effective filters 


±7A 


±7A 


±7A 


±25A 



+ including the uncertainties on the conversion of magnitudes to fluxes obtained with BD+17 4708 spectrum, see Table 12 of R09 
* potential bias due to the uncertainty on the weighted average position from r M and i M fit. 
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Fig. 14. Contours at 68.3% confidence level (accounting only 
for statistical uncertainties) for the fit of ACDM cosmology to 
the SNLS third year sample (without additional SNe samples) 
for SALT2 (thin red curve), SiFTO (thin blue) and the combi- 
nation of the two (thick black), on the light curve parameters. 
The d otted inclined ellips e represents the 68.3% confidence level 
from lAstier etail (12006?) which includes 44 low redshift SNe 
along with the 71 SNe from the SNLS first year sample. 



5.4.2. Light curve fitter systematic differences 

There are three sources of uncertainties associated with the light 
curve fitters. First, there are the intrinsic differences of model- 
ing given by the differences of SALT2 and SiFTO outcomes as 
discussed previously. 
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Fig. 15. Impact of the calibration uncertainties of Table[8]on the 
SALT2 values of m* B as a function of redshift (gM and r^-light 
curves are not used on the whole redshift range because of the 
limited rest-frame wavelength range of SALT2). 

Second, due to finite training samples, there is some statisti- 
cal uncertainty associated with the determination of the param- 
eters of the light curve fitters. For SALT2, those uncertainties 
are provided by the covariance matrix of the model that results 
from the training fit (taking into account the residual scatter). 
The transformation to a covariance matrix of distance moduli is 
detailed in jHA.3l as is the computation of the uncertainty on the 
average distance modulus in a given redshift bin. As for SiFTO, 
only the uncertainties on the a,b,c parameters of the colour rela- 
tions are considered. The uncertainties of the spectral sequence 
(and especially its evolution with time) are ignored, because they 
induce negligible correlations of parameters derived for different 
supernovae: an "error snake" is used to fit light curves, but the 
covariance between different SNe is ignored. 

Finally, as discussed in 34.41 there is some uncertainty on 
the wavelength dependence of the residual scatter about the best 
fit model. We have considered two functional forms for this de- 
pendence which provide us with an estimate of the associated 
systematic uncertainty. 

An additional systematic uncertainty specific to the SALT2 
technique is the regularisation weight. Regularisation is required 
to train the model because of the lack of spectra in the near UV 
for some SN phases. Its weight should be minimal to prevent 
biases, but a weight that is too low generates unphysical high 
frequency wiggles in the resulting spectral sequence (this is a 
typical feature of all deconvolution processes). So the choice of 
the regularisation scale is partly arbitrary. Its impact on distance 
moduli is evaluated by the difference of results obtained with 
two sets of SALT2 trainings, one with the nominal weight and 
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the other one with a weight five times smaller. However as shown 
in Fig. [16] the differences obtained are typically of order of 0.005 
and hence can be neglected. 

5.4.3. Resulting uncertainties on distance moduli 

Figure [16] presents the impact of light curve fitter and calibra- 
tion uncertainties on distance moduli. The uncertainty induced 
by the choice of light curve fitter is estimated by the differences 
of distance moduli obtained with SiFTO and SALT2 averaged 
over redshift bins of size 0.2. The values of M, a and B ob- 
tained in 35.21 have been used. Because of a finite supernova 
sample, the uncertainty induced by the differences of the light 
curve fitters cannot be precisely estimated. We can however see 
that they are of the same order as the calibration uncertainties, 
with A// < 0.03. 

The statistical uncertainty of the light curve model and the 
choice of parametrisation of the residual scatter only play a sig- 
nificant role at the high redshift end of the survey. This is due 
to the combination of Malmquist bias and a limited wavelength 
coverage of the imaging survey. Because of Malmquist bias and 
the colour-luminosity relation, only the bluest supernovae are 
measured with enough signal to noise at high redshift, which 
results in an average colour decrease with increasing redshift. 
In addition to this, the rest-frame (B - V) colour is not directly 
measured, it is rather given by measurements in rest-frame wave- 
length shorter than that of the V-band so that one has to rely on 
the colour relations. As a consequence, any uncertainty on the 
colour relations will impact distance moduli, and especially a 
change of slope of the (U - B) or (t/02 - B) colour relation due 
to a modification of the residual scatter model will modify the 
average colour and impact distance moduli. 

Along with the uncertainties related to the light curve param- 
eters, several other systematic uncertainties alter the precision of 
the cosmological constraints resulting from the Hubble diagram 
fit. They are studied in detail in C10. Among those, we only 
consider the most relevant one for the present analysis (based 
on SNLS SN sample only): a 39% uncertainty on the Malmquist 
bias correction (see C10 for details). The other systematic un- 
certainties presented in C10 are negligible for this analysis. As 
an example, C10 found that the bias on distance moduli due to a 
contamination of the SN la sample by core collapse supernovae 
reaches a value of 0.005 mag at z — 1. This is much smaller than 
the other uncertainties we consider in this paper (see Fig. [TBI). 

The choice of the colour cut value of 0.2 applied to the SN 
sample is partly arbitrary. Choosing instead a cut of 0.3 for in- 
stance adds four SNe to the sample. With those additional su- 
pernovae the best fit Q„, value is shifted by 0.01. The statistical 
fluctuation of this difference of Q.,„ is +0.005 (based on a jack- 
knife analysis), so the shift we obtain is marginally consistent 
with a statistical fluctuation. 

5.5. Evolution of the colour-luminosity relation with redshift 

Recently, iKessler et ail (12009b (hereafter K09) found an evolu- 
tion of the B parameter with redshift, based on light curve pa- 
rameters resulting from SALT2 fits of the SNLS first year SNe 
data set (A06). Before presenting the final estimate of Qm (in- 
cluding all previously listed systematic uncertainties), we have 
to address this issue as our cosmological constraints are based 
on the assumption that B is constant. 

The measurement of the colour-luminosity relation at high 
redshift is difficult because of three combined effects: 
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Fig. 16. Uncertainties on the average distance modulus fj. in red- 
shift bins of 0.2: impact of the statistical uncertainty of the train- 
ing (for SALT2, thin solid curve), calibration uncertainties (dot- 
ted curve), residual scatter model (dotted short dashed curve), 
systematic uncertainty due to SALT2 regularisation (dotted long 
dashed curve) and differences between results obtained with the 
two light curve fitters (thick solid curve). Values of or and B that 
minimise residuals from the Hubble diagram were used. 

i) The fitted value of B depends on the assumed uncertain- 
ties of_J^6_^_J^]_(kj20wn_as a total least-squares problem, 
see lMarkovskv & Hu ffel 2007 for a review). 

ii) Because of a limited detection/selection efficiency, we miss 
red SNe at high z, and hence the range of colours to fit for the 
relation is reduced. This makes the determination of B both more 
uncertain and more sensitive to the assumed uncertainties. 

iii) SNe rest-frame (B—V) colour is not directly measured at high 
Z. We measure a colour at shorter rest-frame wavelength, and 
rely on colour relations to estimate (B - V). Incorrect colour re- 
lations cause the apparent/? to change with redshift. Furthermore 
the (B - V) estimate depends on assumed uncertainties of the SN 
model and intrinsic scatter of SN colours, which are difficult to 
evaluate (see 94 .41 . 

Hence, the fit of B not only requires an unbiased estimate of 
(B—V) but also of the uncertainties of this estimate (including all 
sources of scatter): underestimated uncertainties on (B—V) result 
in a negative bias on B. The uncertainties of colour relations and 
colour scatter should be propagated. 

With the improved modeling of the residual scatter presented 
in this paper with respect to the first version of SALT2 used by 
K09, we obtain larger uncertainties on the colour parameter and 
hence larger B at high z than those presented in K09. We re- 
peat their analysis which consists of fitting for M, a and B in 
redshift bins assuming a fixed cosmological model. The results 
are shown in Fig.[17]using both SALT2 and SiFTO parameters. 
Although our redshift-/? relation is much shallower than that ob- 
tained by K09, we still see a B evolution with the S ALT2 param- 
eters, B being smaller at higher z. 

This change in the estimate of B at high redshift between this 
analysis and that of K09 comes from the combination of two 
effects. 

i) The SALT2 colour variation law used in this paper is 
steeper than that of G07 for A < 3600 A mostly because of a 
higher number of parameters in the model training (see Fig. [6}. 
At z — 0.8, a change of the observed (tm - iia) colour translates 
into a change of the colour parameter 8C t hiswork — 0-66 x 6(tm - 
im) with the version of SALT2 used in this paper whereas, using 
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that of G07, one obtains 5Ccm - 0.81 x 6(r M - i'mH Hence, 
8C this work /dCcxniz = 0.8) = 0.8 and because of this, we expect 

ii) For supernovae at z > 0.7, which correspond to an aver- 
age redshift z - 0.8, we obtain an average uncertainty on SNe 
colours of cr c = 0.06 in this paper whereas one obtain an uncer- 
tainty of 0.05 when fitting the same data with the SALT2 version 
presented in G07 and used in K09. If we apply the scaling factor 
on the colours of 0.8 coming from the change of the colour vari- 
ation law, the uncertainty on colours reduces to <x c =a 0.04. This 
difference of colour uncertainties (0.06 versus 0.04) comes from 
differences in the model uncertainties that are shown in Fig.|8]in 
this paper and Fig. 7 in G07. Using a Monte Carlo simulation of 
a linear fit of the colour-luminosity relation for SNe at z > 0.7, 
where we consider in the fit that cr c = 0.04 when the value is set 
to 0.06 in the simulation, we find a bias /?/,- ( -fitrue — -10 

Combining the two effects (change in the estimate of colours 
and their uncertainties cr c ) results in a change of B of /3 t hi S wor k - 
Pcoiiz = 0.8) =2 1.6. At lower z, where we measure directly the 
rest-frame (B - V) colour and benefit from a larger range of SNe 
colours to fit the slope of the correlation with luminosity, we do 
not expect any significant difference. 

Using SiFTO parameters, we see in Fig. [17] the opposite 
trend than that obtained with SALT2: larger B at higher z (at a 
3<x significance based on a linear fit, but again, the statistical un- 
certainties are not meaningful), which illustrates how sensitive 
this measurement is to the light curve fitter modeling and error 
propagation. 

As a consequence, given the current uncertainties in the em- 
pirical light curve models, we are not able to conclude on an 
evolution of B with redshift. To account for this, we have added 
a systematic uncertainty on the evolution of B parametrised by its 
first derivative with redshift, for which we assume cr(dB/ dz) = 1 . 

5.6. Measurement of Cl M 
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Fig. 17. a and B estimates in redshift bins of 0.2 using SALT2 
(filled circles), SiFTO (open circles) and the combined light 
curve parameters (gray triangles, see 35.3b . 



We propagate the sources of uncertainties listed above with the 
technique described in C10 and summarised in 35.31 Except for 
the uncertainty associated with the difference of the two light 
curve fitters, the derivatives in the matrix H are the average of 
those obtained with SALT2 and SiFTO. For instance, the two 
sets of parameters T] a (p) result from fits with different weighting 
of the gMfMiMZM light curves, so that they have different depen- 
dences on the calibration systematics. 

Recently ISullivan et alj dsubmittedl) found with a Act evi- 
dence that the distance moduli (rest-frame magnitude corrected 
for the correlation with stretch and colour) of SNe in massive 
host galaxies are smaller than those of SNe in faint galaxies. In 
order to account for this effect in the cosmology fits, a second 
normalisation parameter Mb is introduced in order to treat inde- 
pendently SNe in low or high mass galaxies (with a split value at 
a galaxy mass of 10 10 M o ). When we apply this procedure to the 
SNLS sample, the retrieved value of Q.m is shifted by 0.0035. So 
this effect is negligible for the present analysis based on SNLS 
data only. 



At z - 0.8, the z M -band light curve, which have a poor signal to 
noise ratio, have a low weight for the colour estimate. 

16 We consider in the Monte Carlo simulations a value p = 3, 96 SNe, 
an intrinsic dispersion of distance moduli of 0.87, and a dispersion of 
the true SNe colours of 0.05 (the RMS of the observed colour distribu- 
tion is of 0.08, so that one expects an RMS of the true colours of order 
of 0.05 after subtracting the contribution of uncertainties). 



The contours in the Dm. plane are presented in Fig. Q~8] 
With a flatness prior (Qm + £2a = 1), we obtain Qm = 0.211 ± 
0.077, the uncertainty of which can be approximately decom- 
posed as 

Hubble diagram statistical uncertainty ±0.034 
Choice of empirical model (SALT2 vs SiFTO) +0.026 
Model training statistical uncertainty ±0.034 
Photometric calibration ±0.048 
Uncertainty on Malmquist bias correction (40%) ±0.004 
Uncertainty on the residual scatter (see 34. 4t ±0.010 
Potential B evolution (cr(dB/dz) = 1) ±0.022 
SNe brightness vs host galaxy mass ±0.003 

When the potential /? evolution is accounted for, the best fit 
(£2m, Qa) is shifted. This is due to a modification of the co vari- 
ances of the distance moduli and not their values. Modifying the 
covariances alters the weight of SNe as a function of redshift and 
hence the best fit cosmological parameters. 

Of course, when SNe samples at lower redshifts are added 
to the cosmological fit, both the statistical and systematic un- 
certainties are significantly reduced (though other systematics 
have to be included), and the relative impact of each uncertainty 
changes. This is studied in detail in C10. 
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Fig. 18. Contours at 68.3% confidence level for the fit to a 
ACDM cosmology from the SNLS third year sample (without 
additional SNe samples) accounting for: Hubble diagram statis- 
tical uncertainties and the difference between SALT2 and SiFTO 
results (dotted curve), photometric calibration (dashed curve), 
model training statistical uncertainties and Malmquist bias un- 
certainties (thin solid), and finally the potential /3 evolution (thick 
solid curve). 



6. Conclusion 

We have presented photometric properties and distances mea- 
surements of SNe discovered during the first 3-year of the SNLS 
project. This paper is the first of a series of three: the two oth- 
ers present the cosmological constraints obtained by combining 
this data set with other SNe la samples (CIO) and taking into ac- 
count co nstraints from other cosmological probes dSullivan et alJ 
lin prep.h . The various steps of the analysis were cross-checked 
using significantly different techniques, in order to identify and 
minimise systematic uncertainties. The photometric reduction 
was done using two independent pipelines. The calibration trans- 
fer from tertiary stars to supernovae light curves was performed 
with systematic uncertainties of 0.002 mag. Two light curve fit- 
ters were considered in this paper. Both are intended to esti- 
mate three parameters (magnitude, shape and colour) that can 
be subsequently linearly combined to determine the luminosity 
distances of the SNe la. Despite a common goal, the two meth- 
ods are significantly different in their implementation, allowing 
us to estimate, through the comparison of their outcome, the 



uncertainties associated with the choice of modeling. Whereas 
SALT2 accounts for the SNe variability with a principal compo- 
nent estimation altered by a colour variation law, SiFTO assumes 
a pure time stretching of an average spectral sequence, and relies 
on a limited set of linear relations between rest-frame colours 
of SNe la. The two techniques lead to systematically different 
results (up to 0.02-0.03 on distance moduli), at a level compa- 
rable to other sources of uncertainties such as calibration. The 
quantitative impact of systematic uncertainties, and their relative 
weight with respect to the purely statistical uncertainty on dis- 
tance moduli, depends on the subsequent usage of the distance 
estimates to constrain cosmological parameters, as they intro- 
duce large correla tions from supernova to supernova. 

In contrast to iKessler et al.l (120091) . we find no convincing 
evidence for an evolution of the colour-luminosity relation with 
redshift (i.e. evolution of /?). This is traced to improvements 
made in the empirical modeling of the SNe la spectral sequence 
and to a more precise estimate of their remaining variability be- 
yond that addressed by the models. Nevertheless, we account for 
a possible evolution of /3 in the evaluation of systematic uncer- 
tainties. 

We obtain Q M = 0.211 + 0.034(stat) + 0.069(sys) when fit- 
ting a flat ACDM cosmology to the SNLS data set only. The 
uncertainty is dominated by the calibration uncertainties (0.048) 
despite a large improvement in the precision of this calibration 
with respect to previous supernovae surveys. However, the im- 
pact of systematic uncertainties depends on the data set consid- 
ered. In C10, we add external supernova samples to this one, and 
constrains further the cosmological models while accounting for 
other sources of systematic uncertainties. 
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Table 9. List of SNLS 3 years type la supernova sample. 
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TT~\ O/ICAAAA 

JD-2450000 


lelescope(s) 


Type 


7 


Ker. 


t(B — V) 


Phot. 


03Dlar 


2:27:14.669 
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TJZ. JZ.tl .z / 


J JWZ 


T^ppIt 


OlMd 


7AS0 -t- 001 fTH 
U. /*+ JU ± U.UUIU 


Pn 1 

rdiu 


00Q 
U.UU7 


griz 


uJ-L* jjii 


14-17-9S "384 

It. 1 / .ZJ.JOt 


4-S9-^7-7 SS 

TJZ. J / . / .JJ 


JJuJ 


T^ppIt 

IVCLIv 


OlMd 


7 1 80 -l 00 1 TH^ 

U. / 10U _t \J.\J\Jl\J 


Pn 1 
rdiu 


00Q 
U.UU7 


griz 


UJl-J JJlv 


1 4-1 6-47 4*^ 


4- c ,9-^V^9 Q6 

TJZ.JJ.JZ.7U 




KVpV 

IVCOIv 


CMIo 

OlM^d 


7^60 -4- 001 fT-h 
u. / Juu _t u.uuiu 


Pn 1 
rdiu 


010 

U.U1U 


gnz 


osn^in 

\JJYJ JJ(J| 


14-91 -4S 441 
It.Zl .t J.tt 1 


TJJ. 1 .t / . / 1 


J JZJ 


vJcIIl-l> 


OlM^d 


S7Q0 -1- 001 

u.j / y\j ± u.uuiu ^n^ 


Rr08 

jz>ruo 


010 
U.U1U 


griz 


osrnir 
uji-j Jjr 


1 4-1 Q-98 748 

It. 17.Z0. / to 


TjZ.J 1 .JJ.tJ 


JJuZ 




OlM^d 


^700 -4- 01 00 (*s\ 
U.J /UU ± U.UIUU \o) 


Pn 1 

rdiu 


008 

u.uuo 


griz 


U J l--' J Ivlll 


14-99-^8 97S 

It.ZZ.JO.Z / J 


-i-S^-4-1 90 

TJJ.t. 1 .ZU 




npm M I^ppV 

vJC111 _ 1M ,IVCLIv 


01> Id K 


Q600 + 01 00 
U.7UUU _t vj.vj ivjyj \o ) 


Rr08 

DlUo 


010 
U.U1U 


griz 


UJl-OKp 


1 4-90-9 Q44 
lt.ZW-Z.7tt 


.CO. 1 6-1 4 Q8 
TJZ. ID. It. 70 




JvcCK,vJclIl-lM 


CMIn 

OlMd 


SSOO -4- 01 00 fSft 
U.OJUU ± U.UIUU ^0^ 


Pn 1 

rdiu 


01 1 
U.Ul 1 


griz 


05D3kt 


14:19:53.720 


+52:44:34.75 


3524 


Gem-N 


SNIa 


0.6480 + 0.0010 (H) 


Br08 


0.007 


sriz 


05D3kx 


14:21:50.019 


+53:8:13.63 


3530 


Keck 


SNIa 


0.2190 + 0.0010(H) 


FalO 


0.011 


griz 


05D31a 


14:21:25.418 


+52:21:28.70 


3530 


Keck 


SNIa* 


0.9360 + 0.0010 (H) 


FalO 


0.011 


griz 


05D31b 


14:17:31.758 


+53:10:4.54 


3526 


Gem-N 


SNIa 


0.6470 + 0.0010 (H) 


Br08 


0.011 


griz 


05D31c 


14:22:22.920 


+52:28:43.92 


3527 


Gem-N 


SNIa 


0.4610 + 0.0010(H) 


Br08 


0.009 


griz 


05D31r 


14:22:12.142 


+53:11:3.39 


3530 


Keck 


SNIa 


0.6000 ± 0.0100 (S) 


FalO 


0.011 


griz 



Table 9. continued. 



l\ O.H1C 


R A n9000"> 


DFP H90001 


O JJCLLl Lllll LJalCyoJ 
JL'-z.tJUUuu 


TV* 1 c p r\ r\ f 1 1 c 1 
iClCaLUlJC^o ) 


Type 




l\C 1 . 


F( R — V\ 


Phot 


UJ-Lv Jlllil 


1 4-1 8-SQ S^H 

lH-.iO.J7.OJJ 


4-^9-ZLO-^ 91 

TJZ.4U. J ■ Z 1 


JJU 1 


VJC111-1M 




^700 + 01 00 f*£\ 
u.u /uu _t yj.vjivjyj 


R1-O8 

DIUO 


00Q 

U.UU7 


gnz 


\JJLJ jnm 


1/1.1 8-ZLS 9ZL9 
io.tj.zhz 


4-S9- 1 Q-9^ ^9 

TJZ. 17.ZJ.JZ 


J JO J 


vjclll-l> 




7^00 -+- 001 fTH 
u. /uuu ± \J.\JVJ\-VJ \tn) 


Ur08 

oruo 


00Q 
U.UU7 


griz 


VJJLJ Jllll| 


1 ZL-1 QO ZL98 
14. ly.U.HZo 


_i_S9-9^-6 A'X 
4-JZ.Zj. 0.4J 


J Jo J 


vJClll-lM 


CMIn 


!Af\(\ -t- 001 fTH 
U.Z^fUU ± \J.\J\Ji-VJ 


Ur08 

oruo 


00Q 
U.UU7 


griz 


UJ J--' J 111 A 


1 ZL-99-Q 101 

lt.ZZ.7. 1 U 1 


4-S9- 1 VQ 0^ 

TJZ. 1 J.7.UJ 


JJO 1 


vJcill _ l> 


ol> la 


4700 + 01 00 

U.t- / UU _t U.ULUU \^JJ 


Ur08 
0IU0 


019 
U.U1Z 


gnz 




1 ZL-9 1 -9 Q87 
14. Z 1 .Z.70 ! 


4-S9-9Q-zn S8 
TjZ.Z7.4jJ0 




Oclll-i> 


olMcl 


1 ^QO -1- 001 fV\\ 


oruo 


01 1 
U.U1 1 


gnz 


UJl^T-al 


99-1 6-^ 1 ZL1 
ZZ. ID. J J. 141 


1 8-0-17 79 
-lo.U.l / . / z 


JJJZ 


VT T 

V L 1 


OlMcl 


AQQCi -+- 001 fTH 
U.'H-WU ± U.UULU 


r>aU7 


09Q 
U.UZ7 


griz 


UJl-^ag 


99-1 ^-ZL7 ZL9S 
ZZ. Ij.t/ .^tzj 


1 8-Q-SZL QS 

-10.7.J4.7J 




VI T 

V L 1 




^AOO -4- 01 00 
U.ULUU ± \J.\J1-\J\J \o ) 


r>aU7 


09Q 
U.UZ7 


griz 


UJlyT-av 


99-1ZL-10 SzLzl 

ZZ. 14. lu.J'rT 


1 7-SzL-ZH 00 

-1 / .J4-.4J.UU 


JJU 1 


vJClll-o 




SOQO -4- 00 1 

U.JU"U _t \J.\J\Jl\J 


RrOS 

DIUO 


098 
u.uzo 


griz 


UJlyT-ay 


99-1 zL-^n 1 Q1 

ZZ.lt.JJ.l7l 


1 7-/L6-^ 69 

1 / .4U.J.UZ 


j> J J J 


VI T 

V LI 


Oli la k 


zLOSO -+- 001 (H\ 

U.H-UOU -t \J.\J\JL\J 


RnOQ 
DaU7 


096 
u.uzo 


gnz 


UJl-JH-Dc 


99-16-S^ ^88 
ZZ. ID.jj.joo 


17- 1/1. Q OA 
- 1 / . 14.7.70 


jjjO 


VT T 

V LI 


OlMd 


S^70 -1- 001 
U.JJ /U ± \J.\J\J1-\J 


RnOQ 
r>aU7 


098 
U.UZO 


griz 


UJ L'H-Ul 


99-16-96 ^zL6 

ZZ. 1U.ZU.JHU 


1 8-1 'VSS zlzL 
10. 1j.jj.44 


J> JO 1 


vJclll-o 


CIVTn 


SSQO -+- 001 

U.J0"U _t \J.\J\JL\J 




09Q 

U.UZ7 


gnz 


UJ-LvH-Ul 


99-1 S-S6 S^O 
zz. ij.ju.jju 


1 7-SQ-Q 8S 
1 / .J7.7.0J 




VI T 

V L 1 


olMa 


77S0 -+- 001 (TH 

U. / / JU ± U.UULU 


RnOQ 
DaU7 


098 
u.uzo 


griz 


\JJLJ'+Dj 


99-1ZL-SQ 66ZL 
ZZ. IH.jy.OuH 


1 8-8-0 96 
-lO.O.U.ZU 


JJO 1 


VI T 

V L 1 


olMa 


701 -1- 001 fTH 


RnOQ 
r>aU7 


098 
U.UZO 


griz 


UJJ-Vt-DIII 


ZZ. 1 / .4.0 1 1 


1 7-ZLO-^Q ZLZL 
-1 / .4-U.J7.44- 


J JDZ 


Oclll-o 


olMa 


^790 -t- 001 fTH 
u.j /zu ± \j.\jvj\-\j \tn) 


RrOS 
oruo 


097 
u.uz / 


griz 


UJ L'T-t^ll 


99-1 ^-^ 1 ZL^ZL 

ZZ. 1J.J1 .4J4 


17-17-1Q Szl 
1 / . 1 / . 1 7. 04 


J JOH-, J Jy l 


Ham M VI T 

vJClll _ l> , V LI 


CIVTn 


7^^0 -1- 001 (Th 
u. / u ju _t i\j \n ) 


RrOS RnOQ 

DlUO,DaU7 


096 
u.uzo 


griz 


UJ L'T-t^tJ 


99-1/L-Q 6S0 
ZZ. l^j.OJU 


1 8-1 ^-^6 1 7 

10. 1 J. JU. 1 / 


J Jo J 


VI T 

V L 1 


CMTn 


701 -1- 001 (W\ 

U. /ULU _t U.UULU 


RnOQ 
DaU7 


0^1 

U.UJ 1 


gnz 


UJl-JH-CS 


99. 17. 17 818 
zz. 1/.1/.0I0 


1 7-S9-/L8 SZL 
-1 / .JZ.4o.j4 


J JO J 


VT T 

V Li 


OlMd 


7Q00 -4- 01 00 

U. / VUU ± U.ULUU \o) 


RnOQ 
r>aU7 


098 
U.UZO 


griz 


UJ-LvH-L. W 


99-1ZL-S0 080 
zz.14.ju.uou 


1 7-ZLZL- 1 Q S6 

1 / .44. 1 7. JO 


JJO" 


VT T 

V LI 


olMa 


^7S0 -1- 001 


RnOQ 
DaU7 


096 
u.uzo 


griz 


UJ-LJt-UL 


99-1ZL-9S 8ZL7 
zz. i^.zj.o^t / 


1 7-ZL0- 1 6 00 
-1 / .4U. ID. UU 


JD 1 J 


Oclll-IM 


olMa 


A070 -1- 001 fTH 
u.^fu/u ± u.uulu \tn) 


RrOS 

oruo 


096 
U.UZO 


griz 




99.1 6-SS 878 
ZZ. 1O.JJ.0 / 


1 8-^-S 96 
-lu.J.J.ZO 


JOZZ 


VT T 

V Li 


olMa 


8 S SO -1- 001 ^T4^ 
u.ojju ± \j.\jvj\-VJ \tn) 


RnOQ 
r>aU7 


098 
U.UZO 


griz 


OSDzLHv 


99-1 S-^0 1 1 Q 

ZZ. IJ.JU. 117 


1 8-1 81 
- 1 0. 1 Z.J J. 1 


JO 1 " 


vJClll _ l> 




81 00 + 01 00 

U.O LUU _t U.ULUU \^JJ 


RrOS 
DlUo 


096 
u.uzo 


gnz 




99.1 'J-SQ ^0^ 

ZZ. 1J.J7.JUJ 


18-1 9-zLzL ZLZL 
- 1 0. 1 Z. 44. 44 


JQ JQ 


VT T 

V Li 


OlMa* 


f^OSO -1- 001 
U.UUJU ± \J.\J\J1-\J 


RnOQ 
r>aU7 


0^0 
U.UJU 


griz 


UJLJ'H-cJ 


99.1 ^.<9 S 17 
ZZ. 1 J.jZ.J 1 / 


18-11 -ZLZL SQ 

-lo.ll .44.J7 


ju jy 


VT T 

V Li 


OlMa 


S8S0 -i- 001 *rn 

U.JOJU ± U.UULU 


RnOQ 
r>aU7 


09Q 
U.UZ7 


gnz 




99-1 6-97 S1 7 

ZZ. 1 U.Z / . J 1 / 


1 7-zLzL- 1 7ZL 

1 / .44. 1 U. / 4 


JO JO 


VT T 

V L 1 


JlMu 


S^^O -1- 001 (TH 
u.jjuu _t u.uulu 


RnOQ 
DaU7 


098 
u.uzo 


griz 


UJ-L-Jt-cv 


99.1 -3-^0 ZL9ZL 

ZZ. 1 J . JU.4Z4 


1 7-ZL9-9 8^ 
-1 / .4Z.Z.OJ 


JUZ^f 


VI T 

V Li 


OlMa 


7990 -t- 001 ^T4^ 
U. /ZZU ± \J.\JVJ1-VJ 


RnOQ 
r>aU7 


097 
U.UZ / 


griz 


UJ-LJt-Ic 


99.1 S-ZLO 968 
ZZ. Ij.H-U.ZOu 


17-^1 -ZL9 ^8 
- 1 / . J 1 .4Z.JO 


Ju JV 


VI T 

V Li 


OlMax 


QRAO -i- 001 fTH 
u.vo^fu ± vj.vvji-Vj \sn) 


RnOQ 
r>aU7 


09ZL 
U.UZ4 


griz 


UJJ-^H-ll 


99-1 6-90 1 60 
zz. 10. zu.iuu 


1 8-9-^ 90 
-lo.z.jj.zu 


JUt-Z 


VI T 

V LI 


oi> la k 


4070 -1- 001 (Th 
u.t-uzu -t u.uulu 


RnOQ 
DaU7 


09Q 

U.UZ7 


griz 




99-16-ZL1 ^zLS 
ZZ. 10.41 . jto 


1 7-^^-zL^ f\f\ 
-11 .JJ.4J.00 


JU'H-l 


VT T 

V Li 


olMa 


S^QO -i- 001 fV\\ 
U.OJVU ± \J.\J\J1-\J 


RnOQ 
r>aU7 


09Q 
U.UZ7 


gnz 




99.1 S-90 Q1 9 
ZZ. lj.ZW.7lZ 


1 7. 1 f..c 07 
-1 / .ID.J.U/ 


JUJW 


Ocill-o 


OlMa 


^7^0 -1- 001 
U.J / JU ± \J.\J\J \.\J yn ) 


RrOS 
oruo 


096 
U.UZO 


griz 




99-1ZL-ZL7 ^S6 
ZZ.14.4/.JJU 


1 "7.-3 1 -SzL 8ZL 

1 / . J 1 . J4.04 


JU / " 


VJClll-lM 


OlMa 


RORO -i- 001 (TH 

U.OUOU _t U.UULU 


RrOS 

D1UO 


09ZL 

U.UZ4 


gnz 


UJ-LJt-1111 


99.17.1 "2 

ZZ. 1 / . 1 J.JD4 


1 7-SZL-ZLS 7ZL 
-1 / .J4.4J. /4 


J / UJ 


Oclll-o 


OlMa* 


8490 -t- 001 fTH 
u.o^fzu ± \j.\j\jL\j \tn) 


RrOS 

oruo 


09Q 
U.UZ7 


griz 


\JDLJ 1 aD 


9-9ZL- SQ 0^1 

Z.Z4.J7.UJ 1 


ZL-ZL0-S1 90 
-4.4U.J 1 .ZU 


J /^Hl 


VI T 

V Li 


OlMa 


1 890 -1- 001 fTH 

\J.\.OZ.\J ± \J.\J\J1.\J 


RnOQ 
r>aU7 


096 
U.UZO 


griz 


UUlVZdg 


1 0-1 -zH ^6^ 

1U.1 .TJ.JUJ 


-1-1 -S1 -^7 07 

Tl.Jl.J/.U/ 


J / u 1 


VI T 

VLl 


OlMa 


^ 1 00 -1- 00 1 (H\ 

U.JLUU _t \J.\J\JL\J 


RnOQ 
DaU7 


01 Q 
U.UI7 


gnz (,s; 


UULJZDK 


Q-^8-zL9 87 1 
7.J0.4Z.0 1 1 


4-9-1 0-1 Q 1zL 

4-Z. Iu.l7.14 


^77A 


VT T 

V Li 


OlMa 


4QQ0 -1- 001 


RnOQ 
r>aU7 


01 Q 
U.U1 7 


griz 




Q-^Q-Q SS7 

7.J7.7.JJ / 


4-9-1 8-Q 06 
tZ. I0.7.UO 


J / 7J 


VT T 

V Li 


OlMa 


S^10-i-0 0010 fV\\ 
U.JJLU ± \J.\J\J1-\J 


RnOQ 
r>aU7 


090 

u.uzu 


gnz [s) 


UUL'Z,L U 


1 0-1 -9^ ZLQ^ 

1U.1 .ZJ.47J 


4-9-97-99 

TZ.Z / .ZZ. JJ 


J / "J 


VI T 

V LI 


oi> la k 


1 0000 + 01 00 
l .\j\J\J\J _t vj.vj ivjyj v o ) 


RnOQ 
DaU7 


01 Q 
U.UI7 


griz 


UU-LJZLC 


Q-SQ-96 887 
7.J7.ZO.00 / 


4-9-97-ZL 88 
+z.z / .4.00 


j / y j 


VI T 

V Li 


CMIn 

OlMa 


S^90 -t- 001 ^T4^ 
U.JJZU ± \J.\JVJ\-VJ \tn) 


RnOQ 
r>aU7 


01 Q 
U.U1 7 


griz 


06D2cd 


10:0:57.961 


+2:32:2.54 


3797 


VLT 


SNIa 


0.9300 + 0.0010 (H) 


Ba09 


0.020 


sriz (s) 


06D2ce 


10:1:44.219 


4-2:42:47.93 


3796 


VLT 


SNIa 


0.8200 + 0.0100 (S) 


Ba09 


0.020 


griz (s) 


06D2ck 


10:1:28.016 


+ 1:51:49.07 


3802 


VLT 


SNIa* 


0.5520 + 0.0010 (H) 


Ba09 


0.019 


griz (s) 


06D2ez 


10:1:38.665 


+ 1:56:38.72 


3853 


Keck 


SNIa 


0.0820 ± 0.0010 (H) 


FalO 


0.018 


griz (s) 


06D2fb 


9:59:32.147 


+ 1:49:40.52 


3855,3854 


Keck 


SNIa 


0.1240 + 0.0010(H) 


FalO 


0.021 


griz (s) 


06D2ff 


10:0:16.345 


+2:14:35.14 


3853 


Keck 


SNIa 


0.3450 ± 0.0010 (H) 


FalO 


0.019 


griz (s) 



Table 9. continued. 



1 1 tllllC 


R A H9000) 


DFP H70001 


O JJCLLl Hill LJalCyoJ 

m-74S0000 

J i-J Z.'-r J UU UU 


Tip 1 c c r\ r\ f 1 1 c 1 


Type 




l\C 1 . 


F( R — V\ 


Phnt 




7.J7.0. iu~ 


+ 1 -SfV1 8 S8 
T 1 .Ju. 1 0.J0 


jOOI 


VI T 
V LI 


Ol> IdT 


U.Ot-UU It U.U1UU ^0 ) 


RnOQ 


070 

U.UZU 


gnz 




1 0-0-ZL7 zL1 Q 
1U.U.+ / .^4-17 


-1-1 -SQ-7 ZL6 

T 1 .J7. / .+0 


^88zL 






U.T-fZU it W.WW1U 


1 


001s 


griz 


\J\jLJ jDZ 


1/1.1 7. 1 OA? 
1*+. 1 / . IU.U^tZ 


4-S3- 1 -9Q 37 
Tjj. 1 .Ay .j 1 


j l l L r 


vJcIIl-lM 




7970 -t- 001 


Ur08 


00Q 

v.wy 


griz 




1 ZL-90-ZL3 S7zL 


-1-S9- 1 1 -97 83 

TJZ. 1 1 .Z / .OJ 


^76Q 
J 1 DV 


r^^-m 1ST 


olMa* 


1 0^70 -1- 001 fV\\ 


r>rwo 


0010 


griz 


uulj jcc 


izi.17.a1 ^nn 


-1-S9-SZL-ZLZL QO 
tjz. j i +.'+'+. yv 


^768 
J / Do 


r^^m 1ST 
vJCill-lN 


CM|,,4. 
olMa* 


^8^0 -1- 001 (W\ 
u.Doju it U.UU1U 


RrOS 


008 

u.uuo 


gnz 


UUl>OCIl 


1 /I.1 Q-9C C/IC 


4-S9-38-97 SzL 

TjZ.JO.Z / .J^t 


^800 


vJClIl-l> 


OlMa 


7^70 -i- 001 (W\ 


RrOS 

orwo 


00Q 

v.wy 


gnz ,s; 


UOLV JU1 


1/1 .99. 1 7 900 
l^t.ZZ. 1 / .LVy 


4-S9-S7-97 6S 
TJZ.J / .Z / .OJ 






OlMa 


AA90 -t- 001 (W\ 
U.T-fZU it W.WW1U 


1 


0010 

U.U1U 


griz 




1/1 .99. 1 9 C77 
l^t.ZZ. 1Z.J / / 


4-S9-38-97 OS 

TJZ.JO.Z / .UJ 


jOjj 




OlMclK 


^S70 -i- 001 


Pn 1 
FtllU 


00Q 

v.wy 


gnz ^s; 


UUlvJUU 


1 zL-1 fv^Q 3Q? 


4_c9.c9.c7 1£ 

TJZ. JZ.J / . 1U 




IVCOIv 


OlMa 


7760 -i- 001 (H\ 
\J. 1 ZOW it U.UUIU 


Pn 1 


010 
V.V LV 


gnz 


UOjLOUL 


1/1- 17-91 fifvd 
It. 1 / .z 1 .OOt 


. c9.97.1n 1zL 
TjZ.Z / .iu.ih 


^860 ^oca 




OlMa 


7870 -1- 001 fV\\ 


Pn 1 
rdiu 


008 
U.UUO 


griz 


o^rnpn 

UUL' JCLl 


1 zL- 1 8-ZL7 80? 


4-S7-ZLS-93 61 

TjZ.tJ.Zj.U 1 


^881 


rjpm N 


OlMa 


dOAO -1- 001 

W.H-Wt-W it U.UUIU 


W1 

VV 1U 


008 
u.uuo 


gnz 




14-17-1 087 


-t-S7- 1 3S6 6zL 
tjz. IJ.JO.OH 


JOOI ,jOOJ 


vJC111 _ 1M ,JVCOIv 


OlMa 


S 1 QO -1- 00 1 
U.JI7U it \j.\j\Ji\j 


W1 

VV 1U 


013 

U.Ul J 


griz 


UO-LJ JcIIl 


1 zL- 1 Q-93 zL1 Q 
L^+. iy .Aj iy 


1 ca. 1 .99 oc 

TJJ. 1 .ZZ. OJ 






OlMa 


^QOO -+- 01 00 fSft 
U.DVUU ± U.U1UU ^0 ^ 


Pn 1 
FtllU 


01 1 
U.Ul 1 


griz 


WDLV Jell 


1/1-91 • 1 3 


. c9.97.91 00 
Tjz.z /.zi.oo 






OlMclK 


1 0/^00 -+- 01 00 f<5^ 


Pn 1 
FdlU 


01 1 

U.Ul 1 


griz 


o^rnpt 


1 zL-??-zL 770 


4-S7-70-1 8 1 7 

TJZ.ZU.IO.IZ 


JOO\J 




OlMa 


S7fS0 -1- 001 (W\ 

U.J / UU it U.UUIU 
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06D4bw 
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VLT 


SNIa 


0.7320 + 0.0010 (H) 


Ba09 


0.025 


2riz 


06D4ce 


22:17:17.245 


-17:14:9.37 


3935 


VLT 


SNIa 


0.8500 + 0.0100 (S) 


Ba09 


0.028 


griz 


06D4cl 


22:16:57.053 


-17:27:2.28 


3937 


VLT 


SNIa 


1.0000 + 0.0100 (S) 


Ba09 


0.027 


griz 


06D4co 


22:15:26.518 


-17:52:9.07 


3946 


VLT 


SNIa 


0.4370 + 0.0010 (H) 


Ba09 


0.026 


griz 


06D4cq 


22:16:55.523 


-17:42:43.13 


3946 


VLT 


SNIa 


0.41 10 + 0.0010(H) 


Ba09 


0.027 


griz 


06D4dh 


22:14:19.693 


-17:35:6.03 


3946 


Gem-N 


SNIa 


0.3030 + 0.0010 (H) 
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Notes. Photometric status, (a) : no photometry because of variable host or second SN in host, (b) close to a bright star, (c) edge of field; (s): missing epochs for light curve fit. 



Refe rences. H05:lHowell et alj (12005?). H06: lHowell et all (2006). Ba09: iBalland et al.1 d2009t) . Br08: iBronder et al.1 dzOul) E08: lEllis et al.1 ( 120081) . F10: lFakhouri & SNLS collaboration! JUT 
W10: IWalker et all din prep.1) . 
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Table 10. Photometry of the SNLS 3-year SNe la. 



Name 


Filter" 


MJD 


Flux" 


Flux uncertainty 6 


(BDlar 


g 


52880.58 


1.3035E+01 


5.4246E+00 


(BDlar 


i 


52881.50 


1.6157E+01 


1.8629E+01 


(BDlar 


r 


52881.54 


7.0282E+00 


9.9105E+00 


(BDlar 


z 


52881.56 


2.2086E+01 


7.1951E+01 


03Dlar 


i 


52886.60 


1.8944E+02 


1.9855E+01 


03Dlar 


g 


52900.50 


1.4280E+02 


1.4970E+01 


03Dlar 


i 


52900.53 


6.5402E+02 


2.3202E+01 


(BDlar 


r 


52900.56 


4.5705E+02 


1.6919E+01 


(BDlar 


z 


52900.59 


5.7204E+02 


8.9594E+01 


(BDlar 


i 


52904.53 


6.941 1E+02 


1.7710E+01 


(BDlar 


r 


52904.57 


3.9937E+02 


1.2907E+01 


(BDlar 


g 


52904.60 


U473E+02 


7.5546E+00 


(BDlar 


g 


52908.53 


8.4923E+01 


7.5163E+00 


(BDlar 


r 


52908.56 


3.8457E+02 


1.2797E+01 


(BDlar 


i 


52908.60 


6.0466E+02 


2.0464E+01 


(BDlar 


z 


52909.61 


6.6974E+02 


8.0965E+01 


(BDlar 


i 


52912.49 


5.7083E+02 


1.9573E+01 


(BDlar 


r 


52914.57 


2.0934E+02 


2.1674E+01 


(BDlar 


i 


52915.62 


5.0601E+02 


2.5140E+01 


(BDlar 


r 


52916.46 


2.7066E+02 


1.2909E+01 



Notes. An electronic version of this table is available at the Centre de Donnees astronomiques de Strasbourg (CDS). Light curves can also be 
downloaded at the University of Toronto's Research Repository https : //tspace . library .utoronto . ca/handle/ 1807/245 12. 
(n) The average focal plane coordinates of SNe are available at CDS. Those are needed to estimate the filter transmission function, see § 13.21 
ib) Fluxes and uncertainties are given for a fiducial zero point of 30, a magnitude is mag = -2.5 log 10 (Flux) + 30. The interpretation of those 
magnitudes for supernovae study is detailed in § 13.31 
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Table 11. SNLS Type la supernovae parameters. 



Name 


m*(SALT2) 


Xi(SALT2) 


C(SALT2) 


m*(SIFTO) 




s(SIFTO) 


C(SIFTO) 


(BDlar 


23 


941 


+ 





033 


-0.945 ± 0.209 


0.266 ± 0.035 


23 


969 


+ 





036 





878 


± 





026 


0.289 ± 0.024 


03Dlau 


23 


009 


+ 





026 


1.273 ±0.198 


-0.015 ± 0.030 


23 


004 


+ 





031 


1 


144 


± 





019 


0.015 ± 0.025 


03Dlaw 


23 


584 


+ 





030 


1.161 ±0.342 


-0.021 ±0.037 


23 


562 


+ 





040 


1 


065 


± 





030 


-0.015 ± 0.042 


03Dlax 


22 


961 


+ 





026 


-0.574 ±0.1 14 


-0.112 ±0.031 


23 


011 


+ 





032 





916 


± 





013 


-0.048 ± 0.025 


03Dlbk 


24 


296 


+ 





047 


0.211 ±0.274 


-0.249 ± 0.054 


24 


363 


+ 





067 


1 


016 


± 





030 


-0.155 ±0.044 


03Dlbm 


24 


335 


+ 





039 


0.567 ± 0.347 


0.358 ± 0.056 


24 


330 


+ 





047 


1 


061 


± 





042 


0.350 ± 0.052 


03Dlbp 


22 


387 


+ 





032 


-1.145 ±0.135 


-0.051 ±0.028 


22 


453 


+ 





030 





893 


± 





013 


0.030 ±0.019 


03Dlco 


24 


076 


+ 





046 


0.833 ± 0.445 


-0.047 ± 0.064 


24 


075 


+ 





049 


1 


082 


± 





031 


-0.017 ±0.041 


03Dlew 


24 


323 


+ 





050 


0.502 ± 0.365 


-0.060 ± 0.064 


24 


362 


+ 





074 


1 


023 


± 





037 


-0.056 ± 0.048 


03Dlfc 


21 


855 


+ 





030 


0.599 ±0.172 


-0.026 ± 0.025 


21 


867 


+ 





028 


1 


037 


± 





019 


0.025 ±0.018 


03Dlfq 


24 


503 


+ 





045 


-0.752 ± 0.366 


-0.061 ±0.063 


24 


493 


+ 





059 





828 


± 





039 


-0.075 ± 0.045 


03Dlgt 


24 


042 


+ 





065 


-1.032 ±0.598 


0.333 ± 0.084 
























03D3af 


23 


497 


+ 





039 


0.046 ± 0.254 


0.068 ± 0.069 
























03D3aw 


22 


631 


+ 





045 


0.620 ±0.281 


-0.091 ±0.043 


22 


664 


+ 





046 


1 


068 


± 





036 


-0.054 ± 0.038 


03D3ay 


22 


269 


+ 





044 


0.511 ±0.230 


-0.058 ± 0.035 


22 


310 


+ 





045 


1 


058 


± 





032 


-0.008 ± 0.036 


03D3ba 


21 


959 


+ 





042 


0.891 ±0.172 


0.101 ±0.039 


22 


018 


+ 





042 


1 


112 


± 





040 


0.164 ±0.025 


03D3bl 


22 


932 


+ 





028 


0.266 ± 0.236 


0.210 ±0.031 


22 


958 


+ 





028 





972 


± 





036 


0.237 ±0.019 


03D3cd 


22 


583 


+ 





070 


1.882 ±0.828 


-0.037 ± 0.049 


22 


587 


+ 





037 


1 


211 


+ 





058 


0.023 ± 0.025 


03D4ag 


21 


273 


+ 





033 


0.938 ±0.149 


-0.070 ± 0.024 


21 


264 


+ 





026 


1 


115 


± 





018 


-0.049 ± 0.017 


03D4at 


23 


732 


+ 





033 


0.203 ± 0.350 


-0.058 ± 0.065 


23 


713 


+ 





046 


1 


014 


± 





044 


0.012 ± 0.042 


03D4au 


23 


789 


+ 





033 


0.481 ±0.346 


0.120 ±0.044 


23 


832 


+ 





029 


1 


054 


± 





039 


0.160 ±0.036 


03D4cj 


21 


054 


+ 





032 


1.240 ±0.134 


-0.088 ± 0.025 


21 


038 


+ 





023 


1 


114 


± 





014 


-0.069 ± 0.014 


03D4cx 


24 


439 


+ 





062 


-0.108 ±0.582 


0.017 ± 0.063 


24 


433 


+ 





083 





905 


± 





084 


-0.011 ±0.057 


03D4cy 


24 


669 


+ 





081 


0.857 ±0.686 


-0.073 ± 0.068 


24 


702 


+ 





099 


1 


119 


± 





066 


-0.017 ±0.056 


03D4cz 


24 


016 


+ 





051 


-1.800 ±0.433 


-0.147 ±0.089 


24 


049 


+ 





059 





837 


± 





037 


-0.038 ± 0.047 


03D4dh 


23 


387 


+ 





030 


1.277 ±0.234 


-0.074 ± 0.048 


23 


370 


+ 





041 


1 


121 


± 





017 


-0.047 ± 0.033 


03D4di 


24 


296 


+ 





056 


1.347 ±0.450 


-0.089 ± 0.060 


24 


298 


+ 





074 


1 


151 


± 





040 


-0.033 ± 0.047 


03D4dy 


23 


255 


+ 





032 


1.354 ±0.250 


-0.094 ± 0.036 


23 


263 


+ 





042 


1 


112 


± 





018 


-0.055 ±0.031 


03D4fd 


24 


222 


+ 





041 


1.019 ±0.489 


-0.052 ± 0.062 


24 


215 


+ 





053 


1 


045 


± 





049 


-0.033 ± 0.045 


03D4gf 


23 


325 


+ 





028 


0.719 ±0.356 


-0.034 ± 0.039 


23 


329 


+ 





038 


1 


120 


± 





033 


0.022 ± 0.042 


03D4gg 


23 


422 


+ 





036 


0.906 ± 0.474 


0.007 ± 0.042 


23 


384 


+ 





044 


1 


098 


± 





052 


0.053 ± 0.047 


04Dlag 


22 


982 


+ 





032 


0.316 ±0.345 


-0.165 ±0.050 
























04Dlaj 


23 


897 


+ 





034 


0.507 ± 0.323 


-0.038 ± 0.062 


23 


884 


+ 





048 


1 


010 


± 





051 


-0.006 ± 0.045 


04Dlak 


23 


641 


+ 





048 


-1.678 ±0.350 


0.015 ± 0.055 


23 


611 


+ 





054 





815 


± 





030 


-0.155 ±0.100 


04Dldc 


21 


072 


+ 





030 


-1.137 ±0.087 


-0.000 ± 0.026 


21 


084 


+ 





018 





833 


± 





031 


0.011 ±0.014 


04Dlde 


24 


124 


+ 





035 


0.889 ± 0.262 


-0.132 ± 0.052 


24 


139 


+ 





047 


1 


111 


± 





025 


-0.064 ± 0.038 


04Dlff 


24 


219 


+ 





042 


0.738 ±0.301 


0.029 ± 0.055 


24 


237 


+ 





057 


1 


087 


± 





029 


0.037 ± 0.043 


04Dlhd 


22 


162 


+ 





025 


0.756 ± 0.096 


-0.096 ± 0.023 


22 


159 


+ 





026 


1 


065 


± 





012 


-0.059 ±0.015 


04Dlhx 


23 


703 


+ 





026 


0.474 ± 0.206 


0.107 ±0.035 


23 


715 


+ 





037 


1 


038 


± 





024 


0.149 ±0.030 


04Dlhy 


24 


299 


+ 





041 


1.324 ±0.325 


-0.053 ± 0.053 


24 


284 


+ 





054 


1 


086 


± 





027 


-0.027 ± 0.043 


04Dliv 


24 


600 


+ 





059 


1.510 ±0.418 


-0.121 ±0.052 


24 


583 


+ 





079 


1 


140 


± 





032 


-0.054 ± 0.047 


04Dljd 


24 


407 


+ 





038 


0.250 ± 0.299 


0.060 ± 0.065 


24 


419 


+ 





051 





999 


± 





030 


0.005 ± 0.043 


04Dljg 


23 


264 


+ 





028 


0.337 ±0.182 


-0.1 13 ±0.035 


23 


263 


+ 





039 


1 


019 


± 





023 


-0.079 ± 0.039 


04Dlkj 


23 


341 


+ 





025 


0.306 ±0.145 


-0.080 ± 0.029 


23 


331 


± 





037 


1 


033 


± 





021 


-0.054 ± 0.029 


04Dlks 


24 


130 


+ 





038 


0.874 ± 0.287 


0.080 ± 0.058 


24 


132 


± 





050 


1 


069 


± 





028 


0.049 ± 0.043 


04Dloh 


23 


385 


+ 





032 


0.011 ±0.240 


-0.075 ± 0.033 


23 


376 


+ 





041 


1 


020 


± 





031 


-0.050 ± 0.036 


04Dlow 


24 


317 


+ 





046 


-0.041 ± 0.334 


-0.180 ±0.051 


24 


376 


+ 





064 


1 


017 


± 





032 


-0.125 ±0.043 


04Dlpc 


24 


545 


+ 





036 


0.014 ±0.340 


0.069 ± 0.060 


24 


537 


+ 





051 





944 


± 





032 


0.036 ± 0.041 


04Dlpd 


24 


681 


+ 





059 


0.371 ±0.476 


0.009 ± 0.058 


24 


729 


+ 





077 


1 


034 


± 





049 


0.001 ± 0.049 


04Dlpg 


23 


566 


+ 





033 


1.041 ±0.224 


0.085 ± 0.037 


23 


565 


+ 





037 


1 


076 


± 





024 


0.120 ±0.029 


04Dlpp 


23 


993 


+ 





031 


-1.166 ±0.244 


-0.118 ±0.051 


23 


975 


+ 





047 





861 


± 





024 


-0.062 ± 0.036 


04Dlpu 


23 


988 


+ 





053 


-1.502 ±0.388 


0.074 ± 0.084 


24 


035 


+ 





060 





848 


± 





035 


0.081 ±0.056 


04Dlqd 


24 


239 


+ 





033 


0.337 ±0.272 


-0.002 ±0.051 


24 


191 


+ 





049 





984 


± 





029 


-0.010 ± 0.040 


04Dlrh 


22 


570 


+ 





025 


1.111 ±0.253 


-0.058 ± 0.026 


22 


581 


± 





028 


1 


058 


+ 





025 


-0.011 ±0.016 


04Dlrx 


24 


732 


+ 





060 


0.692 ± 0.500 


-0.107 ±0.059 


24 


745 


+ 





087 


1 


085 


± 





053 


-0.047 ± 0.050 


04Dlsa 


23 


563 


+ 





034 


-0.468 ± 0.307 


-0.092 ± 0.035 


23 


533 


+ 





042 





926 


± 





030 


-0.072 ± 0.044 


04Dlsi 


23 


854 


+ 





032 


0.159 ±0.294 


-0.030 ± 0.052 


23 


852 


+ 





046 





974 


± 





033 


-0.008 ± 0.036 


04Dlsk 


24 


060 


+ 





031 


0.364 ±0.381 


0.075 ± 0.053 


24 


030 


+ 





045 


1 


017 


± 





042 


0.098 ± 0.037 


04D2cc 


24 


236 


+ 





107 


0.133 ±0.788 


-0.059 ± 0.076 


24 


046 


+ 





139 


1 


196 


+ 





063 


-0.126 ±0.111 


04D2cf 


22 


474 


± 





050 


-0.733 ±0.215 


-0.015 ± 0.036 


22 


508 


+ 





069 





835 


± 





063 


0.019 ±0.022 
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Table 11. continued. 



Name 


m*(SALT2) 


Xi(SALT2) 


C(SALT2) 


m*(SIFTO) 




s(SIFTO) 


C(SIFTO) 


04D2fp 


22 


552 


+ 





027 


0.375 + 0.145 


-0.039 ± 0.028 


22 


551 


+ 





030 


1 


033 


+ 





020 


-0.022 ±0.021 


04D2fs 


22 


431 


+ 





027 


0.136 + 0.118 


0.058 ± 0.026 


22 


427 


+ 





028 


1 


002 


+ 





019 


0.068 ±0.019 


04D2gb 


22 


900 


+ 





031 


-1.527 ±0.202 


0.017 ±0.038 


22 


916 


+ 





033 





830 


+ 





019 


0.028 ± 0.027 


04D2gc 


23 


322 


+ 





031 


1.193 ±0.271 


0.012 ±0.037 


23 


322 


+ 





034 


1 


143 


+ 





024 


0.026 ±0.031 


04D2gp 


24 


207 


+ 





052 


-1.261 ±0.436 


-0.150 ±0.082 


24 


200 


+ 





060 





837 


+ 





035 


-0.090 ± 0.047 


04D2iu 


24 


233 


+ 





052 


-1.757 ±0.459 


0.005 ± 0.084 


24 


233 


+ 





055 





805 


+ 





034 


-0.003 ± 0.045 


04D2ja 


24 


114 


+ 





053 


0.530 ±0.421 


-0.213 ±0.066 


24 


114 


+ 





058 


1 


026 


+ 





045 


-0.104 ±0.043 


04D2kr 


23 


863 


+ 





032 


0.875 ± 0.306 


-0.046 ± 0.050 


23 


841 


± 





047 


1 


026 


+ 





038 


-0.027 ± 0.038 


04D2mc 


22 


517 


+ 





039 


-1.650 ±0.149 


0.085 ± 0.033 


22 


629 


± 





035 





869 


+ 





017 


0.166 ±0.022 


04D2mh 


23 


383 


+ 





031 


1.599 ±0.235 


-0.004 ± 0.035 


23 


405 


+ 





039 


1 


150 


+ 





017 


0.056 ± 0.027 


04D2mj 


23 


764 


+ 





033 


1.521 ±0.263 


0.127 ±0.037 


23 


791 


+ 





037 


1 


148 


+ 





022 


0.181 ±0.025 


04D3co 


23 


761 


+ 





036 


-0.664 ± 0.234 


-0.044 ± 0.049 


23 


734 


+ 





050 





943 


+ 





024 


0.045 ± 0.054 


04D3cp 


24 


093 


+ 





042 


0.714 ±0.293 


-0.255 ± 0.055 


24 


095 


+ 





057 


1 


027 


+ 





028 


-0.164 ±0.042 


04D3cy 


23 


779 


+ 





036 


-0.087 ± 0.247 


-0.029 ± 0.062 


23 


799 


+ 





046 





964 


+ 





025 


-0.026 ± 0.041 


04D3dd 


25 


016 


+ 





096 


0.678 ± 0.733 


-0.034 ± 0.072 


25 


366 


+ 





184 


1 


153 


+ 





074 


-0.036 ± 0.069 


04D3df 


23 


532 


+ 





029 


-1.780 ±0.173 


0.110 ±0.034 


23 


499 


+ 





030 





770 


+ 





015 


0.072 ± 0.022 


04D3do 


23 


573 


+ 





028 


-0.723 ±0.172 


-0.120 ± 0.034 


23 


561 


+ 





041 





892 


+ 





019 


-0.119 ±0.034 


04D3ez 


21 


668 


+ 





029 


-0.975 ±0.1 12 


0.033 ± 0.026 


21 


715 


± 





019 





887 


+ 





021 


0.111 ±0.013 


04D3fk 


22 


508 


+ 





026 


-0.364 ± 0.093 


0.067 ± 0.024 


22 


554 


+ 





026 





960 


+ 





015 


0.115 ±0.017 


04D3fq 


24 


105 


+ 





037 


-0.143 ±0.353 


-0.054 ± 0.064 


24 


101 


+ 





049 





935 


+ 





042 


-0.046 ± 0.042 


04D3gt 


23 


226 


+ 





026 


-0.182 ±0.122 


0.177 ±0.029 


23 


279 


+ 





015 





980 


+ 





017 


0.234 ±0.019 


04D3gx 


24 


620 


+ 





060 


-0.474 ± 0.401 


-0.159 ± 0.057 


24 


672 


+ 





080 





933 


+ 





047 


-0.124 ±0.048 


04D3hn 


23 


484 


+ 





028 


-0.460 ± 0.144 


0.067 ± 0.037 


23 


512 


± 





036 





923 


+ 





017 


0.093 ±0.031 


04D3kr 


21 


941 


+ 





027 


1.155 ±0.122 


-0.042 ± 0.024 


21 


961 


+ 





026 


1 


130 


+ 





011 


0.003 ±0.016 


04D3ks 


23 


818 


+ 





043 


0.610 ±0.311 


-0.074 ± 0.060 


23 


865 


+ 





052 


1 


089 


+ 





031 


-0.017 ± 0.040 


04D31p 


24 


978 


+ 





103 


-1.491 ±0.517 


-0.066 ± 0.077 


25 


003 


+ 





151 





799 


+ 





061 


-0.079 ± 0.063 


04D31u 


24 


357 


+ 





039 


-0.459 ± 0.248 


-0.123 ± 0.057 


24 


361 


+ 





054 





923 


+ 





026 


-0.104 ±0.042 


04D3mk 


24 


275 


+ 





038 


-0.201 ± 0.264 


-0.143 ± 0.054 


24 


280 


+ 





052 





934 


+ 





032 


-0.116 ±0.042 


04D3ml 


24 


542 


+ 





058 


1.038 ±0.426 


-0.128 ±0.053 


24 


526 


+ 





074 


1 


124 


+ 





045 


-0.049 ± 0.047 


04D3nc 


24 


275 


+ 





048 


1.455 ±0.579 


-0.058 ± 0.058 


24 


280 


+ 





054 


1 


131 


+ 





039 


-0.025 ± 0.043 


04D3nh 


22 


135 


+ 





025 


0.595 ±0.1 17 


-0.024 ± 0.023 


22 


137 


± 





025 


1 


052 


+ 





017 


0.008 ±0.016 


04D3nq 


21 


000 


+ 





028 


0.620 ±0.108 


-0.015 ± 0.026 


21 


007 


± 





021 


1 


054 


+ 





022 


0.015 ±0.016 


04D3nr 


24 


574 


+ 





066 


0.312 ±0.492 


-0.007 ± 0.055 


24 


549 


+ 





076 





950 


+ 





050 


-0.011 ±0.047 


04D3ny 


24 


221 


+ 





054 


-0.071 ±0.618 


-0.025 ± 0.066 


24 


272 


+ 





064 


1 


061 


+ 





054 


0.006 ± 0.045 


04D3oe 


24 


060 


+ 





033 


-0.771 ±0.312 


-0.247 ± 0.057 


24 


035 


+ 





046 





906 


+ 





027 


-0.136 ±0.039 


04D4an 


24 


027 


+ 





042 


0.074 ± 0.435 


0.012 ± 0.043 


24 


044 


± 





053 





967 


+ 





063 


0.018 ±0.042 


04D4bk 


24 


312 


+ 





046 


1.841 ±0.498 


-0.094 ± 0.060 


24 


309 


± 





058 


1 


168 


+ 





040 


-0.057 ± 0.045 


04D4bq 


23 


339 


+ 





034 


0.823 ±0.282 


0.100 ±0.040 


23 


342 


+ 





040 


1 


101 


+ 





027 


0.138 ±0.036 


04D4dm 


24 


384 


+ 





044 


-0.176 ±0.337 


-0.006 ± 0.060 


24 


388 


+ 





059 





985 


+ 





035 


0.005 ± 0.045 


04D4dw 


24 


551 


+ 





076 


1.129 ±0.474 


-0.143 ± 0.063 


24 


482 


+ 





114 


1 


131 


+ 





043 


-0.038 ± 0.056 


04D4fx 


23 


506 


+ 





033 


1.021 ±0.243 


-0.041 ± 0.049 


23 


472 


+ 





042 


1 


116 


+ 





020 


0.033 ± 0.033 


04D4gg 


22 


723 


+ 





031 


1.204 ±0.170 


0.074 ± 0.028 


22 


769 


+ 





030 


1 


139 


+ 





016 


0.138 ±0.016 


04D4gz 


23 


392 


+ 





038 


0.136 ±0.236 


0.153 ±0.031 


23 


405 


+ 





037 





990 


+ 





035 


0.145 ± 0.028 


04D4hf 


24 


786 


+ 





082 


0.879 ± 0.628 


0.026 ± 0.078 


24 


786 


+ 





102 


1 


070 


+ 





065 


0.014 ±0.059 


04D4ht 


22 


656 


+ 





030 


0.765 ±0.152 


0.437 ± 0.027 


22 


683 


+ 





020 


1 


059 


+ 





024 


0.476 ± 0.016 


04D4hu 


23 


921 


+ 





032 


0.202 ± 0.255 


-0.131 ±0.054 


23 


891 


± 





046 





992 


+ 





027 


-0.065 ± 0.038 


04D4ib 


23 


590 


+ 





030 


1.050 ±0.238 


-0.144 ±0.050 


23 


571 


+ 





044 


1 


088 


+ 





025 


-0.080 ± 0.037 


04D4ic 


24 


111 


+ 





044 


-0.935 ± 0.424 


-0.022 ±0.071 


24 


106 


+ 





054 





904 


+ 





041 


-0.032 ± 0.044 


04D4id 


24 


205 


+ 





046 


0.725 ± 0.452 


-0.185 ±0.058 


24 


192 


+ 





057 


1 


078 


+ 





043 


-0.070 ± 0.045 


04D4ih 


24 


387 


+ 





055 


0.185 ±0.307 


-0.237 ± 0.049 


24 


428 


+ 





075 


1 


041 


+ 





041 


-0.144 ±0.046 


04D4ii 


24 


382 


+ 





050 


1.515 ±0.500 


0.002 ±0.061 


24 


384 


+ 





060 


1 


169 


+ 





042 


0.043 ± 0.046 


04D4im 


23 


814 


+ 





035 


0.808 ± 0.264 


-0.012 ± 0.056 


23 


863 


+ 





047 


1 


212 


+ 





029 


0.076 ± 0.039 


04D4in 


22 


894 


+ 





030 


1.660 ±0.209 


-0.079 ± 0.030 


22 


900 


+ 





033 


1 


165 


+ 





015 


-0.024 ±0.021 


04D4jr 


22 


642 


+ 





027 


1.521 ±0.187 


-0.057 ± 0.029 


22 


631 


+ 





029 


1 


159 


+ 





016 


-0.027 ±0.019 


04D4ju 


23 


756 


+ 





027 


0.537 ± 0.225 


0.146 ±0.035 


23 


774 


+ 





032 


1 


038 


+ 





030 


0.187 ±0.029 


04D4jw 


24 


769 


+ 





081 


-0.676 ± 0.376 


-0.233 ±0.061 


24 


875 


+ 





179 





863 


+ 





050 


-0.187 ±0.065 


04D4jy 


24 


717 


+ 





074 


0.030 ± 0.473 


-0.120 ±0.069 


24 


771 


+ 





129 


1 


047 


+ 





058 


-0.068 ± 0.058 


05Dlaz 


24 


226 


+ 





044 


1.730 ±0.417 


-0.003 ± 0.058 


24 


248 


+ 





056 


1 


187 


+ 





038 


0.001 ± 0.044 


05Dlcb 


23 


699 


+ 





027 


-0.203 ± 0.214 


-0.015 ± 0.041 


23 


708 


+ 





041 





958 


+ 





026 


-0.018 ±0.027 


05Dlcc 


23 


505 


+ 





029 


-0.023 ±0.167 


-0.021 ±0.032 


23 


474 


+ 





040 





968 


+ 





023 


-0.017 ±0.026 


05Dlck 


24 


066 


+ 





033 


0.109 ±0.258 


0.100 ±0.039 


24 


060 


+ 





044 





997 


+ 





031 


0.099 ± 0.034 
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Table 11. continued. 



Name 


m*(SALT2) 


Xi(SALT2) 


C(SALT2) 


m*(SIFTO) 




s(SIFTO) 


C(SIFTO) 


05Dlcl 


24 


350 


+ 





048 


2.358 + 0.436 


-0.040 ± 0.056 


24 


323 


+ 





054 


1 


236 


± 





037 


-0.009 ± 0.043 


05Dldn 


23 


300 


+ 





028 


1.336 ±0.257 


-0.016 ± 0.036 


23 


319 


+ 





038 


1 


117 


± 





021 


0.020 ± 0.037 


05Dldx 


23 


293 


+ 





025 


0.716 ±0.153 


-0.066 ± 0.030 


23 


299 


+ 





037 


1 


081 


± 





017 


-0.015 ±0.027 


05Dlee 


23 


524 


+ 





030 


-0.560 ± 0.202 


-0.018 ±0.036 


23 


574 


+ 





040 





966 


± 





027 


0.033 ± 0.034 


05Dlej 


22 


775 


+ 





032 


-1.264 ±0.136 


0.265 ± 0.028 


22 


875 


+ 





031 





894 


± 





018 


0.339 ± 0.020 


05Dlem 


24 


265 


+ 





044 


0.342 ± 0.303 


-0.141 ±0.055 


24 


271 


+ 





059 


1 


003 


± 





030 


-0.1 14 ±0.043 


05Dleo 


24 


310 


+ 





037 


-1.343 ±0.313 


-0.076 ± 0.064 


24 


296 


+ 





050 





857 


± 





031 


-0.047 ± 0.040 


05Dler 


24 


606 


+ 





052 


0.516 ±0.465 


0.062 ± 0.069 


24 


599 


+ 





065 


1 


030 


± 





044 


0.026 ± 0.047 


05Dlhk 


21 


160 


+ 





036 


1.443 ±0.186 


-0.038 ± 0.028 


21 


200 


+ 





023 


1 


185 


± 





015 


-0.005 ± 0.014 


05Dlhm 


24 


087 


+ 





037 


1.220 ±0.339 


0.129 ± 0.041 


24 


101 


+ 





046 


1 


121 


± 





030 


0.153 ±0.038 


05Dlhn 


21 


155 


+ 





042 


0.654 ±0.198 


0.306 ± 0.032 


21 


171 


+ 





030 


1 


051 


± 





020 


0.344 ± 0.024 


05Dlif 


24 


031 


+ 





036 


0.191 ±0.258 


-0.074 ± 0.050 


24 


061 


+ 





049 


1 


019 


± 





030 


-0.030 ± 0.037 


05Dlix 


22 


874 


+ 





027 


0.626 ±0.162 


-0.072 ± 0.027 


22 


876 


± 





034 


1 


048 


± 





020 


-0.031 ±0.020 


05Dliy 


21 


685 


+ 





041 


-1.433 ±0.169 


-0.037 ± 0.033 


21 


720 


± 





030 





843 


± 





013 


0.011 ±0.023 


05Dliz 


24 


390 


+ 





071 


1.346 ±0.727 


-0.124 ±0.076 


24 


363 


+ 





082 


1 


052 


± 





064 


-0.089 ±0.051 


05Dlju 


24 


569 


+ 





048 


-0.241 ± 0.480 


0.247 ± 0.074 


24 


545 


+ 





060 





950 


± 





047 


0.238 ± 0.058 


05Dlkl 


24 


148 


+ 





040 


0.744 ±0.331 


0.134 ±0.042 


24 


147 


+ 





045 


1 


070 


± 





037 


0.162 ±0.041 


05D2ab 


21 


985 


+ 





026 


-0.131 ±0.104 


-0.047 ± 0.023 


22 


004 


± 





026 





972 


± 





023 


-0.012 ±0.017 


05D2ac 


22 


657 


+ 





026 


1.083 ±0.146 


-0.049 ± 0.027 


22 


685 


+ 





031 


1 


157 


± 





015 


-0.009 ±0.017 


05D2ah 


20 


749 


+ 





028 


-0.001 ± 0.080 


-0.010 ± 0.025 


20 


770 


+ 





019 





988 


+ 





026 


0.013 ±0.015 


05D2ay 


24 


661 


+ 





070 


-0.077 ± 0.726 


0.009 ± 0.092 


24 


640 


+ 





082 





976 


± 





064 


-0.004 ± 0.055 


05D2bt 


23 


505 


+ 





034 


0.606 ± 0.295 


-0.172 ± 0.047 


23 


509 


+ 





045 


1 


069 


± 





026 


-0.094 ± 0.033 


05D2bv 


22 


717 


+ 





026 


0.030 ±0.141 


-0.125 ±0.025 


22 


705 


+ 





034 





975 


± 





024 


-0.107 ±0.018 


05D2bw 


24 


348 


+ 





069 


-0.438 ± 0.505 


-0.173 ±0.064 


24 


413 


+ 





087 





994 


± 





063 


-0.091 ±0.052 


05D2by 


24 


556 


+ 





052 


1.470 ±0.548 


-0.031 ±0.061 


24 


543 


+ 





065 


1 


102 


± 





042 


-0.020 ± 0.046 


05D2cb 


23 


392 


+ 





034 


0.850 ± 0.209 


0.157 ± 0.030 


23 


409 


+ 





033 


1 


115 


± 





024 


0.196 ±0.022 


05D2ci 


23 


566 


+ 





050 


-1.048 ±0.298 


0.039 ± 0.066 


23 


639 


+ 





053 





915 


± 





028 


0.026 ± 0.042 


05D2ck 


24 


474 


+ 





042 


-2.312 ±0.471 


-0.041 ±0.071 


24 


448 


+ 





051 





719 


± 





041 


-0.031 ±0.043 


05D2ct 


24 


379 


+ 





058 


0.434 ±0.5 11 


0.070 ± 0.067 


24 


364 


+ 





060 


1 


042 


± 





047 


0.088 ± 0.048 


05D2dt 


23 


650 


+ 





034 


0.265 ± 0.296 


0.022 ± 0.037 


23 


639 


+ 





046 


1 


007 


± 





049 


0.034 ± 0.040 


05D2dw 


22 


469 


+ 





031 


1 .094 ± 0.203 


-0.018 ±0.027 


22 


489 


± 





029 


1 


133 


± 





022 


0.023 ±0.015 


05D2dy 


22 


921 


+ 





029 


1.252 ±0.293 


-0.136 ±0.031 


22 


891 


+ 





034 


1 


068 


± 





019 


-0.099 ± 0.024 


05D2eb 


23 


007 


+ 





038 


1.151 ±0.348 


-0.062 ± 0.034 


22 


993 


+ 





041 


1 


104 


± 





035 


-0.043 ± 0.029 


05D2ec 


23 


675 


+ 





032 


-0.061 ±0.226 


-0.113 ±0.051 


23 


648 


+ 





047 





992 


± 





026 


-0.053 ± 0.038 


05D2ei 


22 


143 


+ 





030 


0.640 ±0.169 


0.090 ± 0.026 


22 


228 


+ 





020 


1 


098 


± 





023 


0.162 ±0.017 


05D2fq 


24 


000 


+ 





034 


0.914 ±0.334 


-0.096 ± 0.054 


23 


965 


± 





047 


1 


076 


± 





034 


-0.011 ±0.039 


05D2hc 


22 


683 


+ 





027 


-0.520 ±0.179 


0.031 ± 0.025 


22 


692 


± 





027 





918 


± 





023 


0.047 ± 0.017 


05D2he 


23 


933 


+ 





032 


0.383 ±0.438 


0.061 ± 0.043 


23 


954 


+ 





043 


1 


045 


± 





040 


0.066 ± 0.035 


05D2ie 


22 


223 


+ 





030 


-0.136 ±0.164 


-0.082 ± 0.027 


22 


267 


+ 





031 


1 


000 


± 





020 


-0.045 ± 0.020 


05D21e 


23 


955 


+ 





035 


0.578 ± 0.326 


-0.047 ± 0.049 


23 


941 


+ 





049 


1 


048 


± 





035 


-0.012 ±0.037 


05D2mp 


22 


415 


+ 





040 


1.395 ±0.277 


0.025 ± 0.028 


22 


405 


+ 





030 


1 


122 


± 





020 


0.060 ±0.018 


05D2my 


24 


649 


+ 





063 


1.376 ±0.524 


-0.053 ± 0.058 


24 


667 


+ 





079 


1 


114 


± 





050 


-0.024 ± 0.050 


05D2nn 


24 


340 


+ 





077 


-1.346 ±0.540 


-0.242 ± 0.068 


24 


414 


+ 





081 





904 


± 





060 


-0.106 ±0.051 


05D2nt 


24 


091 


+ 





034 


1.592 ±0.362 


-0.039 ±0.051 


24 


080 


+ 





048 


1 


148 


± 





029 


-0.002 ± 0.036 


05D2ob 


24 


780 


+ 





058 


1 .004 ± 0.678 


0.019 ± 0.066 


24 


822 


+ 





074 


1 


101 


± 





062 


0.021 ±0.050 


05D3ax 


23 


605 


+ 





052 


0.624 ± 0.477 


-0.078 ± 0.063 


23 


604 


+ 





052 


1 


060 


± 





051 


-0.021 ± 0.040 


05D3cf 


22 


949 


+ 





028 


-0.154 ±0.181 


0.011 ± 0.028 


22 


967 


+ 





030 





963 


± 





031 


0.042 ± 0.021 


05D3ci 


23 


612 


+ 





064 


2.120 ±0.559 


0.146 ± 0.045 


23 


515 


+ 





049 


1 


116 


± 





034 


0.147 ±0.034 


05D3cq 


24 


179 


+ 





045 


-0.818 ±0.379 


-0.225 ± 0.057 


24 


225 


+ 





061 


1 


008 


± 





050 


-0.137 ±0.046 


05D3cx 


23 


899 


+ 





036 


0.724 ± 0.458 


-0.153 ±0.049 


23 


896 


+ 





051 


1 


029 


± 





045 


-0.112 ±0.044 


05D3dd 


22 


929 


+ 





026 


-0.135 ±0.129 


-0.048 ± 0.025 


22 


938 


+ 





036 





983 


± 





031 


-0.019 ±0.021 


05D3dh 


24 


195 


+ 





039 


0.252 ± 0.287 


0.100 ±0.062 


24 


181 


+ 





055 


1 


000 


± 





031 


0.038 ± 0.046 


05D3gp 


23 


508 


+ 





051 


-0.536 ±0.421 


-0.094 ± 0.049 


23 


523 


+ 





056 





961 


± 





064 


-0.046 ± 0.050 


05D3gv 


23 


997 


+ 





037 


-0.746 ± 0.308 


-0.081 ±0.061 


23 


976 


± 





051 





908 


± 





031 


-0.049 ± 0.039 


05D3na 


24 


384 


+ 





057 


—0.023 ± 0.425 


A AO O i A A£ C 

0.088 ± 0.065 


24 


358 


+ 





067 


A 


948 


± 





039 


AACO i AA/1C 

0.058 ± 0.045 


05D3hh 


24 


295 


+ 





054 


0.678 ± 0.438 


-0.014 ± 0.062 


24 


313 


+ 





058 


1 


064 


± 





042 


0.014 ± 0.043 


05D3hq 


22 


012 


+ 





030 


-0.152 ±0.117 


-0.033 ± 0.025 


22 


080 


+ 





033 


1 


028 


± 





022 


0.005 ± 0.020 


05D3hs 


23 


503 


+ 





034 


0.560 ± 0.229 


-0.192 ±0.047 


23 


477 


+ 





048 


1 


032 


± 





026 


-0.134 ±0.033 


05D3ht 


24 


371 


+ 





058 


0.901 ± 0.404 


-0.149 ± 0.056 


24 


432 


+ 





069 


1 


147 


± 





039 


-0.073 ± 0.046 


05D3jb 


23 


918 


+ 





032 


0.968 ± 0.261 


-0.072 ± 0.050 


23 


932 


+ 





045 


1 


100 


± 





022 


-0.010 ±0.035 


05D3jh 


23 


737 


+ 





031 


-0.527 ±0.197 


-0.145 ±0.049 


23 


737 


+ 





048 





913 


± 





021 


-0.094 ± 0.034 
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Table 11. continued. 



Name 


m*(SALT2) 


Xi(SALT2) 


C(SALT2) 


m*(SIFTO) 




s(SIFTO) 


C(SIFTO) 


05D3jk 


23 


699 


+ 





028 


0.923 ± 0.202 


-0.130 ±0.046 


23 


704 


+ 





044 


1 


060 


± 





024 


-0.064 ± 0.033 


05D3jq 


23 


304 


+ 





026 


1.456 ±0.280 


-0.002 ± 0.030 


23 


321 


+ 





036 


1 


162 


± 





020 


0.041 ± 0.022 


05D3jr 


22 


654 


+ 





027 


-0.829 ±0.138 


0.063 ± 0.024 


22 


660 


+ 





028 





897 


± 





019 


0.097 ±0.016 


05D3km 


24 


706 


+ 





058 


0.518 ±0.337 


-0.163 ±0.049 


24 


782 


+ 





080 





990 


± 





038 


-0.144 ±0.047 


05D3kp 


24 


117 


+ 





036 


1.335 ±0.275 


-0.128 ± 0.049 


24 


128 


± 





052 


1 


122 


± 





020 


-0.083 ± 0.041 


05D3kt 


23 


961 


+ 





030 


-0.347 ± 0.204 


0.039 ± 0.048 


23 


950 


+ 





043 


1 


001 


± 





028 


0.102 ±0.032 


05D3kx 


20 


866 


+ 





029 


0.759 ±0.108 


-0.040 ± 0.025 


20 


855 


+ 





017 


1 


055 


± 





015 


-0.031 ±0.014 


05D31a 


24 


470 


+ 





051 


-0.163 ±0.282 


-0.124 ±0.051 


24 


471 


+ 





068 





952 


± 





035 


-0.082 ± 0.043 


05D31b 


23 


887 


+ 





028 


0.437 ±0.197 


0.017 ± 0.048 


23 


882 


+ 





042 


1 


035 


± 





021 


0.027 ± 0.032 


05D31c 


22 


973 


+ 





024 


-0.597 ±0.1 19 


-0.052 ± 0.025 


22 


979 


+ 





028 





899 


± 





016 


-0.026 ±0.016 


05D31r 


23 


850 


+ 





043 


0.283 ±0.412 


0.088 ± 0.041 


23 


839 


+ 





047 


1 


004 


± 





039 


0.083 ± 0.035 


05D3mh 


24 


073 


+ 





039 


1.132 ±0.490 


0.064 ± 0.054 


24 


109 


+ 





053 


1 


043 


± 





041 


0.023 ± 0.035 


05D3mn 


24 


032 


+ 





029 


0.234 ± 0.298 


-0.044 ± 0.050 


24 


028 


± 





045 





955 


± 





035 


-0.033 ± 0.037 


05D3mq 


21 


498 


+ 





037 


-0.760 ± 0.288 


-0.006 ±0.031 


21 


541 


± 





023 





915 


± 





021 


0.045 ± 0.017 


05D3mx 


23 


028 


+ 





030 


-1.424 ±0.150 


-0.092 ± 0.035 


23 


046 


+ 





033 





822 


± 





015 


-0.059 ± 0.025 


05D3ne 


20 


246 


+ 





049 


-1.660 ±0.294 


-0.181 ±0.039 


20 


244 


+ 





030 





801 


± 





017 


-0.146 ±0.022 


05D4af 


23 


095 


+ 





028 


0.204 ± 0.245 


-0.044 ± 0.030 


23 


111 


+ 





037 


1 


005 


± 





039 


-0.016 ±0.032 


05D4ag 


23 


877 


+ 





046 


0.505 ± 0.475 


0.075 ± 0.059 


23 


888 


+ 





053 


1 


051 


± 





069 


0.045 ± 0.047 


05D4av 


23 


547 


+ 





028 


0.920 ±0.199 


0.154 ±0.031 


23 


559 


+ 





035 


1 


096 


± 





019 


0.183 ±0.022 


05D4ay 


24 


244 


+ 





054 


-2.069 ±0.319 


0.388 ±0.051 


24 


368 


+ 





057 





863 


± 





040 


0.424 ± 0.038 


05D4be 


22 


909 


+ 





026 


1.070 ±0.168 


-0.167 ±0.028 


22 


911 


+ 





035 


1 


085 


± 





017 


-0.094 ± 0.023 


05D4bf 


23 


621 


+ 





031 


0.295 ± 0.260 


-0.016 ± 0.040 


23 


616 


+ 





044 


1 


025 


± 





032 


0.020 ± 0.041 


05D4bi 


24 


044 


+ 





041 


1.160 ±0.402 


-0.135 ±0.062 


24 


035 


+ 





050 


1 


077 


± 





041 


-0.092 ± 0.044 


05D4bj 


24 


096 


+ 





033 


0.187 ±0.283 


0.034 ± 0.055 


24 


084 


+ 





047 


1 


007 


± 





036 


0.054 ± 0.040 


05D4bm 


22 


221 


+ 





026 


0.284 ± 0.087 


-0.070 ± 0.023 


22 


214 


+ 





028 


1 


019 


± 





013 


-0.042 ±0.018 


05D4cn 


24 


076 


+ 





034 


1.178 ±0.269 


-0.028 ± 0.056 


24 


102 


+ 





047 


1 


100 


± 





025 


0.029 ± 0.037 


05D4cq 


23 


740 


+ 





031 


0.837 ± 0.289 


-0.1 18 ±0.052 


23 


695 


+ 





045 


1 


069 


± 





027 


-0.036 ± 0.038 


05D4cs 


23 


966 


+ 





033 


1.406 ±0.244 


-0.120 ± 0.048 


23 


942 


± 





048 


1 


123 


± 





021 


-0.039 ± 0.041 


05D4cw 


22 


146 


+ 





026 


-0.753 ±0.109 


-0.140 ±0.025 


22 


133 


+ 





029 





911 


± 





015 


-0.134 ±0.018 


05D4dt 


22 


807 


+ 





026 


-0.899 ±0.1 15 


-0.049 ± 0.025 


22 


790 


+ 





029 





885 


± 





014 


-0.034 ±0.017 


05D4dw 


24 


420 


+ 





043 


0.574 ±0.373 


-0.010 ± 0.059 


24 


423 


+ 





058 


1 


034 


± 





036 


-0.019 ± 0.044 


05D4dy 


24 


596 


+ 





044 


-0.267 ± 0.406 


-0.129 ±0.063 


24 


606 


+ 





057 





945 


± 





043 


-0.095 ± 0.046 


05D4ef 


23 


819 


+ 





036 


-1.419 ±0.233 


-0.087 ± 0.040 


23 


825 


+ 





045 





832 


± 





020 


-0.063 ± 0.035 


05D4ej 


23 


717 


+ 





030 


0.331 ±0.195 


-0.036 ± 0.036 


23 


759 


+ 





040 


1 


033 


± 





024 


0.044 ± 0.032 


05D4ek 


23 


285 


+ 





030 


0.312 ±0.194 


0.048 ±0.031 


23 


309 


+ 





038 


1 


086 


± 





022 


0.084 ± 0.026 


05D4ev 


24 


227 


+ 





043 


-1.060 ±0.321 


-0.063 ± 0.062 


24 


263 


± 





054 





913 


± 





031 


-0.003 ± 0.040 


05D4fe 


24 


907 


+ 





078 


3.119 ±0.811 


-0.067 ± 0.063 


24 


863 


+ 





085 


1 


339 


± 





057 


0.021 ±0.050 


05D4ff 


22 


589 


+ 





028 


-0.635 ±0.128 


-0.008 ± 0.026 


22 


626 


+ 





032 





938 


± 





018 


0.029 ±0.019 


05D4fg 


24 


177 


+ 





040 


0.374 ± 0.279 


-0.155 ±0.054 


24 


181 


+ 





056 


1 


015 


± 





029 


-0.102 ±0.042 


05D4fo 


22 


440 


+ 





026 


-0.674 ± 0.099 


-0.063 ± 0.023 


22 


475 


+ 





028 





926 


± 





013 


-0.015 ±0.017 


05D4gw 


24 


455 


+ 





053 


0.188 ±0.456 


0.023 ± 0.065 


24 


476 


+ 





062 


1 


034 


± 





047 


0.034 ± 0.048 


05D4hn 


24 


184 


+ 





055 


1.576 ±0.907 


0.068 ± 0.072 


24 


175 


+ 





057 


1 


079 


± 





049 


0.008 ± 0.045 


06Dlab 


21 


437 


+ 





065 


-2.066 ± 0.290 


0.148 ±0.050 


21 


606 


+ 





030 





812 


± 





016 


0.246 ± 0.025 


06D2bk 


23 


285 


+ 





031 


0.748 ± 0.335 


0.010 ±0.034 


23 


253 


+ 





036 


1 


036 


± 





042 


0.033 ± 0.038 


06D2cb 


24 


621 


+ 





120 


0.563 ± 1.002 


-0.127 ±0.074 


24 


649 


+ 





159 


1 


077 


± 





118 


-0.058 ± 0.075 


06D2cc 


23 


454 


+ 





044 


-0.463 ± 0.372 


0.059 ± 0.050 


23 


471 


± 





048 





941 


± 





040 


0.087 ± 0.042 


06D2ga 


24 


285 


+ 





080 


1.676 ± 1.010 


-0.031 ±0.076 


24 


263 


+ 





071 


1 


176 


± 





073 


0.013 ±0.049 


06D2gb 


22 


997 


+ 





041 


-1.566 ±0.258 


-0.000 ± 0.043 


23 


002 


+ 





045 





831 


± 





026 


0.021 ±0.030 


06D3bz 


23 


962 


+ 





042 


-0.345 ± 0.363 


-0.106 ±0.064 


23 


929 


+ 





055 





888 


± 





042 


-0.051 ±0.043 


06D3cb 


24 


392 


+ 





109 


0.696 ± 0.542 


-0.121 ±0.067 
























06D3cc 


24 


055 


+ 





066 


0.856 ± 0.557 


-0.023 ± 0.080 


24 


050 


+ 





066 


1 


072 


± 





054 


0.011 ±0.045 


06D3df 


22 


675 


+ 





025 


1.128 ±0.217 


-0.024 ± 0.024 


22 


679 


+ 





028 


1 


127 


± 





027 


0.028 ±0.015 


06D3do 


23 


900 


+ 





035 


0.713 ±0.553 


-0.086 ±0.051 


23 


865 


+ 





052 


1 


117 


± 





067 


-0.033 ± 0.052 


06D3dt 


22 


161 


+ 





027 


0.146 ±0.113 


0.092 ± 0.023 


22 


168 


+ 





020 





945 


± 





031 


0.112 ±0.013 


06D3ed 


22 


614 


+ 





025 


-0.236 ± 0.085 


-0.064 ± 0.023 


22 


599 


+ 





028 





957 


± 





015 


-0.043 ±0.016 


06D3el 


22 


895 


+ 





027 


0.880 ±0.164 


-0.145 ±0.029 


22 


922 


+ 





032 


1 


069 


+ 





014 


-0.055 ± 0.020 


06D3em 


24 


366 


+ 





035 


0.071 ±0.384 


0.150 ±0.058 


24 


361 


+ 





050 


1 


009 


± 





039 


0.130 ±0.039 


06D3en 


24 


674 


+ 





087 


-1.143 ±0.473 


-0.195 ±0.061 


24 


753 


+ 





157 





839 


± 





048 


-0.154 ±0.061 


06D3et 


23 


455 


+ 





031 


-1.329 ±0.188 


-0.103 ±0.037 


23 


555 


± 





038 





870 


± 





014 


-0.009 ± 0.030 


06D3fp 


21 


735 


+ 





027 


0.049 ±0.081 


0.061 ±0.023 


21 


751 


+ 





020 





992 


± 





018 


0.115 ±0.014 


06D3gh 


23 


922 


+ 





044 


0.411 ±0.380 


-0.075 ± 0.063 


23 


897 


+ 





053 


1 


019 


± 





044 


-0.006 ± 0.042 
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Table 11. continued. 



Name 



m^(SALT2) 



Xi(SALT2) 



C(SALT2) 



w*(SIFTO) s(SIFTO) 



C(SIFTO) 



06D3gn 

06D4ba 

06D4bo 

06D4bw 

06D4ce 

06D4cl 

06D4co 

06D4cq 

06D4dh 



21.843: 
23.733: 
23.227: 
23.898: 
24.183: 
24.494 : 
22.513: 
22.561 : 
21.458: 



0.032 
0.049 
0.028 
0.036 
0.065 
0.062 
0.024 
0.024 
0.027 



-0.531+0.136 
0.722 + 0.414 
0.959 + 0.226 
0.759 + 0.269 
1.631 ±0.636 
1.182 ±0.419 
-0.306 ±0.100 
0.481 ±0.111 
0.592 ± 0.086 



0.102 ± 
-0.134 d 
-0.047 d 
-0.005 d 
-0.109 d 
-0.099 d 
-0.053 d 
-0.041 d 
-0.151 d 



0.028 
: 0.067 
: 0.030 
: 0.058 
: 0.061 
: 0.054 
: 0.023 
: 0.022 
: 0.022 



21.931 
23.760 : 
23.221 : 
23.883: 
24.195 : 
24.598 : 
22.513: 
22.550 : 
21.437: 



0.020 
0.058 
0.037 
0.050 
0.065 
0.089 
0.028 
0.028 
0.027 



0.963 ±0.015 
1.053 ±0.045 
1.097 ±0.027 
1.057 ±0.029 
1.165 ±0.047 
1.132 ±0.038 
0.959 ±0.016 
1.033 ±0.015 
1.049 ±0.017 



0.185 ±0.018 
-0.063 ± 0.045 
-0.038 ± 0.029 
-0.008 ± 0.042 
-0.046 ± 0.047 
-0.071 ±0.050 
-0.041 ±0.015 
-0.003 ±0.017 
-0.129 ±0.017 



Notes. An electronic version of this table including covariances is available at the Centre de Donnees astronomiques de Strasbourg 
(CDS). The SALT2 and SiFTO fit parameters can also be downloaded at the University of Toronto's Research Repository 
https ://tspace. library .utoronto. ca/handle/1807/245 12. 
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Table 12. The SALT2 low-z training sample of Type la supemovae light curves 
and spectra 



SN name 


z a 


Bands* 


Number of spectra* 


snl981b 


0.0060 


UBVR (B83) 





snl981d 


0.0059 


UBV (H91) 





snl989b 


0.0024 


UBVRI (W94) 


3 (AC) 


snl990af 


0.0503 


BV (H96) 





snl990n 


0.0034 


UBVRI (L91,L98) 


12 (AC,IUE,CFA) 


snl991t 


0.0060 


UBVRI (F92,L98,A04,K04b) 


32 (AQCFA) 


snl992a 


0.0063 


UBVRI (S92,A04) 


19 (CFAJUE) 


snl992al 


0.0146 


BVRI (H96) 





snl992bc 


0.0202 


BVRI (H96) 





snl992bo 


0.0185 


BVRI (H96) 





snl993o 


0.0510 


BVI (H96) 





snl994ae 


0.0043 


UBVRI (R99,A04) 


18 (CPA) 


snl994d 


0.0015 


UBVRI (R95,P96,M96,A04) 


24 (AC) 


sn 1994s 


0.0152 


BVRI (R99) 


3 (AC) 


snl995ac 


0.0499 


BVRI (R99) 


5 (CFA) 


snl995al 


0.0051 


BVRI (R99) 


2 (CFA) 


snl995bd 


0.0146 


BVRI (R99) 


9 (CFA) 


snl995d 


0.0066 


BVRI (R99,A04) 


12 (CFA) 


snl996bl 


0.0360 


BVRI (R99) 





snl996bo 


0.0173 


BVRI (R99,A04) 





snl996x 


0.0069 


UBVRI (R99) 


9 (AC) 


snl997bp 


0.0083 


UBVRI (J05) 





snl997bq 


0.0094 


UBVRI (J05) 





snl997do 


0.0101 


UBVRI (J05) 


12 (M08) 


snl997e 


0.0135 


UBVRI (J05) 





snl998ab 


0.0271 


UBVRI (J05) 


1 1 (M08) 


snl998aq 


0.0037 


UBVRI (R05) 


29 (M08) 


snl998bu 


0.0030 


UBVRI (S99) 


35 (M08) 


snl998dh 


0.0089 


UBVRI (J05) 


10 (M08) 


snl998ef 


0.0177 


UBVRI (J05) 





snl998es 


0.0106 


UBVRI (J05) 


26 (M08) 


snl999aa 


0.0144 


UBVRI (J05) 


28 (M08) 


snl999ac 


0.0095 


UBVRI (J05) 


19 (M08) 


snl999aw 


0.0380 


BVRI (St02) 





snl999bp 


0.0770 


BVRI (K08) 





snl999cc 


0.0313 


UBVRI (J05,K06) 





snl999cp 


0.0095 


BVRI (J05,K00) 





snl999dk 


0.0150 


UBVRI (K01,A04) 





snl999dq 


0.0143 


UBVRI (J05) 


23 (M08) 


snl999ee 


0.0114 


UBVRI (S02) 


12 (H02) 


snl999ek 


0.0175 


UBVRI (J05,K04b) 





snl999gp 


0.0267 


UBVRI (J05,K01) 


10 (M08) 


sn2000cn 


0.0235 


UBVRI (J05) 





sn2000dk 


0.0174 


UBVRI (J05) 





sn2000e 


0.0047 


UBVRI (V03) 


5 (V03) 


sn2000fa 


0.0213 


UBVRI (J05) 


14 (M08) 


sn2001ba 


0.0294 


BVI (K04a) 





sn2001bt 


0.0146 


BVRI (K04b) 





sn2001cz 


0.0154 


UBVRI (K04b) 





sn2001el 


0.0039 


UBVRI (K03) 





sn2001ep 


0.0130 


UBVRI (H09) 





sn2001fh 


0.0130 


UBVRI (H09) 





sn2001v 


0.0150 


UBVRI (Vi03,H09) 





sn2002bf 


0.0244 


UBVRI (H09) 





sn2002bo 


0.0042 


UBVRI (Z03,B04,K04b,H09) 


15 (B04) 


sn2002de 


0.0281 


UBVRI (H09) 





sn2002dj 


0.0094 


UBVRI (H09) 


16 


sn2002dp 


0.0116 


UBVRI (H09) 





sn2002er 


0.0086 


UBVRI (P04) 





sn2002ha 


0.0140 


UBVRI (H09) 
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Table 12. continued. 



SN name 


z" 


Bands* 


Number of spectra* 


sn2002he 


0.0246 


UBVRI (H09) 





sn2002hu 


0.0300 


UBVRI (H09) 





sn2003du 


0.0064 


UBVRI (A05) 


7 (A05,H09) 


sn2003fa 


0.0060 


UBVRI (H09) 





sn2003kf 


0.0074 


UBVRI (H09) 





sn2003u 


0.0263 


UBVRI (H09) 





sn2003w 


0.0201 


UBVRI (H09) 





sn2004as 


0.0310 


UBVRI (H09) 





sn2004eo 


0.0157 


UBVRI (P07) 





sn2004fu 


0.0092 


UBVRI (T06,H09) (H09) 





sn2005cf 


0.0064 


UBVRIubv (H09) 





sn2005el 


0.0149 


UBVri (H09) 





sn2005eq 


0.0290 


UBVri (H09) 





sn2005eu 


0.0349 


BVri (H09) 





sn2005hc 


0.0459 


UBVri (H09) 





sn2005iq 


0.0340 


UBVri (H09) 





sn2005kc 


0.0151 


UBVri (H09) 





sn2005ki 


0.0192 


UBVri (H09) 





sn20051z 


0.0410 


UBVri (H09) 





sn2005ms 


0.0252 


UBVri (H09) 





sn2006ac 


0.0231 


UBVri (H09) 





sn2006ax 


0.0167 


UBVri (H09) 





sn2006az 


0.0309 


UBVri (H09) 





sn2006bt 


0.0322 


UBVri (H09) 





sn2006cc 


0.0325 


UBVri (H09) 





sn2006d 


0.0085 


UBVri (H09) 





sn2006gr 


0.0346 


UBVri (H09) 





sn2006kf 


0.0213 


BVri (H09) 





sn20061e 


0.0174 


UBVri (H09) 





sn20061f 


0.0132 


UBVri (H09) 





sn2006mp 


0.0230 


BVri (H09) 





sn2006oa 


0.0600 


UBVri (H09) 





sn2006ob 


0.0592 


UBVri (H09) 





sn2006s 


0.0321 


UBVri (H09) 





sn2006sr 


0.0241 


UBVri (H09) 





sn2007af 


0.0055 


BVubvn (H09) 





sn2007au 


0.0206 


BVri (H09) 





sn2007bc 


0.0208 


BVri (H09) 





sn2007bd 


0.0310 


BVri (H09) 





sn2007ci 


0.0181 


BVri (H09) 





sn/UU/co 


U.UZ /U 


UrJVUDVll (HUV) 


n 

u 


sn2007cq 


0.0259 


BVubvri (H09) 





sn2007f 


0.0236 


BVri (H09) 





sn2007qe 


0.0240 


BVri (H09) 





sn2008bf 


0.0240 


UBVri (H09) 





sn2006qo 


0.0305 


BVri (H09) 





sn2007ca 


0.0141 


BVri (H09) 






" Heliocentric redshift. 
b Ref erences. B83: 
(1199 I) . L91: 
R95: 
R99: 



Richmond 



Leibund 
id et alJdl 
al.l(ll999h" 



Buta & Turne"rl dl983h . 

199 lh . F92 



Leibundgut et all 
" 19951) . H96: 



Br83: iBranch et al 



Filiorjenko et al 



.. (119831) . P 87: iPhillips et al l dl987l) . H91: iHamuv et alJ 
"(Il992h S92: ISuntzeffl d!992h W9 4: IWeUs et alJ 



Hamuv et al.Nl996l) . M96:lMeikle et all d 1 9961) . P9 6: iPatat et al.l d 1 9961). L9 8~ lLira et al 



1994), 



1998), 



. JRiess et 
d2002l), S02:"lStritzinger et al 



S99: lSuntzeffetal 
d20O2h. ~ 



19991). K00:lKrisciunas et alJd2000h. K01:lKrisciunas et al.ld200ll). H02:lHamuv et al 



1-11199' 

., St02: IStrolger et all d2002h. K03 : iKrisciunas et alJ d2003h. V03: IValentini et all d2003h 
Vi03: IVinko et alJ d2003l). W03: IWang et all d2003ah. Z03: IZapataetal] d2003l). A04: lAltavilla et ap d2004l). B04: iBenetti et al.l 
d200- ' 



G04: iGaravini etaD d2004l) , K04a: IKrisciunas et alJ (2004al). K0 4 b: IKrisciunas et alJ d2004bl P04: IPignata et al.l d20O 
A05: jAnupama et alJ d2005h J05: " 



Jha et al l d2006h. R05: iRTess et alJ d2005h. K06 
K08: iKowalski et al.1 d2008h . M08: iMatheson et al.1 



(200 6J), P07: |P astorello^eLalJ (12007 ) 

AQ lBarbon et al.l d 19991) . IUE: llNEsld2006l) . CFA: www.cfa.harvard.edu/oir/Research/supernova/. 



Krisciunas et"aD d2006l). T06: iTsvetkovl 
d2008l) . H09: iHicken et al.1 d2009al) . 



36 



J. Guy et al., SNLS Collaboration: SNLS-3: SNe la photometric distances and cosmological constraints 



References 

Adelman-McCarthy, J. K., Agiieros, M. A., Allam, S. S., et al. 2007, ApJS, 172, 
634 

Altavilla, G. et al. 2004, Mon. Not. Roy. Astron. Soc, 349, 1344 

Anupama, G. C, Sahu, D. K., & Jose, J. 2005, A&A, 429, 667 

Appenzeller, I., Fricke, K., Fiirtig, W., et al. 1998, The Messenger, 94, 1 

Astier, P., Guy, J., Regnault, N., et al. 2006, A&A, 447, 31 

Balland, C, Baumont, S., Basa, S., et al. 2009, A&A, 507, 85 

Barbon, R., Buondf, V., Cappellaro, E., & Turatto, M. 1999, A&AS, 139, 531 

Baron, E. & Hauschildt, P. H. 1998, ApJ, 495, 370 

Baumont, S., Balland, C, Astier, P., et al. 2008, A&A, 491, 567 

Benetti, S., Meikle, P., Stehle, M., et al. 2004, MNRAS, 348, 261 

Bertin, E. & Arnouts, S. 1996, A&AS, 117, 393 

Bohlin, R. C. & Gilliland, R. L. 2004, AJ, 128, 3053 

Boulade, O., Chariot, X., Abbon, P., et al. 2003, in Instrument Design and 
Performance for Optical/Infrared Ground-based Telescopes. Edited by lye, 
Masanori; Moorwood, Alan F. M. Proceedings of the SPIE, Volume 4841, 
pp. 72-81 (2003)., 72-81 
Branch, D., Lacy, C. H., McCall, M. L., et al. 1983, ApJ, 270, 123 
Bronder, T. J., Hook, I. M., Astier, P., et al. 2008, A&A, 477, 717 
Buta, R. J. & Turner, A. 1983, PASP, 95, 72 
Cardelli, J. A., Clayton, G. C, & Mathis, J. S. 1989, APJ, 345, 245 
Conley, A., Carlberg, R. G, Guy, J., et al. 2007, ApJ, 664, L13 
Conley, A., Goldhaber, G, Wang, L., et al. 2006, ApJ, 644, 1 
Conley, A., Guy, J., Sullivan, M., et al. in prep., to ApJ 
Conley, A., Sullivan, M., Hsiao, E. Y., et al. 2008, ApJ, 681, 482 
Ellis, R. S., Sullivan, M., Nugent, P. E., et al. 2008, ApJ, 674, 51 
Faber, S. M., Phillips, A. C, Kibrick, R. I., et al. 2003, in Society of 
Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 
4841, Society of Photo-Optical Instrumentation Engineers (SPIE ) Conference 
Series, ed. M. Iye & A. F. M. Moorwood, 1657-1669 
Fakhouri, H. K. & SNLS collaboration, in prep., to ApJ 
Filippenko, A. V, Richmond, M. W., Matheson, T, et al. 1992, ApJ, 384, L15 
Freedman, W. L., Bums, C. R., Phillips, M. M., et al. 2009, ApJ, 704, 1036 
Garavini, G, Folatelli, G, Goobar, A., et al. 2004, AJ, 128, 387 
Goldhaber, G, Groom, D. E., Kim, A., et al. 2001, ApJ, 558, 359 
Guy, J., Astier, P., Baumont, S., et al. 2007, A&A, 466, 11 
Guy, J., Astier, P., Nobili, S., Regnault, N, & Pain, R. 2005, A&A, 443, 781 
Hamuy, M., Maza, J., Pinto, P. A., et al. 2002, AJ, 124, 2339 
Hamuy, M., Phillips, M. M., Maza, J., et al. 1991, AJ, 102, 208 
Hamuy, M., Phillips, M. M., Suntzeff, N. B., et al. 1996, AJ, 1 12, 2408 
Harville, D. A. 1977, JASA, 358, 320 

Hauschildt, P. H. & Baron, E. 1999, J. Comp. Applied Math., 109, 41 

Hauschildt, P. H., Baron, E., & Allard, F. 1997, ApJ, 483, 390 

Hicken, M„ Challis, P., Jha, S., et al. 2009a, ArXiv e-prints 

Hicken, M., Wood-Vasey, W. M., Blondin, S., et al. 2009b, ApJ, 700, 1097 

Holtzman, J. A., Marriner, J., Kessler, R., et al. 2008, AJ, 136, 2306 

Hook, I. M., J0rgensen, I., Allington-Smith, J. R., et al. 2004, PASP, 116, 425 

Howell, D. A., Sullivan, M., Nugent, P. E., et al. 2006, Nature, 443, 308 

Howell, D. A., Sullivan, M., Perrett, K., et al. 2005, ApJ, 634, 1190 

Hsiao, E. Y., Conley, A., Howell, D. A., et al. 2007, ApJ, 663, 1187 

INES. 2006, IUE Newly Extracted Spectra, Website, 

http://ines.vilspa.esa.es/ines/ 
Jha, S., Kirshner, R. P., Challis, P., et al. 2006, AJ, 131, 527 
Jha, S., Riess, A. G, & Kirshner, R. P. 2007, ApJ, 659, 122 
Kasen, D., Ropke, F. K., & Woosley, S. E. 2009, Nature, 460, 869 
Kessler, R„ Becker, A. C, Cinabro, D., et al. 2009, ApJS, 185, 32 
Knop, R. A., Aldering, G, Amanullah, R., et al. 2003, ApJ, 598, 102 
Kowalski, M., Rubin, D., Aldering, G, et al. 2008, ArXiv e-prints, 804 
Krisciunas, K., Hastings, N. C, Loomis, K, et al. 2000, ApJ, 539, 658 
Krisciunas, K., Phillips, M. M., Stubbs, C, et al. 2001, AJ, 122, 1616 
Krisciunas, K., Phillips, M. M., Suntzeff, N. B., et al. 2004a, AJ, 127, 1664 
Krisciunas, K., Prieto, J. L., Garnavich, P. M., et al. 2006, AJ, 131, 1639 
Krisciunas, K., Suntzeff, N. B., Candia, P., et al. 2003, AJ, 125, 166 
Krisciunas, K., Suntzeff, N. B., Phillips, M. M., et al. 2004b, AJ, 128, 3034 
Landolt, A. U. 1992, AJ, 104, 340 

Leibundgut, B., Kirshner, R. P., Filippenko, A. V, et al. 1991, ApJ, 371, L23 
Lira, P., Suntzeff, N. B., Phillips, M. M., et al. 1998, AJ, 116, 1006 
Magnier, E. A. & Cuillandre, J.-C. 2004, PASP, 116, 449 

Markovsky, I. & Huffel, S. V. 2007, Signal Processing, 87, 2283 , special Section: 

Total Least Squares and Errors-in- Variables Modeling 
Matheson, T, Kirshner, R. P., Challis, P., et al. 2008, AJ, 135, 1598 
Meikle, W. P. S., Cumming, R. J., Geballe, T. R., et al. 1996, MNRAS, 281, 263 
Monet, D. G, Levine, S. E., Canzian, B., et al. 2003, AJ, 125, 984 
Nugent, P., Kim, A., & Perlmutter, S. 2002, PASP, 1 14, 803 
Oke, J. B., Cohen, J. G, Carr, M., et al. 1995, PASP, 107, 375 
Pastorello, A., Mazzali, P. A., Pignata, G, et al. 2007, MNRAS, 377, 1531 



Patat, F, Benetti, S., Cappellaro, E., et al. 1996, MNRAS, 278, 111 

Perlmutter, S., Aldering, G, Goldhaber, G, et al. 1999, ApJ, 517, 565 

Perrett, K., Balam, D., Sullivan, M., et al. in prep., to ApJ 

Phillips, M. M. 1993, Astrophysical Journal Letters, 413, L105 

Phillips, M. M., Phillips, A. C, Heathcote, S. R., et al. 1987, PASP, 99, 592 

Pignata, G, Patat, F, Benetti, S., et al. 2004, MNRAS, 355, 178 

Regnault, N., Conley, A., Guy, J., et al. 2009, A&A, 506, 999 

Richmond, M. W., Treffers, R. R., Filippenko, A. V, et al. 1995, AJ, 109, 2121 

Riess, A. G, Filippenko, A. V, Challis, P., et al. 1998, AJ, 116, 1009 

Riess, A. G, Kirshner, R. P., Schmidt, B. P., et al. 1999, AJ, 117, 707 

Riess, A. G, Li, W., Stetson, P. B., et al. 2005, ApJ, 627, 579 

Riess, A. G, Strolger, L„ Tonry, J., et al. 2004, ApJ, 607, 665 

Riess, A. G, Strolger, L.-G., Casertano, S., et al. 2007, ApJ, 659, 98 

Stetson, P. B. 1987, PASP, 99, 191 

Stone, R. C, Pier, J. R., & Monet, D. G. 1999, AJ, 118, 2488 

Stritzinger, M., Hamuy, M., Suntzeff, N. B., et al. 2002, AJ, 124, 2100 

Strolger, L.-G, Smith, R. C, Suntzeff, N. B., et al. 2002, AJ, 124, 2905 

Sullivan, M., Conley, A., Howell, D. A., et al. submitted, MNRAS 

Sullivan, M., Guy, J., Conley, A., et al. in prep., to ApJ 

Sullivan, M., Howell, D. A., Perrett, K., et al. 2006, AJ, 131, 960 

Suntzeff, N. 1992, IAU Colloq., ed. R. McCray (Cambridge:Cambridge 

University Press), 145 
Suntzeff, N. B., Phillips, M. M., Covarrubias, R., et al. 1999, AJ, 117, 1175 
Tsvetkov, D. Y. 2006, astro-ph/060605 1 

Valentini, G, Di Carlo, E., Massi, F, et al. 2003, ApJ, 595, 779 
Vinko, J., Biro, I. B., Csak, B., et al. 2003, A&A, 397, 115 
Walker, E. S., Hook, I. M., Sullivan, M., et al. in prep., to MNRAS 
Wang, L., Baade, D., Hofiich, P., et al. 2003a, ApJ, 591, 1110 
Wang, L., Goldhaber, G, Aldering, G, & Perlmutter, S. 2003b, ApJ, 590, 944 
Wells, L. A., Phillips, M. M., Suntzeff, B., et al. 1994, AJ, 108, 2233 
Wood-Vasey, W. M., Miknaitis, G, Stubbs, C. W., et al. 2007, ApJ, 666, 694 
Zapata, A. A., Candia, P., Krisciunas, K, Phillips, M. M., & Suntzeff, N. B. 
2003, American Astronomical Society Meeting Abstracts, 203, 



Appendix A: Technicals details on the training of 
SALT2 

A. 1. Parameters and regularisation 

The SALT2 SED sequence is evaluated for phases between -20 
and +50 days, in the wavelength range from 200 to 920 nm. 
The model is addressed using a basis of third order b-splines 
for which the phase and wavelength axis have been distorted in 
order to allow higher resolution in some areas of the parameter 
space. Those distortions have been tuned to minimise residuals 
from the training data set for a given number of spline nodes. 
For wavelength, the distortion is almost equivalent to a con- 
stant resolution in log space, except in the UV range where we 
lack spectroscopic constraints. For this paper, the number of pa- 
rameters in phase x wavelength is 14 x 100 both for the aver- 
age SED sequence and the first component that is correlated to 
stretch. The magnitude of the colour variation law is modelled 
as C x P(A) where C is the colour parameter and P{A) is a poly- 
nomial with linear extrapolations for wavelength smaller than 
280 nm or larger than 700 nm, with four free coefficients, the 
two others being fixed so that P(Ab) = and P(Ay) = 1. 

The choice of the number of parameters is a trade off be- 
tween the statistical noise of the model, and the minimisation 
of biases introduced by a forced smoothness of the model. For 
instance, SALT2 has many more parameters than SiFTO which 
assumes a pure stretching of the light curves. Also, in some ar- 
eas of the parameter space, we lack data to train the model. For 
instance in the UV range, we have a lot of photometric informa- 
tion from the SNLS dataset at all phases but are missing early- 
and late-time spectroscopy. In this situation, the training fit is 
trying to deconvolve the SN SED from broad-band integrated 
flux measurements, which introduces high frequency noise as 
for any deconvolution procedure. As a consequence, when inte- 
grated in filters other than those used for the training, the model 
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light curves in the UV are quite noisy. This is only a cosmetic 
issue when the statistical uncertainty of the model is propagated 
to the estimates of the light curve parameters (amplitude, shape 
and colour, see 9A.3l l. Still, we can apply priors on the properties 
of the SED sequence to cure part of the deconvolution noise in 
order to improve the statistical accuracy of the model. This how- 
ever systematically introduces biases that have to be controlled 
and minimised. 

A regularisation term for each of the two components (SED 
sequences) is added to the standard x 2 to be minimised of the 
following form: 



w tt' J> [ m iJ+i ~ m uf 



Uj,l 

+ ^ w(i, j, I) [m i+k j+i my - m i+k j my+/] 

iJXl 

where my is the flux value of the SED sequence on a grid of 
phase and wavelength indexed by i and j respectively (similarly 
k and I index offsets in phase and wavelength) . w(i, j, I) is a 
weighting function whose value is inversely proportional to the 
amount of spectroscopic data at the phase and wavelength po- 
sition (/, j); it also decreases for increasing / (wavelength dif- 
ference) so that the regularisation has no effect on the relative 
values of the SED for wavelength distant by 100 nm or more (in 
order to limit the biases on colours). The first term minimises 
the gradient as a function of wavelength when there is no spec- 
troscopic information, and the second is null if the model can be 
factorised as a dyadic product my = a, x bj. Hence, the effect 
of the regularisation is to force the model into a simple inter- 
polation of spectra at different phases when there is no spectro- 
scopic information at intermediate phases, with a minimal im- 
pact on regions that are well constrained by data and a small 
effect on broad-band colours. The strength of this regularisation 
is adjusted via the normalisation of w{i, j, I) so that the value of 
the colour difference (U - B) — (B - V) at maximum luminos- 
ity is altered by less than 0.01. The impact of this regularisation 
on distance moduli obtained with SALT2 is however controlled 
(see {gA2j. 

A.2. Error modeling 

For an optimal light curve fit, we need to evaluate the intrin- 
sic SNe la scatter about the model in order to weight data ac- 
cordingly. A perfect evaluation of the covariance of this remain- 
ing variability as a function of phase and wavelength would be 
equivalent to a complete modeling of SNe la. This is of course 
beyond reach so some approximations have to be made. As 
in G07, we evaluated the intrinsic scatter of a light curve assum- 
ing that the residuals were not correlated as a function of phase, 
and independently broad-band colour uncertainties defined by 
the scatter of light curve amplitudes as a function of the effective 
rest-frame wavelength of the observation pass-bands. 

The intrinsic scatter of light curves as a function of phase, 
effective wavelength and Xi, is evaluated using the covariance 
matrix of the best fit model as described in G07 §6.1. 

The broad-band colour uncertainties (Zf-correction uncer- 
tainties is G07) can be evaluate d by minimisin g the following 
restricted log likelihood (see e. g. lHarvillell 1 9771) : 

F = R T WR - logdet(W) + logdet(A T WA) (A.l) 

where the first term is the usual x 2 (W being the inverse of the 
covariance matrix of the data, and R is the vector of residuals), 



and A is the matrix of derivatives of the model with respect to all 
parameters for each measurement (A T WA is the inverse of the 
covariance matrix of the best fit parameters). The second term is 
expected for a standard log-likelihood, the last one accounts for 
the correlation among residuals introduced by the fitted model. 

In practice, we have to assume a functional form for the 
broad-band colour uncertainties as a function of wavelength 
since only a limited number of parameters can be fit. The fol- 
lowing two functional forms which exhibit different asymptotic 
behaviour have been considered (exponential of polynomial and 
combination of sigmoids): 



/ 4 



k[(A) = exp 

MX) = 



i=0 



a -ai 



+ ai 



1 + exp [(A - A )/AA 

«2 — fli 

f 1 + exp [(-A + A 2 )/AA 2 ] 



(A.2) 



An additional scatter for the observed U -band data is eval- 
uated at the same time, since it suffers from sizable calibration 
systematic uncertainties. 

Instead of minimising the restricted log likelihood of Eq. IA.ll 
on the full training set of SALT2, which comprises light curves 
and spectra, we apply the minimisation to a compressed data set. 
For each light curve of the SALT2 training, a flux scale is fitted 
to improve the match to the best fit model. An effective wave- 
length of the corresponding rest-frame filter is also computed 
(accounting for the spectrum of the SN in the filter integration 
window). Those flux scales with their associated uncertainties 
and the effective wavelength provide us the data needed for an 
estimate of the broad-band colour uncertainties. In this process, 
we fit simultaneously corrections to the colour variation law and 
the colours of the average SED. 

This technique has been successfully tested on Monte Carlo 
simulations. The corrections to the SED and colour variation law 
turn out to be negligible at the end of the minimisation: about one 
milli-magnitude for the colours of the SED and 4% of the colour 
variation amplitude in the [/-band (or 0.005 mag difference in 
[/-band for C = 0.1), with the caveat that the [/-band dispersion 
(Tu has been accounted for in the training of SALT2F1. 

The remaining scatter as a function of wavelength for the two 
parameterizations of Eq. IA.2l is shown in Fig. [8] We also obtain 
o-u = 0.11+0.02. 

A3. Propagation of the model statistical uncertainties 

As discussed in §4.11 the statistical uncertainties of the SED se- 
quence model have to be propagated to the estimates of luminos- 
ity, shape and colour of the SNe. Indeed, they introduce redshift- 
dependent uncertainties on distances that do not cancel with an 
infinite number of SNe. 

For SALT2, we use the covariance matrix of the training, 
including the residual scatter evaluated in the previous section. 
This latter is accounted for by adding to the covariance matrix of 
each light curve a term k(A)MM T , where M is a vector contain- 
ing the model values corresponding to each photometric point, 
and A is the effective rest-frame wavelength of the filter used for 
the observations. Let 9 be the parameters of the model, Cg their 



17 The [/-band calibration uncertainty is accounted for with an ad- 
ditional term crjjMM 7 in the covariance matrices of the light curves, 
where M is a vector of the model values for each data point. 
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covariance matrix, and 77, the parameters (m* B ,Xi,C) of a given 
SN, we may write the^ 2 as: 

* 2 = Z [fi ~ m{e > m)]T Wi [fi ~ mi{9 > " )] (A - 3) 

i=l,iV 

where 

- fi are the flux measurements of SN i 

- nij is the model. 

- Wi incorporates both the measurement uncertainties and the 
error model discussed previously. 

For small deviations (66, 6rjj) to a set of parameters (do, ^o), the 
normal equations provide us with the following relation between 
drjj and 66: 

rji = (BiWiBjy 1 [BiWiffi - A] 66)] (A.4) 

where r, = fi - 111,(60, ?7i,oX A, = dm-JctO and B, = dmi/dr//. This 
can be inserted into the "global" normal equation: 

J] AiWiAj - DjE; l Di 66= [ A ' W < ~ ^Wi] n 

i=l,N J i=lJSt 

(A.5) 

with: 

D t = BWAj 
E t = BiWiBj 

This defines 56, and the parameters £77, are obtained from Eq. 
IA.4I The covariance matrix blocks can be computed: 



Cov(6,6T l = Cg 1 = J] A^Aj -DjEi'Di 

i=l,N 

Cov(6,j li ) = -CgDjE; 1 
Cov^jjj) = dijE; 1 +Er l D i C g D^E] 1 



(A.6) 



In the expression JA.6b for Cov(r]i, rjj), the first term accounts for 
the photometric noise, the second accounts for the model noise. 
Note that when considering Cov( 77;, 77,), the second term does not 
vanish when the photometric data of this SN becomes infinitely 
precise, i.e. W, — > 00. 

These expressions are valid for SNe in the training sample. 
We can as well compute the contribution of model uncertainty to 
the fit of a single SN which is not in the training sample. The^- 2 
is the same as above, but 6 is fixed (and 66 = 0). Using the same 
notations as above, the normal equation for 6rji now reads: 

BiWiBjdTji = BfWiTi 

and the uncertainty of the model affects the 77,- estimator 77,- by 
offsetting r,\ 



drji df\i dr- { 
dd dr-, dd 



[(BiWiBjy^iWi] [aJ] = e; 1 d, 



So the covariance of SN parameters due to model uncertainties 
reads: 

Cov 9 (j li ,j lj ) = E7 1 D i CgD T j Er l 



which is identical to the second term in expression JA.6b above. 
So the expression of the covariance of SN parameters' estimators 
does not depend on their belonging to the training sample. Note 



however that when the training sample increases, this covariance 
decreases, as we expect. 

We can now compute a covariance matrix of distance mod- 
uli (jii) of the whole sample due to light curve model statisti- 
cal uncertainties. This covariance matrix is added to other noise 
sources and used in the cosmological fit. It reads: 

C^UJ) = Cov(Mi,fij) = V T EJ y D i C e D' j E j V 

where V is defined by fi, = V r 77, (V = {I, a, -fi} from Eq.[5]). 

The uncertainty of the difference between the average dis- 
tance modulus in a redshift bin [21,22] an d the average on the 
whole sample (see Fig.fToTi is given by 



1/2 



>r: . .. -(V' :: (),V. : ) 



where V ZlZ2 (i) = (l/N ZlZ2 - 1 /N) if z\ < Zi < Zi, and -l/N if not, 
z, being the redshift of the SN number i, N ZlZt the number of SNe 
in the bin, and N the total number of supernovae in the sample. 

A.4. Public release of the light curve fitter 

The SALT2 model parameters used in this paper, the code 
to fit light curves, the calibration files and models of 
the instrument responses are public. They can be found 
at |http : / /supernovae . in2p3 . £r/~guy/salt/ along with 
some documentation. Note however that the usage of this pack- 
age requires some care, in particular the rest-frame wavelength 
validity range of the model is limited. Also, the photometric cal- 
ibration of the SNe light curves of the training sample is im- 
printed in the spectral model broadband colours. Hence it is bet- 
ter suited to retrain the model when more precise calibrations of 
those samples become available, or when new data sets that may 
improve the model can be used. SALT2 is indeed more a method 
than a final spectral model. Finally, as described in detail in this 
paper, the uncertainties of the model parameters have to be taken 
into account (training sample statistics, calibration uncertainties, 
residual scatter model), and the codes to propagate those are not 
easy to distribute. For all these reasons, users are welcome to 
discuss with the author of their usage of the light curve fitter. 



Appendix B: Flux bias on PSF photometry 

We consider the standard PSF photometry when flux and posi- 
tion are fitted simultaneously, using a PSF P(x, y), normalised to 
1. We call the flux and position parameters (/, 6 X , 6 y ). Assuming 
the noise is stationary and the object faint, the weight matrix of 
their estimators reads: 



W = w 



(jP 2 fj Pd x P f f PdyP 

fj(d x P) 2 f 2 Jd x Pd x P 
f 2 j(d y P) 2 



where w is the inverse of the noise variance per unit area. 

Assuming that the PSF is symmetric in x and y, the non- 
diagonal terms vanish, and the flux and position estimates are 
thus independent. The sums over pixels have been replaced by 
integrals which is adequate for our well sampled images. The 
flux estimator reads: 



J I(x, y)P(x - 8 X , y - 6y)dxdy 
J P 2 (x,y)dxdy 
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where both sums run over pixels, the true position is assumed to 
be (0,0), the fitted position is (d x , 5 V ), and I is the sky-subtracted 
image. Assuming that the image reads: 

I(x,y) = fP(x,y) + e(x,y) 

(with e representing the measurement noise: E[e] = 0, and its 
variance reads per unit area), the flux estimator expectation 
value can be written as: 



E[f] = f 



E[ JP(x,y)P(x - 6 x ,y - 5 y )dxdy] 
J P 2 (x,y)dxdy 

1 £[<5/] / Pd\P + E\5y\ j Pd 2 y P 
J P 2 (x, y)dxdy 




+ 0{5 A )\ 



The flux and position variances read: 



Var[f] = |w J P 2 J 

Var[5A = E[5l] = [wf j (5,P) 2 J = (-wf 2 j Pd 2 p) 



Hence, 



E[f\*fU 



VarU} 

f 2 



So, the PSF flux is underestimated at low S/N, due to the inac- 
curacy of the position, because position shifts yield on average 
smaller PSF fluxes. Note that expression IB. 1 1 assumes that flux 
and position are measured from the same data, and also assumes 
a stationary noise and a faint object. The neglected higher order 
terms may become important below S/N of a few. For S/N above 
1, the bias is smaller than the uncertainty, but should anyway be 
studied if averages of faint fluxes are used in the analysis, such 
as at the high redshift end of a Hubble diagram. 
For a Gaussian PSF, we have: 

Var[f] = w^Ancr 2 
Var\5 x \ = w-^JTtr 4 // 2 

where <x stands for the Gaussian RMS This, combined with 
Eq. IB.ll gives 



where the sum runs over images. The amplitude of the light 
curve will be fitted to the individual image measurements fi, and 
is a linear combination of those. The least squares estimator of 
the amplitude A of the light curve reads: 



A — k ^ a,fj, with a; 



fi 



Var[fi] 



where k is a global constant depending on the chosen definition 
of A. We now propagate the flux bias at individual epochs into 
the lightcurve amplitude: 




Using that: 



(B.l) We finally have: 



Var[A] 



E[A] = A 



Z 



1 - 



ft 



Var[fi] 



Var[A] 
~A 1 ~ 



and we can hence use the relative uncertainty of the lightcurve 
amplitude as an indicator of its own bias, when position and am- 
plitude are measured from the same data. 



E[f]=f\l~ 



E[5 X ' + 8j] 
4cr 2 



(B.2) 



Schematically, the measurement of a supernova light curve 
consists of a fit of a common position and varying fluxes to an 
image series. As for a single image, the flux estimators are bi- 
ased: 

- ( £[(^ 2 +(£ 2 ]] 

*M=/,|i--4^-j 

where f, is the supernovae flux at epoch ; and cr, is the RMS of 
the PSF (assumed Gaussian) at the same epoch. One may note 
that the images with the poorest seeing are the less affected. We 
assume the position to be measured from the images series. Its 
variance reads: 



E[6 X 2 + c5 v 2 ] 



2Var[6 x ] 



z 



Jf 



\o- 2 Var\ft) 
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